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Abstract 

We review the method of symplectic invariants recently introduced to solve matrix 
models loop equations, and further extended beyond the context of matrix models. For 
any given spectral curve, one defines a sequence of differential forms, and a sequence of 
complex numbers Fg. We recall the definition of the invariants Fg, and we explain their 
main properties, in particular symplectic invariance, integrability, modularity,... Then, 
we give several examples of applications, in particular matrix models, enumeration of 
discrete surfaces (maps), algebraic geometry and topological strings, non-intersecting 
brownian motions,... 
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1 Introduction 

Recently, it was understood how to solve, order by order in the so-called "topological 
expansion", the loop equations (Schwinger- Dyson equations) for matrix integrals [17]. 
The solution brought an unexpectedly rich structure [60], which, did not only solve 
the 1-matrix model, but which also solved multi-matrix models, as well as their limits. 
Later, it was understood that this structure also appears in other matrix models, and 
in problems of enumerative geometry, not directly related to matrix models. 

Thus, there is an underlying structure which can be defined beyond the context of 
matrix models, and relies only on the intrinsic algebro-geometric properties of a plane 
curve, called the spectral curve. 

In other words, for any regular (to be defined below) complex plane curve £ = 
{y{x)} (whether it is related to a matrix model or not), we can define a sequence 
of numbers Fg{S), g = 0, 1, 2, ... , oo. Those numbers Fg are called the symplectic 
invariants of the spectral curve S (first introduced in ^60j). The reason is because 
two spectral curves S and S which can be deduced from one another by a symplectic 
transformation (i.e. they have the same wedge product dx A dy = dx A dy), have the 
same Fg's. Fg is called the symplectic invariant of degree 2 — 2g , because under a 
rescaling y Xy, Fg scales as Fg — > \^~'^^Fg (except Fi which is logarithmic). 

Moreover, for a spectral curve £ = {y{x)}, we define not only its symplectic invari- 
ants -Fg's, we also define a doubly infinite sequence of symmetric meromorphic forms 
LJn\xi, . . . , x„)[£^], G N, (? G N, and such that Fg = uj\^\ For n > 1, those forms are 
not symplectic invariants, but they have many nice properties. They allow to compute 
the derivatives of the F^'s with respect to any parameter on which £ could depend. 

Those geometric objects are interesting, in particular for their applications to vari- 
ous problems of enumerative geometry (each problem corresponding to a given spectral 
curve £), but also on their own. Indeed they have remarkable properties for arbitrary 
spectral curves, i.e. even for spectral curves not known to correspond to any enumer- 
ative geometry problem. 

In particular, they are related to the Kodaira-Spencer field theory, to Frobenius 
manifolds, to the WDVV special geometry, topological strings and Dijkgraaf-Vafa con- 
jecture. They are expected to be the B-model partition function, and through mirror 
symmetry, the F^'s are thus expected to be the generating functions of Gromov-Witten 
invariants of genus g for some toric geometries. 

As another example of interesting properties, the F^'s have nice modular behaviors, 
and, for instance, they provide a solution to holomorphic anomaly equations. 

They also contain an integrable structure, related to "multicomponent KP" hierar- 
chy, e.g. they satisfy determinantal formulae, Hirota equations,... 
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Another nice property, is that they can be computed by a simple diagrammatic 
method, which makes them really easy to use. For instance the holomorphic anomaly 
equations can be proved only by drawing diagrams. 

Regarding the applications, we will consider the following examples: 

- Enumeration of discrete surfaces, possibly carrying colors on their faces (Ising 
model), as well as the asymptotics of large discrete surfaces. 

- For the curve y = ^ sin {27iy/x), the F^'s compute the Weyl-Petersson volumes. 

- We will consider also the Kontsevich spectral curve, related to the Kontsevich 
integral, for which the F^'s are generating functions for intersection numbers of Chern 
classes of cotangent bundles at marked points, and the wjf^^s are generating functions 
of Mumford k classes (in some sense by "forgetting" some marked points). 

- For the curve y = Argcosh(x), and deformations of that curve, the F^'s are 
generating functions for counting partitions with the Plancherel measure, related to 
the computation of Hurwitz numbers. 

- g-deformed versions of Plancherel measure sums of partitions can also be com- 
puted with symplectic invariants of some appropriate spectral curve, which, not so 
surprisingly, is the (singular locus of the) mirror of a toric Calabi-Yau manifold. This 
is consistent with the conjecture that the F^'s are related to Gromov-Witten invariants. 
Indeed, Gromov-Witten invariants of toric Calabi-Yau 3-folds can be computed, using 
the topological vertex, as sums of partitions, typically q-deformed Plancherel sums, for 
the simplest examples of toric Calabi-Yau 3-folds. 

2 Symplectic invariants of spectral curves 

Symplectic invariants were introduced in [60], as a common framework for the solution 
of loop equations of several matrix models: 1-matrix, 2-matrix, matrix with external 
fields,..., as well as their double scaling limits. Then it was discovered that they have 
many nice properties, in particular symplectic invariance, and that they appear in other 
problems of enumerative geometry, not necessarily related to random matrices. 

Here we only briefly summarize the construction of [00], without proofs, and we 
refer the reader to the original article for more details. 

2.1 Spectral curves 

In this article, we define a spectral curve as follow^]: 

Definition 2.1 A spectral curve £ = {C,x,y), is the data of a compact Riemann 
surface C, and two analytical functions x and y on some open domain in L. 

In some sense, we consider a parametric representation of the spectral curve y[x\ 
where the space of the parameter 2; is a Riemann surface L. 

"'^This definition is not exactly the one usually encountered in integrable systems [14], in fact it 
turns out that the plane curve we are considering here, is the "classical limit" of the full spectral 
curve. We call it spectral curve by abuse of language, and because it has become customary to do so. 
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Definition 2.2 If C is a compact Riemann surface of genus g, and x and y are mero- 
morphic functions on C, we say that the spectral curve is algebraic. If in addition, C is 
the Riemann sphere (C = = C U {00} , i.e. of genus g = 0), we say that the spectral 
curve is rational. 

Indeed, for an algebraic spectral curve, it is always possible to find a polynomial 
relationship between x and y: 

Pol(x,?/)=0. (2-1) 

For a rational spectral curve, the polynomial equation Pol{x,y) = 0, can be parame- 
terized with two rational functions x{z) and y{z) of a complex variable z. 

Definition 2.3 A spectral curve {C,x,y) is called regular if: 

• The differential form dx has a finite number of zeroes dx{ai) = 0, and all zeroes 
of dx are simple zeroes. 

• The differential dy does not vanish at the zeroes of dx, i.e. dy{ai) 7^ 0. 

This means that near x{ai), y behaves locally like a square-root y{z) ~ y(cti) + 
C \/x{z) — x{ai), or in other words, that the curve y{x) has a vertical tangent at ai. 

From now on, we assume that we are considering only regular spectral curves. 
Symplectic invariants are defined only for regular spectral curves, and they diverge 
when the spectral curve becomes singular. Examples of singular spectral curves are 
considered in section 14.81 they play a central role in the double scaling limit in chapter 

El 

Definition 2.4 We say that two spectral curves £ = {C,x,y) and S = {C,x,y) are 
symplectically equivalent if there is a conformal mapping C ^ C, and if under this 
mapping dx A dy ^ dx A dy. The group of symplectomorphisms is generated by: 

• X = X, y = y + R{x), R{x) = rational function of x. 

cx+d' ^ ad— be 

• X = f{x), y = jTf^y, where f is analytical and injective in the image of x. 

• X = y, y = —x. 

All those transformations conserve the symplectic form on C, whence the name: 

dx A dy = dx A dy. (2-2) 

The main property of the F^'s we are going to define, is that they are symplectic 
invariants, i.e. two curves which are symplecticaly equivalent, have the same F^'s. 
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2.1.1 Examples of spectral curves 

Interesting examples of spectral curves may come from several areas of physics or 
mathematics, and are related to some problems of enumerative geometry. We will 
study in details some examples between section [5] and [TH Here, in order to illustrate 
our notion of spectral curve, we give some examples of spectral curves of interest 
extracted from those applications. 

For the readers familiar with matrix models, the spectral curve under considera- 
tion here, can be thought of, as the "equilibrium density of eigenvalues of the random 
matrix". It is not to be confused with the large N density of eigenvalues, although, 
for many simple cases the two may coincid^. In the most simple matrix models, the 
spectral curve is algebraic. For formal random matrix models, designed as combina- 
torics generating funcions for counting discrete surfaces, the spectral curve is shown to 
be rational (see section [7]). 

In the context of string theory, the spectral curve is often given by a transcendental 
equation of the form if(e'^',e^) = 0, where if is a polynomial. It is not an algebraic 
spectral curve, but is closely related to an algebraic curve. In that case, dx and dy are 
abelian meromorphic differentials on the compact Riemann surface C corresponding to 
H (see section [TT]). 

The origin of all the examples below are described with more details in sections \5\ 

• The following curve is a rational spectral curve: 

£ = pi = C U {oo} , x{z) = , y{z) = - 3^. (2-3) 

It satisfies the algebraic equation y"^ — 2 = — 3x. It is an hyperelliptical curve of 
genus ^ = 0. This spectral curve is related to the so-called "pure gravity Liouville field 
theory". It will often be called the "pure gravity" spectral curve, or also the (3,2) 
spectral curve, because pure gravity is the (3, 2) minimal conformal field theory, it has 
central charge c = 0. See sections 14.81 and [HI 

• The curve y = y/x, is also a rational spectral curve which satisfies y'^ = x, and 
which can be parameterized by: 

£ = , x{z) = z^ , y[z) = z. (2-4) 

This spectral curve arises in the study of the extreme eigenvalues statistics of a random 
matrix, i.e. in the study of the Tracy- Widom law and of the Airy kernel [111] . It will 
often be called the "Airy" spectral curve. It will also be called the (1, 2) spectral curve 
in order to match the classification of minimal conformal field theories. The minimal 
model (1,2) has central charge c = —2. See section [HI 

• The following spectral curve is also a rational spectral curve: 

C = F' , x{z)=j(z + -) , y(z) = -± + ^-^ (2-5) 

\ Z J '-)Z z-^ 



^The two notions coincide for example for matrix integrals with a polynomial potential. They do 
not coincide for example when the potential has an explicit dependence on TV. 
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where 7^ = . This spectral curve arises in the enumeration of quadrangulated 

surfaces, i.e. in the formal quartic matrix model. See section [731 



The following spectral curve 

£ = , x{z) 



y{z) = — sin {2ttz) 
Zn 



(2-6) 



is related to the computation of Weyl-Petersson volumes. Notice that it is not algebraic, 
but it can be parameterized by a complex variable, i.e. by a ^ = Riemann surface. 
See section [T0.2[ 



The following rational spectral curve 



- 3tz 



[2-7) 



is singular at t = 0. Indeed at t = 0, the differential dy = 8(2;^ — t)dz vanishes at 2; = 
which is the zero of dx. We will see below that the F^'s diverge for singular curves, 
and thus the function Fg{t) has a singularity at t = 0. Just from homogeneity, and by 
considering the change of variable z ^ \/t z, v^e see that: 

F,(t)=t^(i-^)F,(1) (2-8) 

which indeed diverges at t = 0. See section 14.81 

• The following spectral curve depends on two parameters p G Z and G C*: 



£ = P^ 



x{z) 
viz) 



(1_ 



1 ^2 



(2-9) 



x(z) 



In z + - In 
2 



zjzQ 



1 - Ijzz^ 



It appears in the enumeration of g— deformed Plancherel sums of partitions, i.e. in 
the computation of the Gromov-Witten invariants of the toric Calabi-Yau manifold 
Xp = 0{-v) © 0{v - 2) ^ P^ See section ESI 

This spectral curve is symplectically equivalent to (compute dx f\dy in both cases): 



£ = P^ 



x{z) =ln( (\--){\- — 

Zq ZZq 




z/zo 



(2-10) 



1 - 1/zzq^ 



This last spectral curve is such that e^ and e^ are rational functions of 2, and thus by 
eliminating z, there exists a polynomial if(e^,e^) such that: 



i/(e^,e^) = 0. 



(2-11) 
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This equation is precisely the singular locus of the mirror manifold of Xp. The full 
mirror manifold (not only its singular locus) is the 3 dimensional submanifold of 
locally given by {(x, 0;+, u;_) G C^/ if(e^,e^) = cu+cu-}. See section [TTl 

• The following spectral curve is of genus (7 = 1, it is algebraic but not rational: 

/: = C/(Z + rZ) , x{z) = p{z,t) , y{z) = p\z,T) (2-12) 

where p is the Weierstrass function, and £ is the torus of modulus r. It is algebraic 
because the Weierstrass function obeys the differential equation: 

= 4p3 _ g^^ - ^3. (2-13) 

This spectral curve is called the Seiberg-Witten curve since it first appeared in a 
solution to M = 2 Supersymmetric Yang-Mills theory proposed by Seiberg and Witten 
in [TUT] . 

2.2 Geometry of the spectral curve 
2.2.1 Genus and cycles 

The only compact Riemann surface of genus (7 = is the Riemann sphere = CU{oo}. 
It is simply connected. 

A compact Riemann surface C of genus g can be equipped with a symplectic 
basis (not unique) of 2g non-contractible cycles such that: 

n Bj = 6ij , AnAj = , BiH Bj = 0. (2-14) 

They are such that C\ {UiAi Uj Bi) is a simply connected domain of C, which we shall 
call the fundamental domain. 




The choice of cycles and of a fundamental domain is rather arbitrary, and is not 
unique. Many of the quantities we are going to consider depend on that choice. 
The quantities which do not depend on that choice are called modular invariant. 

On a compact Riemann surface of genus ^ > 1, there exist holomorphic differential 
forms (analytical everywhere on C, in particular with no pole). Those holomorphic 
forms clearly form a vector space (linear combinations are also holomorphic) over C, 
and this vector space has dimension g. 
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When we have a choice of cycles Ai,Bj, it is possible to choose a basis (which is 
unique), which we call dui, . . . ,dug, and normalized such that: 

® duj = 6ij. (2-15) 

Once we have defined those du/s, we can compute the following Riemann matrix of 
periods: 




(2-16) 



This matrix r^j is symmetric, and its imaginary part is positive definite: 

Tij=Tj,i , lmr>0. (2-17) 

2.2.2 Abel map 

Consider an arbitrary origin o in the fundamental domain, and fixed throughout all 
this article. 

For any point z in the fundamental domain, the vector {ui{z), . . . , Ug{z)): 




(2-18) 



where the integration path is in the fundamental domain, is called the Abel map of 
z. It is a vector in C^. It depends on the choice of o by an additive constant, and it 
depends on the choice of the fundamental domain, by a vector in the lattice + r1i^. 
The quotient C^/(Z^ + rZ^) is called the Jacobian. 

The Abel map, sends points of C to points in the Jacobian. 

2.2.3 Bergmann kernel 

Given a choice of cycles, we define the Bergmann kernel: 

B{z^,Z2) (2-19) 

as the unique bilinear differential having one double pole at Zi — z^ (it is called "2nd 
kind") and no other pole, and such that, in any local parameter z: 

B{zuZ2) - -^^^i^ + reg , Vi = l,...,^, (f B{zi, Z2) = 0. (2-20) 

One should keep in mind that the Bergmann kernel depends only on C, and not on 
the functions x and y. 

The Bermann kernel can be seen as the derivative of the Green function, i.e. the 
solution of the heat kernel equation on C 

The Bergmann kernel is clearly unique because the difference of two Bergmann 
kernels would have no pole, and vanishing .4-cycle integrals, therefore it would vanish. 
It is also interesting to note that it is symmetric in its variables zi and Z2- 

Examples: 
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• if £ = = C U {00} =the Riemann Sphere, the Bergmann kernel is a rational 
expression: 

Bi'-'^) = J^r (2-21) 

Most of the applications between section [5] and section [TT|, will be on £ = P^, 
and will use this rational Bergmann kernel. 

• if >C = C/ (Z + rZ) =Torus of modulus r, the Bergmann kernel is 

Z2) = [p{zi - Z2, r) + ^) dzidz2 (2-22) 

where p is the Weierstrass elliptical function. 

• if £ is a compact Riemann surface of genus g > 1, of Riemann matrix of periods 
Tij, the Bergmann kernel is 

Bizi, Z2) = d^^d^^ In i9{u{zi) - u{z2) - c, r)) (2-23) 

where u{z) is the Abel map, c is an odd characteristic, and 6 is the Riemann 
theta function of genus g (cf [66l[65] for theta-functions). 

2.2.4 Generalized Bergmann kernel 

Given an arbitrary symmetric matrix k of size g x g, we consider a "deformed" 
Bergmann kernel: 

9 

B^{zi, Z2) = B{zi, Z2) + 2z7r ^ dui{zi)duj{z2). (2-24) 

If K = we recover the usual Bergmann kernel Bq = B. 

The reason for introducing this k, is that a change of basis of cycles and fundamental 
domain, can be rewritten change of k. 

Indeed, perform a Sp2g{'^) change of symplectic basis of cycles {C,D,C,D have 
coefficients in Z and CD* = DC\ C& = DC\ C& - D& = 1): 

A = E^M^; + E^M^;- ' ^. = Ec'm-4; + $]a,,s;. (2-25) 
j 3 3 3 

The Riemann matrix of periods r' in the new basis A! , B', is related to the old one by 
the modular transformation: 

t' = {D - tD)-\tC - C) (2-26) 

and the Bergmann kernel changes as: 

9 

Bo ^ B'q = Bo + 2i7c ^ dui{zi)duj{z2) , k = {DD^^ - t)'^ (2-27) 
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in other words, the change of cycles can be reabsorbed as a change of k. 

More generally, the kernel in a basis A, B is equal to B'^, in the basis A',B', 
where: 

k' = {& - D^r)K{D - tD) - (D* - D^t)D. (2-28) 

From now on, we will always consider i?^, and we will write B instead of i-e. 
we will omit the n subscript, unless ambiguity. However, for most of the practical 
applications, one often chooses k = 0. 

2.2.5 SchifFer kernel 

In particular if we choose k to be the Zamolodchikov Kahler metric: 

K={T-T)-'='-{\mrr\ (2-29) 
we see that in the new basis A! the matrix k becomes 

k' (Imrr\ (2-30) 

i.e. it takes the same form as in the initial basis. Therefore, with this special value of 
K, the Bergmann kernel -B^ is called the Schiffer kernel [16] and it is modular invariant: 
it does not depend on a choice of cycles. However, the price to pay to have modular 
invariance, is to have a non analytical dependence in r, and thus a non analytical de- 
pendence in the spectral curve. This incompatibility between analyticity and modular 
invariance is the origin of the so-called "holomorphic anomaly equation", see section 
0:21 

Example: if £ = C/(Z + rZ) =Torus of modulus r, the Schiffer kernel is 

B{zi, Z2) = p{zi - Z2, r) dzidz2 (2-31) 
where p is the Weierstrass elliptical function. Compare with eq. fl2-22p . 

2.2.6 Branchpoints 

Branchpoints are the points with a vertical tangent, they are the zeroes of dx. Let 
us write them Oj, i = 1, . . . , #bp: 

Vi, dx{ai) = 0. (2-32) 

Since we consider a regular spectral curve, all branchpoints are simple zeroes of dx, 
the curve y{x) behaves locally like a square root y{z) ~ y{ai) + Ci^/x{z) — x{ai), near 
a branchpoint Oj, and thus, for any z close to a^, there is exactly one point z ^ z in 
the vicinity of such that: 

x{z) = x{z). (2-33) 

z is called the conjugated point of z. It is defined locally near each branchpoint Oj, 
and it is not necessarily defined globally. 

Examples: 



12 



enumeration of maps ^ 1-matrix model in the 1-cut case: 

in this case we have L =Riemann sphere, and (see section [?.3.2p : 

= a + 7(2; + 1/2;) , dx[^z) = X [^z)dz = '^{\ — z~'^) dz. (2-34) 
The zeroes of dxi^z) are z = ±1, and we clearly have z = 1/z: 

ai = 1, 02 = — 1 , z = 1/ z. (2-35) 
In this case z is defined globally, 
pure gravity (3,2): 

in that case we have C =Riemann sphere, and 

x{z) = z^ -2 , dx{z) = 2zdz. (2-36) 
The only zeroe of dx{z) is z = 0, and we have z = —z: 

a = , z = -z. (2-37) 
In this case z is defined globally. 
Ising model (4,3): 

in that case we have C =Riemann sphere, and 

x{z) = z^ -3z , dx{z) = 3(^2 - 1) dz. (2-38) 

The zeroes of dx{z) are = ±1, and near = ±1 we have: 

1 



ai = ±l , z = --{z-aiVl2-3z^). (2-39) 
In this case z is not defined globally, and it depends on a^. 
2.2.7 Recursion Kernel 

For any Zq G C, and any z close to a branchpoint, we define the recursion kernel: 



K{zo,z) 



2 iyiz)-yiz))dxiz) 



(2-40) 

where the integral is taken in a small domain in the vicinity of the concerned branch- 
point. 

K{zq,z) is a meromorphic 1-form in the variable zq, it is defined globally for all 
zq G £, it has simple poles at zq = z and zq = z. 

On the contrary, in the variable z, the kernel K{zq, z) is defined only locally near 
branchpoints z ~ a^, and it is the inverse of a differential. As we shall see below, 
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K{zo, z) will always be used only in the vicinity of branchpoints, and it will always be 
multiplied by a quadratic differential in z, so that the product will be a differential 
form. 

Let us notice that K{zq, z) = K{zo, z), and that K{zo, z) has a simple pole when z 
approaches the branchpoint. Using De L'Hopital's rule, the leading behavior near the 
branchpoint is: 

K{zo, z) ~ . + regular . . . (2-41) 

2.3 Correlation functions 

We start by defining a sequence of meromorphic n-forms Un^ with n = 1,2, .. . and 
g = 0,1,2, .. ., called correlators or correlation functions, by the following recursion: 

Definition 2.5 Given a spectral curve £ = {C,x,y), and a matrix k (see section 
2.2.4\), we define recursively the following meromorphic forms: 



u;f\z) = -y{z)dx{z) (2-42) 

J^\zi,Z2) = Bizi,Z2) (2-43) 
and if 2g — 2 -\- n > , and J is a collective notation for n variables J = {zi, . . . , Zn}: 



9 I 

(9) 



h=o icJ 



(2-44) 

where in the RHS means that we exclude the terms with [h, I) = (0, 0), and {g, J). 

This definition is indeed a recursive one, because all the terms in the RHS have a 
strictly smaller 2g — 2 + n than the LHS. 

The functions Un'^ with 2 — 2g — n < are called stable, the others are unstable (the 
only unstable ones are thus uf*^ and ). 

ujn\zi, . . . , Zn) is a mcromorphic 1-form on C in each variables Zi. It can be proved 
by recursion, that it is in fact a symmetric form. Moreover, if 2 — 2g — n < 0, its only 
poles are at branchpoints Zi — > aj, and have no residues: 

Res 4^)(^i,^2,...,^n) = 0. (2-45) 
Those properties can be proved by recursion, and we refer the reader to |60j . 
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2.4 Free energies 

The previous definition, defines Un'' only if > 1. Now, we define Fg = uj^\ called 
"Free energies" or "symplectic invariant of degree 2 — 2(7" , by the following: 



Definition 2.6 Symplectic invariants 

We define for g > 2: 




(2-46) 



where $ is any function defined locally near branchpoints, such that = ydx ((^ is 
defined up to an additive constant, hut thanks to eq. lli2-45\ ), Fg does not depend on the 
choice of that constant). 

The unstable cases g = and g = I are special. We have: 



Definition 2.7 For g = 1 we define 



F, = ^\n(^TB{{xia,)}) Yly'ia,)^ 



where we define: 



y'iai) 



lim 



y{z) - y{ai 



(2-47) 



(2-48) 



X a,- 



and tb is the Bergmann r-function of Kokotov-Korotkin |g7| /. Ifx{z) is a meromorphic 
function on C, tb depends only on the values of x at its branch points, i.e. Xi = x{ai). 
It is defined by: 



d\nTB{{X,)}) 



dXi z^a, dx{z) 

Notice that B here stands for 5^ with arbitrary k. 



B z,z) 



(2-49) 



The definition of Fq is more involved, and we refer the reader to [60]. A convenient 
way to define Fq, is through its 3rd derivatives, using theorem 14.31 below. In fact, all 
the Fg's with g > 1 are obtained in terms of local behaviors around branchpoints, but 
Fq depends on the whole spectral curve, not only on the vicinity of branchpoints. In 
the context of topological strings, Fq is called the prepotential. 
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3 Diagrammatic representation 

The recursive definitions of u^'' and F^^' can be represented graphically. 

We represent the A;— form u}l^\pi, . . . ,Pk) as a "blob-hke surface" with g holes and 
k legs (or punctures) labeled with the variables Pi, ■ ■ ■ ,Pk, and F'^^^ — ulf^ with legs 
and g holes. 




We represent the Bergmann kernel B{p,q) (which is also , i.e. a blob with 2 
legs and no hole) as a straight non-oriented line between p and q 

B{p,q) := p q . (3-2) 

We represent K{p, q) as a straight arrowed line with the arrow from p towards q, 
and with a tri-valent vertex whose left leg is q and right leg is q 

K{p,q):= P ► . (3-3) 

Graphs 

Definition 3.1 For any k > and g > such that k + 2g > 3, we define: 

Let Gjf^i{p,Pi, ■ ■ ■ ,Pk) be the set of connected trivalent graphs defined as follows: 

1. there are 2g + k — 1 trivalent vertices called vertices. 

2. there is one 1-valent vertex labelled by p, called the root. 

3. there are k 1-valent vertices labelled with pi, . . . ,pk called the leaves. 

4. There are 3g + 2k — 1 edges. 

5. Edges can be arrowed or non-arrowed. There are k -\- g non-arrowed edges and 
2g -\- k — 1 arrowed edges. 

6. The edge starting at p has an arrow leaving from the root p. 

7. The k edges ending at the leaves pi, ■ ■ ■ ,Pk are non-arrowed. 
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8. The arrowed edges form a " spanning plana% binary skeleto"^ tree" with root p. 
The arrows are oriented from root towards leaves. In particular, this induces a 
partial ordering of all vertices. 

9. There are k non-arrowed edges going from a vertex to a leaf, and g non arrowed 
edges joining two inner vertices. Two inner vertices can be connected by a non 
arrowed edge only if one is the parent of the other following the arrows along the 
tree. 

10. If an arrowed edge and a non-arrowed inner edge come out of a vertex, then the 
arrowed edge is the left child. This rule only applies when the non-arrowed edge 
links this vertex to one of its descendants (not one of its parents). 



Example of Ql ' (p) 

As an example, let us build step by step all the graphs of Q[^\p), i.e. g = 2 and 
k = 0. 

We first draw all planar binary skeleton trees with one root p and 2g -\- k — 1 = 3 
arrowed edges: 




Then, we draw g + k = 2 non-arrowed edges in all possible ways such that every vertex 
is trivalent, also satisfying rule 9) of definition 13. 1[ There is only one possibility for the 
first tree, and two for the second one: 




It just remains to specify the left and right children for each vertex. The only 
possibilities in accordance with rule 10) of def J3.1l arj^: 




^It goes through all vertices. 

^Planar tree means that the left child and right child are not equivalent. The right child is marked 
by a black disk on the outgoing edge. 

binary skeleton tree is a binary tree from which we have removed the leaves, i.e. a tree with 
vertices of valence 1, 2 or 3. 

^ Note that the graphs are not necessarily planar. 
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(3-6) 

In order to simplify the drawing, we can draw a black dot to specify the right child. 
This way one gets only planar graphs: 




(3-7) 

Remark that without the prescriptions 9) and 10), one would get 13 different graphs 
whereas we only have 5. 



Weight of a graph 

Consider a graph G G ^^!|li(p,pi, . . . ,Pfc)- Then, to each vertex i — 1, . . . ,2g + k — 1 
of G, we associate a label qi e £, and we associate qi to the beginning of the left child 
edge, and to the right child edge. Thus, each edge (arrowed or not), hnks two labels 
which are points on the spectral curve C 

• To an arrowed edge going from q' towards g, we associate a factor K{q',q). 

• To a non arrowed edge going between q' and q we associate a factor -B(g', q). 

• Following the arrows backwards (i.e. from leaves to root), for each vertex g, we 
take the sum over all branchpoints of residues at g — >■ a^. 



After taking all the residues, we get the weight of the graph: 

w{G) 

which is a multilinear form in p, pi , . . . , . 

Similarly, we define weights of linear combinations of graphs by: 

w{aGi + I3G2) = aw{Gi) + l3w{G2) 

and for a disconnected graph, i.e. a product of two graphs: 

w{G^G2) = w{Gi)w{G2). 

Theorem 3.1 We have: 

( 



(3-8) 



(3-9) 



(3-10) 



E 



G 



(3-11) 
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proof: 

This is precisely what the recursion equations 12-441 of def l2.5l are doing. Indeed, one 
can represent them diagrammatically by 




(3-12) 



□ 

Such graphical notations are very convenient, and are a good support for intuition 
and even help proving some relationships. It was immediately noticed after [57] that 
those diagrams look very much like Feynman graphs, and there was a hope that they 
could be the Feynman's graphs for the Kodaira-Spencer quantum field theory. But they 
ARE NOT Feynman graphs, because Feynman graphs can't have non-local restrictions 
like the fact that non oriented lines can join only a vertex and one of its descendent. 

Those graphs are merely a notation for the recursive definition 12-441 

Lemma 3.1 Symmetry factor: 

The weight of two graphs differing by the exchange of the right and left children of 
a vertex are the same. Indeed, the distinction between right and left child is just a way 
of encoding symmetry factors. 

proof: 

This property follows directly from the fact that K{zo, z) = K{zo, z). □ 
3.1 Examples. 

Let us present some examples of correlation functions and free energy for low orders. 
3.1.1 3-point function. 




= ResK(p,g) [B{q,pi)B{q,p2) + B{q,pi)B{q,p2)] 
-2ResK{p,q) [B{q,p^)B{q,p2)] 
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Res 



B{q,p)B{q,pi)B{q,p2) 
dx{q) dy{q) 
B{ai,p) B{ai,pi) B{ai,p2) 
2dy{ai) dzi{aiY 



(3 - 13) 

where Zi{z) — ■\/ x{z) — z{ai) is a local coordinate near Oj. 
3.1.2 4-point function. 




+ 5 permutations of (1,2,3) 



(3 - 14) 



+ 5 permutations of (1,2,3) 



Res Res K{p,q)K{q,r) [B{q,pi)B{r,p2)B{r,ps) 

q—^a r—^a 

+B{q,Pi)B{f,p2)B{r,ps) + B{q,p2)B{r,pi)B{r,p^) 
+B{q,p2)B{r,pi)B{r,p3) + B{q,ps)B{r,p2)B{f,pi) 
+B{q,P3)B{r,p2)B{r,pi)] 

+ Res Res K{p,q) K{q,r) [B{q,pi)B{r-,p2)Bir,p^) 

+B{q,pi)B{f,p2)B{r,p3) + B{q,p2)B{r,pi)B{r,p.) 
+B{q,p2)B{f,pi)B{r,p3) + B{q,p^)B{r,p2)B{f,pi) 
+B{q,ps)B{f,p2)B{r,pi)]. 



3.1.3 1-point function to order 1. 



ip) 




(3 - 15) 



Res K{p, q)B{q, q). 

q—*a 



3.1.4 1-point function to order 2. 




) + P 
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= Res Res Res K{p,q)K{q,r)K{r, s) B{q,r)B{s, s) 

5— >a r— >a s— >a 

+ Res Res Res K{p,q)K{q,r)K(f, s) B{q,r)B{s,s) 

g— >a r^a s—*a 

+ Res Res Res K{p,q)K{q,r)K{r, s) B{q,s)B{s,f) 

g— >a r— >a s— >a 

+ Res Res Res K{p,q)K{q,r)K{r, s) B{q, s)B(s,f) 

(J— >a ) — >a s— ►a 

+ Res Res Res K{p,q)K{q,r)K{q, s) B{f,r)B(s,s) 

5— >a i — >a s— >a 




where the last expression is obtained using lemma 13.11 

3.1.5 Free energy F2. 

The second free energy reads 

-2F2 = 2 Res Res Res Res <!>{p) K{p,q)K{q,r)K{r, s) B{q,f)B{s,s) 

p— >a q— >a r— >a s— >a 

+2 Res Res Res Res (!>{p) K{p,q)K{q,r)K{r, s) B{q,s)B{s,f) 

p— >a g— >a r— >a s— >a 

+ Res Res Res Res ^{p) K{p,q)K{q,r)K{q, s) B(f,r)B(s,s). 

p^a q^a 1 — >a s^a 

(3 - 17) 



3.2 Teichmuller pants gluings 

Every Riemann surface of genus g with k boundaries can be decomposed into 2g — 2 + k 
pants whose boundaries are 3g — 3 + k closed geodesies (in the Poincare metric with 
constant negative curvature) [73]. The number of ways (in the combinatorial sense) 
of gluing 2g — 2 + k pants by their boundaries is clearly the same as the number of 
diagrams of Qjf \ and each diagram corresponds to one pants decomposition. 

Indeed, consider the root boundary labeled by p, and attach a pair of pants to this 
boundary. Draw an arrowed propagator from the boundary to the first pants. Then, 
choose one of the other boundaries of the pair of pants (there are thus 2 choices, left 
or right), it must be glued to another pair of pants (possibly not distinct from the first 
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one). If this pair of pants was never visited, draw an arrowed propagator, and if it 
was aheady visited, draw a non-arrowed propagator. In the end, you get a diagram of 



qI^\ This procedure is bijective (up to symmetry factors), and to a diagram of 
one may associate a gluing of pants. 
Example with k = 1 and g = 2: 



(9) 

k ' 



(2) 




4 Main properties 

So, for every regular spectral curve £ = (£, x, y) (and matrix k ii C has genus (7 > 0) 
we have defined some meromorphic n-forms ujn'^ and some complex numbers Fg = ulf'^ . 
They have some remarkable properties (see [60 j ) : 

• ujn^ is symmetric in its n variables (this is proved by recursion). 

• If 2(7 — 2 + n > 0, then [zi, . . . , Zn) is a meromorphic form (in Zi for instance) 
with poles only at the branch-points, of degree at most 6(7 — 6 + 2n + 2, and with 
vanishing residue. 

• If two spectral curves S = (£, x, y) and £ = (£, x, y) are symplectically equiva- 
lent, they have the same FgS 01 g > 1 (although they do not have the same tu^f^'s 
in general) 



dx Ady = ±dx Ady 



Fg{£,K)=Fg{£,K). 



(4-1) 



if C is of genus (7 = 0, then r = exp C^^N"^ "^^^g) is a formal tau function, it 

g=0 

obeys Hirota's equation. This theorem can be extended to (7 > 0, with additional 
^-functions, see section 14. 6[ 



Dilaton equation, for 2(7 — 2 + n > 0: 

V Res ^Zn+l)oo['l,{zu...,Zn,Zn+l) = {2-2g-n)uJi;>\zl,...,Zn). (4-2) 



This equation just reflects the homogeneity property, i.e. under a rescaling y — >• 
Xy, we have 0;^ ^ X^-^^-'^Jj^l 

The derivatives of ^ with respect to many parameters on which the spectral 
curve may depend is computed below in section 14.31 
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• The ujIt^s have many other properties, for instance their modular behaviour 
satisfies the Holomorphic anomaly equations. 

Let us study those properties in deeper details. 

4.1 Homogeneity 

If one changes the function y{z) Xy{z), i.e. just a rescaling of the curve, then it is 
clear from eq. f|2-40l) that the kernel K is changed to K/ A and nothing else is changed. 
Thus, uJn^ changes as: 

and in particular 

FgiC,x,Xy) = X'~'<^FgiC,x,y). (4-4) 
This implies that Fg is a homogeneous function of the spectral curve, of degree 2 — 2g. 
In particular, if one choses A = — 1 one gets (for g > 2): 

Fg{C,x,-y,K) = Fg{C,x,y,K). (4-5) 

4.2 Symplectic invariance 

It is clear from the definitions, that Fg and Un'^ depend on the spectral curve only 
through the kernels B and K, and the number and position of branchpoints. 

The Bergmann kernel B depends only on the underlying complex structure of the 
Riemann surface £, thus it remains unchanged if we change the functions x and y, as 
long as we don't change C 

The kernel K, depends on the functions x and y, only through the combination: 

{y{z)-y{z))dx{z). (4-6) 

Therefore K remains unchanged if we don't change this combination. 

In particular, the kernel K, and therefore Fg and remain unchanged, if we 
change: 

• y ^ y + R{x) where R{x) is some rational function of x. 

• y ^ \y, X ^ x/X where A G C*. 

^ ^ _^ ax+b _^ (cx+d)^ 
cx+d^ ^ ad— be 

Those transformations, form a subgroup of the symplectomorphisms. Indeed, in all 
those cases, the symplectic form dx A dy is unchanged. 

In order to have invariance under the full group of symplectomorphisms, we need to 
prove the invariance under the | rotation in the x, y plane, i.e. x ^ y, y ^ —x, which 
also conserves dx A dy. Using homogeneity eq. (14-51) . we see that this is equivalent to 
consider the invariance under x ^ y, y ^ x. 

This transformation however, does not conserve K, it does not conserve the number 
of branchpoints, and it does not conserve the a;i^^'s with n > 1. However, it was proved 
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in [61] that it does conserve the Fg^s. The proof of [61] is very technical. It is inspired 
from the loop equations for the 2-matrix model. It amounts to defining some mixed 
n + m-forms Unjn, where x and y play similar roles, and for which o^^^q coincides with 

the ujn^ for the spectral curve (£, x, y), and cJq^^ coincides with the uii^ for the spectral 

curve {C,y,x). In particular Fg = Wq^q is both the Fg = for the spectral curve 

{C,x,y), and the Fg = uj^^ for the spectral curve {C,y,x). The proof of [61] relies on 

the fact that cj^+i „i + ^njn+i exact form. 

That leads to: 

Theorem 4.1 Symplectic invariance 

The Fg 's are invariant under the group of symplectomorphisms generated by: 

• y y + R{x) where R{x) is some rational function of x. 

• y — > Xy, X x/X where A G C*. 

cx+d' " ad— be 

• X ^ y, y ^ —X. 

In addition, the Fg 's are also invariant under: 

• X ^ X, y ^ —y. 

This theorem is a powerful tool which allows to compare the F^'s of models which 
look a priory very different. We will see examples of applications in section [lO.li 

The cui^'^'s with n > 1 are not conserved under symplectic transformation, instead 
they get shifted by exact forms. 



4.3 Derivatives 

In this section, we study how the F^'s and uon^^s change under a change of spectral 
curve, and in particular under infinitesimal holomorphic changes. 

Consider an infinitesimal change + t5y at fixed x, or in fact it is more 

appropriate to consider the variation of the differential form ydx: 

ydx ydx + eS{ydx) + O(e^) = ydx + tdVl + O(e^) (4-7) 

where dVl is an analytical differential form on an open subset of /fl. If instead of 
working at fixed x, we prefer to work with some local parameter z, we write: 

5y{z)dx{z) - 6x{z)dy{z) = dQ{z). (4-8) 

This shows that the set of holomorphic deformations of the spectral curve is equipped 
with a Poisson structure, but we shall not study it in details in this article, see [131 HI] 
for the Frobenius manifold structure. 

Classification of possible 1-forms dVl: 

The deformation dQ is a 1-form. Here we shall consider only meromorphic deforma- 
tions, and meromorphic 1-forms are classified as 1st kind (no pole), 2nd kind (multiple 
poles, without residues), and 3rd kind (only simple poles). 

^Note that ydx does not need to be a meromorphic form itself to be able to consider such defor- 
mations. 
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First kind deformations are holomorphic forms on £, i.e. they are linear combi- 
nations of the dui^s (see section [2.2.ip : 



dui{z) = ^(f B{z,z') (4-9) 



where Bi = Bi - '^ij^j 



• 2nd kind deformations have double or multiple poles. They can be taken as linear 
combinations of Bergmann kernels or of their derivatives. Choose a point p E C 
If X is regular at p, choose the local parameter ^{z) = x{z) — x{p), and if x has 
a pole of degree d at p, choose ^{z) = x{z)~^/^, and define: 

Bkiz; p) = Res B{z, z') ^{z'yK (4-10) 

z'— >p 

All 2nd kind differentials are linear combinations of such Bk{z\p). 

• 3rd kind differentials have only simple poles, and since the sum of residues must 
vanish, they must have at least 2 simple poles. Choose two points pi and p2 in 
the fundamental domain, and define: 

dSp,,pM= / B{z,z'). (4-11) 

All 3rd kind differentials are linear combinations of such dS. 

Theorem 4.2 A general meromorphic differential form dQ with poles pk's, can be 
written: 

9 

dQ{z) = 2z7r J]5eidui(2;) + ^5tfcci5p„o(z) + J] J]54,j-5j(z;pfc). (4-12) 

i=l k k j 

It can be noticed that the coefficients 6ei, 6tk, St^j are the flat coordinates in the metrics 
of kernel B, of the corresponding Frobenius manifold structure. 

proof: 

Indeed, let pk be the poles of dfl, and write the negative part of the Laurent series 
of dQ near its poles as: 

dniz) ~ 5tk^-yj6t,,^^. (4-13) 

We see that 

dn{z) -J2^ikdSp,,oiz) - ^^54,,5,(2;;pfc) (4-14) 

k k j 

is a 1-form which has no poles, thus it is a holomorphic form, and it is a linear combi- 
nation of the dui^s. □ 



25 



If K = 0, i.e. if B is the Bergmann kernel, it is normalized on the .4-cycles, and 
we have: §^dS = 0, ^^-B^ = 0. Thus, the 6ei are easily computed as Sei = ^ £^ dfl. 
However, if k 7^ 0, this is no longer true, we have 6e = £4 ~ 2^ §b 
variations 5tk,j or 5tk get mixed with the ^e^ through the variations of r. The good 
way to undo this mixing, is by defining a covariant variation: 

Definition 4.1 Covariant variation: 

Ddn = 5dn + tr 5dnT k -^-^ (4-15) 

where Sdn^ is the variation of the Riemann matrix of periods t under ydx ydx + edQ. 
Derivatives with respect to k are studied in details in section 14.4.11 below. 
The important point, is that dQ can always be written as: 

dn{z) = [ B{z, z') ^{z') (4-16) 

where dVL is some continuous path (a chain or a cycle, which is related to the Poincare 
dual of dVl) on £, and h.{z') is an analytical function defined locally in a vicinity of dVL. 

The theorem is then: 

Theorem 4.3 Variation of the spectral curve: 

Under an infinitesimal deformation 6ydx — 6xdy = dfl{z) = fg^ B{z,z') A{z'), 

the ujn^ 's change by: 

Ddnuji'\zu...,zn)= [ ^i^i(^i,...,^„,/) A(z')- (4-17) 

Jan 

For example, in particular with n = we have: 

DdnFg= [ u['\z') A{z'). (4-18) 
Jan 

4.3.1 The loop operator 

This theorem can also be restated in terms of the "loop operator", which corresponds 
to dn{z) = B{z,z'). We define: 

Definition 4.2 The loop operator is: 

D.' = DB(.,^r (4-19) 

It satisfies: 

D,,uji^\z^, ...,Zn)= coi'liz,, . . . , z'). (4-20) 

The loop operator is a derivation, i.e. Dz'{uv) = uD^'V + vD^'U, and it is such that: 
D,,D,, = D,,D,„ and D,,^ = ^D^,. 

In random matrix theory, the loop operator [6] , is most often written as a functional 
derivative with respect to the potential V{x): 

d dcf 



dV{x{z')) 



D,,. (4-21) 
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4.3.2 Inverse of the loop operator 

The loop operator allows to find u^l^li in terms of uJn \ i.e. increases n by 1. The inverse 
operator, which decreases n by 1 can also be written explicitly: 

Theorem 4.4 Let $ be a primitive ofydx, i.e. a function defined on the fundamental 
domain such that d^ = ydx, then we have, if 2 — 2g — n < 0: 

{2-2g- n) ^^^(^i, . . . , = V Res ...,z^,z) $(^). (4-22) 

i 

This theorem is easily proved by recursion on 2g + n — 2. 
This theorem is at the origin of definition 12.61 for n = 0. 



4.4 Modular properties 

In section 12.2.31 we have introduced a deformation of the Bergmann kernel with a 
symmetric matrix n of size g x g. The reason to introduce this deformation, was that 
it encodes the modular dependence of the Bergmann kernel, i.e. how the Bergmann 
kernel changes under a change of choice of cycles A, B. Thus, studying the modular 
dependence of the F^'s and oon^'s amounts to studying their dependence on k. 

Also, in section 1131 we have seen that the covariant derivative involves the compu- 
tation of derivatives with respect to k. 



4.4.1 Dependence on k 



Since the kernels B and K depend linearly on k, all the stable Un^'s and F^'s are 
polynomials in k, of degree 3g — 3 + n. 

Notice that dB{zi, Z2)/dKij = 2in dui{zi) duj{z2) is factorized, i.e. a function of zi 
times a function of Z2. This simple observation, together with 



duiiz) 



1 

2z7r /k 



Biz,z') 



13 i 13 i ^ ^ 



(4-23) 



leads, by an easy recursion, to the following theorem: 
Theorem 4.5 For 2 — 2g — n < : 



2i7idu';f\zi, . . .,Zn)/dKi^. 



Bi 

9 I 



dz (h dz 



{<?-!)/ / 



h=0 ICJ 

(4 - 24) 

where J = {zi, . . . , Zn}, and J2h Y^i "^eans as usual that we exclude (h, I) = (0, 0), (g, J). 
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This theorem can be apphed recursively, to compute higher derivatives, and even- 
tually recover a polynomial of n by its Taylor expansion at k = 0, of the form: 

W='j2 l^i^)' 9'F,\^=o. (4-25) 

fc=0 

According to theorem 14.31 of section fO| the i3-cycle integrals, computed at k = 0, 
are the derivatives with respect to coordinates e, of eg. (14-121) : 

S9)l 



^ ^ Q ^ v-^, ^ . ^(._)^(^^^ . . . , z) (4-26) 



and therefore we have: 



2 dei dej 

9 



h=0 ICJ ^ 



(4 - 27) 



We can thus trade the k dependance into derivatives with respect to the coordinates 
ej. For instance we have at g = 2, and with g = 1: 



-2 , ^ ,^3 



F2{k) = + I {d^F, + dF, + ^ (d'Fo + 4 d'Fo + ^ (lO d^Fo d^Fo) 

_ (4-28) 
where Fg = Fg{K = 0), and (9 = d/de. 

This result is best interpreted graphically. For example with g = 2 we have: 



^ ' 2 2 8 2 



+ (4-29) 

where each line with endpoints is a factor Hij, and each connected piece of Rie- 
mann surface of genus h, with k punctures zi, . . . , Zfc is a d^Fh/dei^ . . . dei^. Each graph 
is a possible " stable" degeneracy of a genus g Riemann surface (imagine each link con- 
tracted to a point), stability means that each connected component of genus h with k 
marked points must have 2 — 2h — k < 0. The prefactor is 1/^Aut, i.e. the inverse 

of the number of automorphisms, for instance in the last graph we have a 

Z2 symmetry by exchanging the 2 spheres, and a cr^ symmetry from permuting the 3 
endpoints of the edges, i.e. 12 = #(^2 x 0-3) automorphisms. 

More generally, by a careful analysis of the combinatorics of the d^-^s, one can see 
(this was done in [53], and coincides with the diagrammatics of [2]) that the Taylor 
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expansion eg. (14-251) . reconstructs the expansion of a formal Gaussian integral (i.e. order 
by order in powers of A^): 



e ^ = / drji . . . drjg e ^ e « 

e e (4-30) 

and the graphical representation above is just the Wick's expansion of the gaussian 
integral. 

This diagrammatic expansion of modular transformations was first derived in ^ in 
the context of topological strings. 



4.4.2 Holomorphic anomaly 

In particular, if we choose n to be the Zamolodchikov Kahler metric k = (r — r)^^, we 
have seen in section 12.2.51 that the Bergmann kernel becomes the Schiffer kernel and 
is modular invariant, which means that it is independent of the choice of cycles A, B. 
Since the only modular dependence of the F^'s and Un^^s is in the Bergmann kernel, 
we have: 

Theorem 4.6 // k is the Zamolodchikov Kdhler metric k = {j — t)~^ , then Fg and 
Un^ 's are modular invariant. 

The price to pay to have modular invariant -F^'s, is that they are no longer analyti- 
cal functions of r, i.e. analytical functions of the spectral curve, and in particular they 
are no longer analytical functions of the e^'s. However, since the only non-analytical 
dependence is polynomial in k, and k"^ is linear in r which is the only non-analytical 
term, and since we have relationships between derivatives with respect to k and deriva- 
tives with respect to e, by a simple chain rule, we obtain the following theorem [63] : 



Theorem 4.7 The Un^ 's satisfy the Holomorphic anomaly equations 



dli 



-1 d^Fp 1 



2, ,(9-1) 



h=0 ICJ 



dt de 



de 



de 



(4-31) 



In particular for n = 0.- 



dF^ 



-1 



(2i7r) 



a^Fo 1 



9e2 



de 



9-1 

E 

h=l 



OF, OF 



g-h 



de de 



(4-32) 
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This equation was first found by Bershadsky, Cecotti, Ooguri and Vafa (which we 
refer to as BCOV |T7]) in the context of topological string theory. Here we see that 
the symplectic invariants Fg always satisfy this equation, and it is tempting to believe 
that the symplectic invariants Fg, should coincide with the string theory amplitudes, 
i.e. the Gromov-Witten invariants. This question is debated below in section [TTl 
Unfortunately, the holomorphic anomaly equations do not have a unique solution, and 
although this conjecture is almost surely correct, no proof exists at the present time, 
apart from a very limited number of case^. 

Let us briefiy sketch the idea of BCOV. String theory partition functions represent 
"path integrals" over the set of all Riemann surfaces with some conformal invariant 
weight. In other words, they are integrals over moduli spaces of Riemann surfaces of 
given topology, and topological strings are integrals with a topological weight, they 
compute intersection numbers of bundles over moduli spaces (see [9T| 1115] for intro- 
duction to topological strings). 

Moduli spaces can be compactified by adding their "boundaries", which corre- 
spond to degenerate Riemann surfaces (for instance when a non contractible cycle gets 
pinched, or when marked points come together). The integrals have thus boundary 
terms, which can be represented by (5-functions, and (5-functions are not holomorphic. 
In other words, string theory partition functions contain non-holomorphic terms which 
count degenerate Riemann surfaces. 

On the other hand, if one decides to integrate only on non-degenerate surfaces, one 
gets holomorphic patition functions, but not modular invariant, because the boundaries 
of the moduli spaces are associated to a choice of pinched cycles. Modular invariant 
means independent of a choice of cycles. 

To summarize, the holomorphic partition function is obtained after a choice of 
boundaries, i.e. a choice of a symplectic basis of non contractible cycles Ai fl Bj = 6ij, 
and cannot be modular invariant. The modular invariance is restored by adding the 
boundaries, but this breaks holomorphicity. 

There is thus a relationship between holomorphicity and modular invariance. 

4.5 Background indepedence and non-perturbative modular 
invariance 

We have seen in the previous section, that the F^'s are not modular invariant, un- 
less we choose k = {t — r)~^, i.e. modular invariance can be restored by breaking 
holomorphicity. 

In fact, there is another way of restoring modular invariance, without breaking holo- 
morphicity. It exploits the fact that the modular transformations of Fg^s, i.e. eq. (14-241) . 
is very similar to the modular transformation of theta-functions. It was shown in [M] , 
that certain combinations of ^-functions and Fg^s, are modular, and reconstruct a non- 

®In [22], this conjecture is proposed as a new definition of the type IIB topological string theory. 
The interested reader may find all the details of this conjecture as well as numerous checks in this 
paper. 
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perturbative, modular partition function, which is also a Tau- function (see section 
and which has a background independence property. 

Consider a characteristics (/i, i^), and a spectral curve £ = {C,x,y), choose k = 0. 
Following [55l [Mj . we introduce a nonperturbative partition function by summing over 
filling fractions, defined by 

Zsin, ly; e) 

o > (2-2hi-li) 



e 



(4 - 33) 



where is the sum of all terms contributing to order A^"-' . In this partition function, 
the Fg^s are the symplectic invariants of the spectral curve S, their derivatives are with 
respect to the background filling fraction e and computed through theorem 14.31 at: 

e = — (h ydx. (4-34) 

Notice that the Fg{Sys and their derivatives depend on the choice of a symplectic basis 
of 2g one-cycles Ai, Bj on C Finally, the theta function 6^,,^ of characteristics (/x, z/) 
is defined by 

and is evaluated at 

NF^, F^ = j^y{x)dx, ^ = ^^0- (4-36) 



u 



In (14-331) . the derivatives of the theta function (I4-35P are w.r.t. u. The derivatives 
3f B and tl 

\®t}vFQ F[ actually means: 



of B and the derivatives of Fg, are written with tensorial notations. For instance 



le(^)F'"F' = ^ V ^^©.. d^F, OF, 

g M,- 1 2!3!1! . ^. du,,du,,du,,du,, de,,de,,de,, de,, ^ ' 

^1 )'^3 5^4 



and the symmetry factor (here | = grin) is the number of relabellings of the indices, 
giving the same pairings, and divided by the order of the group of relabellings, i.e. 
/c! /i! . . . /fc!, as usual in Feynmann graphs. 
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The 6 function above is closely related to the standard theta function, which is 
defined by 

= exp[z7r(n + /x)r(n + ft) + 2m{n + + i^)]- (4-38) 
nezs 

It it easy to see that these two functions are related as follows 

e,,.(«,r) = exp[-iV2(eF^ + ^e^F^')] ^^Ir) (4-39) 



where 



4.5.1 Modularity 



N f / F' \ 



Theorem 4.8 All the terms Zj in eg. lji4-33{ ) are modular, i.e. they transform as the 



characteristics (/i, z/). More precisely, if we make a modular change of cycles A, B 
A, B, we have: 

e^^^«+^^e,,,(ti, f) = Ca e^^^o+^^e,,.(«, r) (4-41) 

and for all j > I: 

Z,(/i,z>) = Z,(/i,i.). (4-42) 



This theorem was proved in [64], mostly using the diagrammatic representation of 
section [3l and the diagrammatic representation of [2]. 

For example the following quantities are modular: 

6' 1 6"' 

Zi = -^F{ + i-^F-, (4-43) 

Z2 = F, + l^F[' + l^F- + ^^^r + l^KF[ + ^^(O^- (4-44) 
2 6^,1, 2(3^^^ 24 6^^1, 6 B^,i, 72 B^_^ 

4.5.2 Background independence 



Theorem 4.9 The partition function eg. li4-33\ ) is independent of the background filling 



fraction e, i.e., for any two filling fractions ei and 62: 

Zeifi, I/, ei) = Z£-(/x, z/, 62). (4-45) 

This theorem follows directly from the definition eg. (14-331) . It has important conse- 
quences which we shall not study here [6l] . 
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4.6 Integrability 



Out of the -Fg's, one can define a "formal tau-function" . In this section, let us assume 
that L = P^, i.e. it has genus g = 0. The higher genus case is discussed in section 14.6.61 
below. 



Definition 4.3 The formal r-function is defined as a formal series in a variable N: 

oo 

luTN = J2^^~^' ^9- (4-46) 
9=0 

Now, we shall explain why it makes sense to call it a r-function. r-functions are 
usually defined in the context of integrable systems, and they have several more or less 
equivalent definitions, see [H]. 

One possible definition of r-functions relies on Hirota equations [72|[Ti]. and another 
one relies on a free fermion representations, i.e. determinantal formulae [75| [86 | [Ti]. 



4.6.1 Determinantal formulae 

In the following of this section, most of the functions have an obvious formal de- 
pendence. For the sake of brevity, we do not write it explicitly as long as it is not 
needed. 

Out of the ci;i^'''s, it is convenient to define the formal series: 

u^nizu ...,.„) = - ^-ff'^f^'^l + Y: iV-^-"4^)(z„ . . . , z„) (4-47) 

and also the "non-connected" correlators: 

n k 
k=l /iU...U4=J i=l 

For example: 

^^2(^1, Z2) = UJ2{zi, Z2) + Ui{zi)uJi{z2), (4-49) 



(^3(^1, 2:2, 23) = i^?.{.Zl, Z2, Z:i) + UJi{zi)uJ2{z2, Z'i) + UJi{z2)uJ2{zi, Z^i) 

+UJi{z-i)uJ2{Zl, Z2) + UJi{Zi)Ui{z2)uJl{z-i). (4-50) 

In other words, the uJn are the cumulants of the a7„'s. 

The following proposition is proved in some cases (hyperelliptical spectral curves 
|15]). and in all matrix models, however, it is expected to hold for any spectral curve: 

Proposition 4.1 There exists a (formal) kernel H{zi, Z2) , such that: 

u^iz) = hm (Hiz,z') - ^^I^] , (4-51) 

z'^z \ Z — Z' I 
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I \ u( \ul \ dx{zi) dx{z2) , . 

0J2[zi, Z2) = -H{zi, Z2)H{z2, zi) - -— (4-52) 

[x[Zi) - x[z2)y 

and if n > 3; 

uJn{zi, ...,Zn)=" det "i,j=i,...,n{H{zi, Zj)) (4-53) 

where the quotation mark " det " means the following: write the determinant as a sum 
over permutations of products of H's: (iei{H{zi, Zj)) = ''^^{—lYY\H{zi, Za-(i)) ■ Then, 

cr i 

every time a permutation has a fixed point a{i) = i we must replace H{zi, Zi) by uJi{zi), 
and every time a permutation has a length 2 cycle a{i) = j, a{j) = i we must replace 
the factor H{zi, Zj)H{zj, Zi) by —uj2{zi, zj). 

This is equivalent to saying that for n > 3, the cumulants are given by: 

n 

Un{Zi,...,Zn)= (-'^yIlH{Zi,Z„^i)). (4-54) 

cyclic cr i=l 

Example: 

^3(2:1, Z2, Zs) = H{zi, Z2)H{Z2, Z3)H{Zs, Zi) + H{zi, Z3)H{Z3, Z2)H{Z2, Zi). (4-55) 

The determinantal formulae for correlation functions were first found by Dyson and 
Mehta [IHl ES] in the context of random matrix theory, and have led to a huge number 
of applications. 

Moreover the kernel H can be written rather explicitly. In all matrix cases, the 
kernel H for the determinantal formulae above, coincides with the kernel H which we 
define below: 

Definition 4.4 We define the formal kernel H as a formal spinor in zi and Z2, given 
by an exponential formula 



n] JZ2 ■Iz2 



This exponential formula for the kernel is to be understood order by order in powers 
of A^, namely: 



-N ydx 

H(zi,Z2) = 



E{Zi, Z2 
jy-l /"^i /"^i /"^i 



(1) 



6 

where E{zi,Z2) is the prime form: 



1 + 

cuf +0{N~'') , (4-57) 



22 



Z2 J Z2. J Z2 



Eiz,,Z2) = ^^. (4-58) 
y/azi az2 
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For example, one of the terms contributing to H to order ^ is: 

/o) - V ~ (4-59) 



■ _ — (^1 - Z2f 

^ ~ ^ - y'{ai) x"{a.i) (a^ - z^f (a^ - 2:2)3 ' 



Z2 J Z2 J Z2 



In all matrix model cases, the kernel H can be written as a Sato-formula, and 
coincides with iY, but this is not proved in general. 

4.6.2 Examples 

For example, if we consider the Airy curve y = ^Jx^ i.e. E = (P^, z'^,z), we find that H 
is the Airy kernel: 

Ai{zl)Bi\ zl) - Ai '{zl)Bi{zl' 

Z-i Zo 



H{zi, Z2) = "^^"^ — ^7'^^' "'^'^ Vzidzi Z2dz2. (4-60) 



The corresponding Baker- Akhiezer function is the Airy function, and the correlators 
uJn, are the correlators given by the determinantal Airy process. 

4.6.3 Sato formula 

The theorem 14. 3[ implies that under an infinitesimal change of spectral curve of the 
3rd kind eg. (14- 11 1) : 6ydx = t rfS^^^^j, we have: 

^i^^^^f^iiiiJil- (4-61) 

and thus: 



22 



■' Z2 J Z2 



•••,<)• (4-62) 



The exponential formula of proposition 14.41 is nothing but the Taylor expansion of 
Fg{t) computed at t = in terms of derivatives taken at t = 0, i.e. Fg{N~^) = 
^^^j^ d^Fg{0). In other words: 

k 

Theorem 4.10 

H(z,, Z2) = r^(/:,a:,y + ^%^) ^ 

^' ^ x{zi)-x{z2) TN{C,X,y) 

This theorem can be interpreted as Sato's formula |106j for integrable systems. In the 
context of random matrix theory, it can be interpreted as Heine's formula |109j . 
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4.6.4 Baker-Akhiezer functions 

Let ai, . . . , am be the poles of the function x{z), of respective degrees di, . . . , dm- Since 
X is a meromorphic form of degree d = "^^di, there are d sheets, i.e. d points on C, 

i 

z^{x), . . . , z'^{x), such that x{z'^) = x. The following matrix: 

n{xuX2) = (h{z^{xi),z'{x2))) (4-64) 

is a square matrix of size d x d. 

The ^f-function of the Lax system [13], is obtained by choosing X2 = oo, i.e. z^{x2) = 
ai. Since some poles a^ are multiple poles, in order to get an invertible matrix, we take 
linear combinations of rows, and we define: 

dU^)) TIE! 

(4-65) 

where o is an arbitrary basepoint, and C,i{z2) = x{z2)^^^^^ is the local parameter in the 
vicinity of a^. 

Those functions are the Baker-Akhiezer functions. 

We also have d couples I = with i = 1, . . . , m, j = 1, . . . ,di, and thus the 
following matrix is a square matrix: 

nx) = {Mz\x)))^^^_^^^^^_^. (4-66) 

It is the \E'-function of the corresponding Lax system |14j . 



i = 1, . . . , m, j = 1,. . . ,di = lim 



4.6.5 Hirota formula 

Notice that H{zi, Z2) has a simple pole at zi = Z2 and behaves like: 



H{z^,Z2)r^^^^^ (4-67) 

Zi - Z2 

near zi = Z2, and this holds for any {g = 0) spectral curve S = {C, x, y). In particular, 
we have, for any two such spectral curves £ and £: 

Theorem 4.11 

Res H{zi, z] S)H{z, Z2; S) = H{zi, Z2; S). (4-68) 

If we consider that H is given by the Sato formula of theorem I4.1UI this theorem is 
precisely the Hirota equation for the r-function tn [H [1]. 

This theorem justifies that we can call tn a Tau-function. By expanding locally 
the dSzi,z2 in the vicinity of poles of x, we can see that it is the multicomponent 
Kadamtsev-Petviashvili (KP) tau-function. There is one set of component for each 
pole ttj of X. In the case where £ is an hyperelliptical curve, of type y"^ = Pol(x), the 
function x has two poles, which are symmetric of oneanother, and everything can be 
written in terms of the expansion near only one pole. In that case tn reduces to the 
Kortweg-de-Vries (KdV) tau-function [14j . 
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4.6.6 Higher genus 

So far, in tliis section, we were considering genus zero spectral curves, i.e. z G C, and 
x{z) and y{z) analytical functions of a complex variable. 

The integrability relied on the Sato formula, which gives the kernel H as the r- 
function of a shifted spectral curve, i.e. the exponential formula. 

For higher genus ^ > 1, the problem is that the exponential formula does not define a 
well-defined spinor on £. Indeed, £ is not simply connected, and the abelian integrals 
. . f^^ ojn are multivalued functions of zi and Z2 because there is not a unique 
integration path between zi and Z2. The exponential formula has to be modified. It 
was proposed to modify it with some theta functions (see section . 

Definition 4.5 Given a characteristics (/U, z/), the "tau- function" is defined by the 
non-perturbative partition function of section \4^ 

Tat (/i, I/, £:) = (/i, z/) . (4-69) 
Then define the spinor kernel through the Sato formula: 

Definition 4.6 

^ /^^^ _ \/dx{z^)dx{z2) rjv(/i, z/, £, X, y + ^ ^ 
''^'"^ ^' ^ x{zi)-x{z2) rAr(/i, z/,£,x,?/) 

With this definition, Hf^^^y) is closely related to the Szego kernel |109] . 

Theorem 4.12 H(^^^y){zi, Z2) is well defined for Zi, Z2 G C. 
proof: 

Integrals of a;i^'''s are in principle defined only on the universal covering of £, and 
one needs to check that after going around an ^-cycle or i3-cycle, H(^^^^){zi, Z2) takes 
the same value. 

Notice that, if zi goes around an ^-cycle, then dSz^^z2 is unchanged, and if zi goes 
around the cycle then dSz^^z2 is shifted by a holomorphic differential: 

dSzi+B,,z2 = dSz^^zi + 2i7idui. (4-71) 

However, it was proved in [55j that the r function above is background independent, 
which exactely means that, for any A: 

T]\f{fi, u, C,x,y + Xdui/dx) = Tj^{fi, £, x, y) (4-72) 

and therefore, we see that H(^^^^){zi, Z2) is unchanged if zi goes around a i3-cycle. □ 
Then, we see that 



Theorem 4.13 H(^^^y^[zi, Z2) obeys the Hirota equation: 

Res H(^^^y){zi, z\ £)H(^^^^){z, Z2\ £) = H(^^u){zi, Z2; £). (4-73) 



*Z2 



If we consider that H(^^^y^ is given by the Sato formula of theorem l4-70[ this theorem 
is precisely the Hirota equation for the Tau-function tn- 
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4.7 Virasoro constraints 



Its has been understood for a long time that the random matrix integrals are funda- 
mentally linked to Virasoro and W-algebras through differentials equations on their 
moduli called Virasoro or W-constraints. The definition of the symplectic invariants 
and of the correlation functions themselves were inspired by these constraints since 
they mimic the solution of the loop equations of random matrix models, the latter 
being considered as equivalent to the Virasoro constraints. 

In a series of papers [3 [HI [9], Alexandrov, Mironov and Morozov go even further 
and propose to generalize the notion of random matrix integrals by defining a general 
string partition function interpolating between different matrix models. This partition 
function is characterized as a "D- module" solution of some Virasoro constraints. 

It is natural to see the symplectic invariants and the r-function built from them as 
a good candidate for this string partition function. It is thus interesting to clarify the 
arising of Virasoro constraints in the theory of symplectic invariants by looking at the 
variations of the latter wrt the moduli of the spectral curve. 

4.7.1 Virasoro at the branch points 

One can slightly rewrite the recursive relations defining the correlation functions eq. fl2-44p 
by moving all the terms to the same side of the equation. One gets: 

9 

,(9-1)/^;. n_LY^Y^,,W n,S9-h) 



= J]Res ir(zo,^) 0.1^2 (^,^,J) + J]J]<V|(^,^)<^^^^^^^^ 

z L 
i h=0 IdJ 

+ydx{z)uj^^l-^{z, J) + ydx{z)uj^Jll^{z, J) 

(4 - 74) 



By summing over the genus g and interpreting the correlation function as the re- 
sult of the loop insertion operator on the symplectic invariants, this equation can be 
rephrased Virasoro constraint: 

Theorem 4.14 For any point z on the spectral curve, the partition function is a zero 
mode of the global Virasoro operator V{z) 

V{z)tn = (4-75) 

with 

V{z) =Y,i K{z, z') : J{z')J (7) : (4-76) 

i ^"-^ 

and the global current is defined by: 

J{z) = Nydx{z) + j^D, (4-77) 
for any point z of the spectral curve. 
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This means that the recursive definition of the correlation functions is nothing but a 
Virasoro constraint on the r-function defined globally on the spectral curve. 

Let us now approach a particular branch point and blow up the spectral curve 
around this point (see section . A rational parametrization of the blown up curve 
can read 



x(z) = 



{i)^k (4-78) 



k ^' 



A;=0 



where the T^^*^'s are the coefficients of the Taylor expansion of ydx around the branch 
points 

oo 

ydx{z) = yTf^^'^\z)dUz) (4-79) 

with the local coordinates 



fc=0 



Uz) = ^/xiz)-x{a^). (4-80) 



From section ITSl one knows that the projection of the correlation functions in the local 
patch around aj built from the local parameter C,i is given by the correlation functions 
of the blown up curve, i.e. the correlation functions of the Kontsevich integral with 
times tI:^\ k = 0, . . . , oo. Since the recursive definition of these correlation functions 
is equivalent to Kontsevich Virasoro constraints this means that the global Virasoro 
operator projects to the continuous Virasoro operator in this local patch of coordinates 
around Oj: 

Theorem 4.15 For any branch point ai, the partition function is a zero mode of a set 
local Kontsevich Virasoro operator Vi{^i{z)) for any point z in a neighborhood of a^: 

V^, V,{i,{z))T^ = Q (4-81) 



where ii{z) = \/ x{z) — x{ai) and the Virasoro operator annihilates the Kontsevich 
T-function: 

m)ZKiTt'^) = (4-82) 

where 

with Kontsevich times 



ZKin") ■= / e " V— -'^ ^"J (4-83) 

J formal 

r«:=iTrA-^ (4-84) 



Indeed, as it is exhibited in section 110.11 the corresponding spectral curve has only one 
branch point and all the moduli of the integral are summed up in the Taylor expansion 
of the differential form xdy at this branch point. 



^See [SI [30] for detail on these continuous Virasoro contraints. 
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4.7.2 Loop equations and Virasoro at the poles 

In the preceding section, we liave translated the recursive definition of the correlation 
functions into a set of Virasoro operators related to the moduli of the spectral curve at 
the branch points. One can proceed in a similar way for the moduli at the poles of the 
one form ydx by building a set of equations solved by the correlation functions called 
"loop equations" since they mimic the loop equations of random matrix theory. 

Theorem 4.16 The correlation functions u;„ are solutions of the loop equations: 

k ^ 

^i+i{z, 7.i,)ujk-i+i{z, zk\l) + ^^^+2(2;, zk) = Pi,k{z, ZK)dx{zY (4-85) 

1=0 

where the function 

k 

y^Ul+l{z', Zi,)Uk-l+l{z', zk\l) + ■W2^k+2{z', z', Zk) 

^-^^ ? i - .iM') 

is a function of z with poles only at the poles of ydx. 



In the matrix model case, these loop equations are often referred to as Virasoro con- 
straints. They indeed encode a set of Virasoro constraints build from the poles of ydx. 
Let us make this assertion clear in the general framework of the symplectic invariants: 

Theorem 4.17 For any point z E C, the r-function satisfies 

V{z)tn = (4-87) 

where one defines the global Virasoro operator 

and ^i{z) = 1^ is a local parameter in the neighborhood of (see eg. li4-iC\ ) for the 
definition of^{z)). 

The r-function can thus be seen as the zero mode of another Virasoro operator 
globally defined on the spectral curve. This new operator, equivalent to the loop 
equations, can be easily projected to a set of local Virasoro operators in the vicinity of 
the poles of ydx instead of the branch points for the first one. In order to follow this 
procedure one has to restrict to ydx which are holomorphic forms with poles ctj such 
that ydx{z) r^^^a, '^kti^k^^{z)d^i{z). 

k 
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Theorem 4.18 For any point z in the, neighborhood of a pole ai of ydx: 
where the local Virasoro operator is defined as the loop operator 



V^\z) :-- 



dUz') 



with the local current 

J">(z) 



E 

k>0 



kti^kUzf'^dii{z) + 



■.J'^\z'f: 



dUz) d 



(4-89) 



(4-90) 



(4-91) 



Remark that these local Virasoro operators are indeed Laurent series in ^i{z) with 
only negative powers whose coefficients are differential operators satisfying Virasoro 
commutation relations: 



v«(^) = 5^v«6(^)-'^(rf6(^))^ 



fc>0 



and 



(4-92) 



(4-93) 



These local operators around the poles have also a natural solution: the one hermitian 
matrix integral 



e * 



Tr 



with a polynomial potential 



V{x) 



formal 



X 



Y^tkx^ 



(4-94) 



(4-95) 



A;=3 



whose coefficients are identified with the moduli at the poles t « jt (see section 15.11 for 
more details). 



4.7.3 Givental decomposition formulae 

Let us suppose in this short section that the spectral curve has genus 0, i.e. £=Riemann 
sphere. In this case, the only moduli of the curve are: 

• either the position of the poles and moduli ^ at these poles; 

• either the position of the branch point and the moduli T^^^ . 

Let us ffist focus on the branch points of the spectral curve. The dependence of 
tn on the moduli at the branch points is constrained by the local Virasoro equations 
eg. (14-8 ip . Thus, this r-function can be decomposed as a product of the zero modes of 
the different local operators at the branch points, i.e. a product of Kontsevich integrals, 
up to a conjugation operator mixing the moduli at the different branch points. 
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Theorem 4.19 ttv can be decomposed into a product of Kontsevich integrals associated 
to the branch points ai : 



(TW,T(2),...) = e" J]Z^(T(')). 



(4-96) 



where the symbol T*-') stands for the infinite family \ T^f^ \ , with the intertwining 



k=0 



operator U defined by 



■ ■ J a-i J a. 



where 



and 



A^'^^\z,z) ■.= B{z,z') - 



dUz)dij{z') 
1 dUz) d 



Nki^{z)dTl^'^' 



(4-97) 



(4-98) 



(4-99) 



One can proceed exactly in the same way by looking at the moduli at the poles: 
this time the decomposition is expressed as a product of 1-hermitian matrix integrals. 

Theorem 4.20 r^v can be decomposed into a product of one hermitian matrix integrals 
associated to the poles of the meromorphic form ydx 



(4-100) 



where ti stands for the infinite set {ti^k}^=Q, with the intertwining operator U defined 
by 

U:=^l (f A^''^\z,z')nj{z')n^{z) 



where 



and 



A^''^\z,z')=B{z, z') 



dUz)di,{z') 



n,{z) ■.= NY,t,,i^{z)dUz) - ^Mrs^ 

^ Nk^f{z)dti^k 



(4-101) 



(4-102) 



(4-103) 



Remark that, in both cases, these decomposition formulae consist in writing a KP 
tau function as a product of KdV tau-functions. Indeed, these formula were already 
derived by Givental in the study of KP tau functions [671 EH] • 
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4.7.4 Vertex operator and integrability 

In this paragraph, we do not consider the Baker- Akhiezer functions as defined in section 
14.61 On the contrary, we define them as the images of the partition function under the 
action of some global operator on the spectral curve. 

Let us define the equivalent of the Baker- Akhiezer (BA) functions: 

Definition 4.7 One defines the x-type global BA and dual BA functions as 

:= exp(fi^(z))Z and := exp (-fi^(z)) Z (4-104) 

where q q q 

fix(g) ■■= j J. = N j ydx{z) j D,. (4-105) 

One also defines the x-type local BA functions as: 

^ ifl{z) := exp {VLi{z)) Z and **o(z) := exp (-(^^(z)) Z (4-106) 

where the operator Vli was defined in eq.^^TOS). 

We finally define the corresponding y-type BA functions: 

^{z) ■.= exp{ny{z))Z and ^*{z) := exp {-ny{z)) Z (4-107) 

where 

^yiq) ■= f Jy = N J'xdyiz) - ^ J' D,. (4-108) 

These functions correspond to deformations of the spectral curve and thus coincide 
with the Baker- Akhiezer functions of eg. (14-651) . 

Lemma 4.1 The BA functions can be written in terms of the partition function as 

*(z) := exp {n,{z)) Z{U) = (4-109) 

*,o(;^) := exp {n,{z)) Z{t,) = ^^^'^^j'^' (4-110) 

^{z) := exp (Qyiz)) Z{U) = ^^^'^-f^^ (4-111) 

and 

*,o(^) := exp {ny,{z)) Zii,) = ^^^'z(tP^ ^^"^^^^ 

where the U are the coefficients of the Taylor expansion of xdy{z) as z ^ ai and [z^^] 
is the usual Hirota symbol 

On the other hand, thanks to the pole structure of the BA functions, one gets 
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Theorem 4.21 The Baker-Akhizer functions satisfy the bilinear Hirota equation 

(4-113) 

where x' and y' are functions on C satisfying another algebraic equation 

£'ix'iz),y'iz)) = (4-114) 

compared to 

£{xiz),yiz)) = 0. (4-115) 

proof: 

The proof rehes on the simple observation that 



^x-N j ydx = VLy-N j xdy = D,. (4-116) 

□ 

Corollary 4.1 If x has q poles and y has p poles, the partition function Z{t) is t- 
function of the multi- component p + q-KP hierarchy since it satisfies the Hirota equa- 
tions: 



E 



Z{t, + [z.']) Z{t', - [z.']) Y 
; — ^ m e 



Z(t. + [.r^])Z(t--[^r-]) ,"?^^-''"^'"'^^ 
Z{U) Z{Q 

(4-117) 

4.8 Singular limits 

The Fg's and oon^'?, can be computed for any regular spectral curve, i.e. as long as the 
branchpoints are simple. When the spectral curve is singular, the F^'s are not defined. 

Nevertheless, consider a 1-parameter family of spectral curves £{t), such that £^(tc) 
is singular, we prove below, that Fg{t) diverges as t — tc, in the following form: 

F,(t)~(t-g^'~'^^'^F,. (4-118) 

The goal of this section is to prove this divergent behavior, and compute the exponent 
IX and the prefactor Fg. These aymptotics are very important in many applications in 
mathematics and physics, for instance Witten's conjecture relates the asymptotics of 
large discrete surfaces, to integrals over moduli spaces of continuous Riemann surfaces. 
Asymptotic formulae play also a key role in the universal limits of random matrix 
eigenvalues statistics, or in the study of universality in the statistics of non-intersecting 
Brownian motions (see section [6]) . 
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In the context of matrix models quantum gravity, the prefactor Fg is called the 
"double scaling limit" of Fg, and the exponent 2 — 2/i is called 

2 — 2fi = 7strmg = " string susceptibility exponent" . (4-119) 

It is such that Fq formally diverges with the exponent — 7string: 



d^F 



° ~ (t-^c)"^''"""- (4-120) 



4.8.1 Blow up of a spectral curve 

Consider a one parameter family of spectral curves S{t) = {C{t),x{z,t),y{z,t)), such 
that S{t) is regular in an interval ]tc,^o]- For the moment we do not assume that S{tc) 
is singular, i.e. it may be either regular or singular. In a small vicinity of tc, we can, 
to leading orders, parameterize £{t) in terms of £ = C{tc). 

Moreover, let a be a branchpoint, and let us study the correlators ujn\zi, . . . ,Zn) 
in the vicinity of a. We choose a rescaled local coordinate ( in the vicinity of a. Let 
us write: 

z = a + {t-tcyC + oiit-tcy). (4-121) 

We want to compute the asymptotic behavior of uJn\a+ {t — tc)" (i, . . . , a + {t — tcy(n', t) 
in the limit t —>■ tc- 

First, let us study the behavior of x and y, by Taylor expansion. Let q be the first 
non-trivial power in the Taylor expansion of x, i.e.: 

x{a + (t - tcYC; t) = x{a; tc) + (t - tc)'''x{C) + o{{t - tc)^ (4-122) 

and similarly, there is an exponent p such that: 

yia +it- tcYC, t) = y{a- tc) + {t - tc^m + o{{t - tcf)- (4-123) 

This means that at t = tc the curve S{tc) behaves like y ~ y{a) + (x — x{a))'^/'^. It is 
regular if ^ = | and singular otherwise. 

The rescaled curve £ = {C, x, y) is called the blow up of the spectral curve near the 
branchpoint a, in the limit t — > tc- 

The choice of the exponent u, must be such that £^ is a regular spectral curve. 
We cannot give a general formula for z/, since it depends on the explicit choice of a 
1-parameter family of spectral curves S{t), and how it is parametrized. Also, here we 
consider only algebraic singularities of type y ~ x^^'^, but the method could certainly 
be extended to other types of singularities. 

Examples: 

• The following spectral curve arises in the enumeration of quadrangulated surfaces 
(see section [73]) : 

£ = P^ , x{z)=j(z + -) , 2/(^) = zl + ^ (4-124) 

\ z J '-)Z z-^ 
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where 



2 1 - VI - 12t4 



y = '-TT. -. (4-125) 

The branchpoints x'{a) = are a = ±1. Consider the branchpoint a = 1, and introduce 
an auxiliary scaling variable t: 

z = l + tC, (4-126) 

such that X and ?/ are independent of t, and let us study the vicinity of t — > tc = 0. 

We whish to study the behaviour of uJn\l+tCi, . . . , l + t(n) in the vicinity of t — » 0, 
i.e. the behaviour of ujn^ in the vicinity of the branchpoint a = 1. 

In the limit t ^ 0, we Taylor expand x and y: 

x{z) = 27 + jt^C + O(t'), (4-127) 

y(^) = — + t47^ + tC(- - 3t47^) + 0{t^). (4-128) 

7 7 

Notice that we have q = 2 and p = 1, which means that our curve is not singular at 
t = (which was expected since it is actually independent of t). 

The blow up is: 

^ (C) = 7 , m = - 3t47') C- (4-129) 

Up to a rescaling, it is the Airy spectral curve (see section [H] and example eq. fl2-4l) ). 

• In the previous example, at ^4 = ^, we have (1 — 8^47^) = 0, and thus, one needs 
to go further in the Taylor expansion. Let us now choose t = 1 — 12^4, tc = 0, and 
a = 1. We have in the limit t — > tc- 



X 



(1 + ^ tk) = 2v^ + ^ (C^ - 2) + o{V~t), (4-130) 



and, in fact, what we really need is the asymptotic behavior of y{z) — y{z), i.e.: 

y(l + ^(t,- t)-^c] - y( ^ —] = -— (C - 30 + o(th. (4-132) 

The blow up of the spectral curve in this limit is thus: 

5(C) = ^(C^-2) , ^(C) = -^(C'-3C). (4-133) 

Notice that it is proportional to the "pure gravity" spectral curve {p,q) = (3,2), see 
section [HI and see the first example in section [2. 1.11 It is the spectral curve which arises 
everytime we have a |/ ~ x^^"^ cusp singularity. 



46 



4.8.2 Asymptotics 



In order to study the asymptotics of the Un^^s, we need to study the asymptotics of 
the kernels B and K, in the hmit t tr- 



We have: 







z near a 


z far from a 


B{zo,z) ~ 


zq near a 


5(Co,C) 


o{{t-m 




Zq far from a 


o((t-t,n 


0(1) 



X (l + 0((t-g^)), (4-134) 



where -B(Co) C) is the Bergmann kernel of the blown up spectral curve {C, x, y): 



5(Co,C) 

Similarly, the kernel K behaves like: 
K{zo,z) 



(C-Co)^' 



(4-135) 





z near a 


z far from a 


zq near a 
Zq far from a 


0{{t- 


0(1) 
0(1) 



(4-136) 



where -ft'(Co; C) is the recursion kernel of the blown up spectral curve (£, a;, y): 



K{Co, C) 



Co-c Co + c; 2y(c)x'(c) 



(4-137) 



Therefore, we see that the leading contribution to LJn\l + 6(q, . . . , 1 + 6(n) is given 
by the terms where all residues are taken near a, and the leading contribution can be 
computed only in terms of B and K. By an easy recursion, on gets: 

Theorem 4.22 Singular limit of Un^ with 2 — 2g — n < and n > 1. 

If a is a branchpoint, the asymptotics of in the vicinity of t t^ and Zi — >■ a, 
are given by: 

ui^\a +{t- t^yc, ...,a + {t- t^rCn) ~ (t - (4.I38) 

where Un"* are the correlators of the blown up spectral curve S = {jC,x,y). 

In this theorem, the exponent u is unspecified, it depends on our choice of 1-parameter 
family of curves, i.e. it depends on each example. We recall that it must be chosen 
such that the blown up curve S is regular. We see examples in section [HI 

This theorem implies in particular, that all correlation functions in the vicinity of 
a regular branchpoint, are, to leading order, the same as the Airy process correlation 
functions, we recover the universals Airy law near regular branchpoints. This is related 
to the universal Tracy- Widom law (lllj . 

One may extend this theorem to uJq^^ = Fg^s: 
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Theorem 4.23 Singular limit of Fg. 

If at t = tc, the branchpoint a becomes singular, and no other branchpoint becomes 
singular, then the asymptotics of Fg,g > 2 in the vicinity of t ^ t^ are given by: 

Fgr^it- tj(2-2s)(P+g)- + - tj(2-29)(P+g)-) (4_i39) 

where Fg are the symplectic invariants of the blown up spectral curve £ = {C,x,y). 

In fact this theorem holds also in the case where the branchpoint a is not singular, 
but it becomes useless. Indeed, if a is not singular, the Blown up spectral curve is the 
Airy curve, and the F^'s of the Airy spectral curve vanish for all g > 1, and therefore 
all what the theorem says in this case, is that Fg does not diverge as {t — tc)^'^~'^^^''^'^'^^'^ . 
And this is obvious since Fg is not divergent at all. 

Also, the condition that a is the only singular branchpoint, is not so necessary. 
In fact, if several branchpoints become singular simultaneously at tc, with the same 
exponents u and (p, g), then the asymptotics of Fg is the sum of contributions of 
most singular branchpoints. The most basic example is a symmetric spectral curve 
y{x) = y{—x), for which branchpoints come by pairs. The leading order of Fg then 
gets a factor 2: ~ 2 (t - Q{^-mp+i> + o{{t - tj(2-2<;){p+g)^)_ 

5 Application to matrix models 

The recursion relations defining the symplectic invariants and their correlation func- 
tions were originally found in the study of the one-Hermitian random matrix model 
[ITJ [32] where they appeared as the solution to the so-called loop equations. It is thus 
interesting to remind to which extent the symplectic invariants give a solution to the 
computation of the free energies and correlation functions' topological expansions in 
different matrix models. It is also interesting to emphasize the special properties of 
the spectral curves obtained from matrix models to remind that they represent only a 
particular subcase in the whole framework for symplectic invariants. 

5.1 1-matrix model 

The formal 1-matrix integral is defined as a formal power series in a variable t. 

Consider a polynomial V{x) (called "potential") of degree c? + 1 > 2 as well as its 
d stationary points ^i. 

Vz = l,...,c?, y'(e.) = (5-1) 
and the non-quadratic part of its Taylor expansion around these points 

5V,{x) = V{x) - V(6) - ^^{x - i,f. (5-2) 
Let us also consider a d-partition of N, i.e. a set of d integers satisfying 

d 

Y^n, = N. (5-3) 

i=l 
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Definition 5.1 Formal 1-hermitian matrix integral. 

The formal 1 -matrix integral is defined as a formal power series in a variable t: 



7 A " ^ Hrt.-, *J Hn.j — 1 \ 



^IMM — « ' 



i Yl det {Mi ^ 1„ . - 1„. ® Mj) '' 



i>j 



(5-4) 



k=l 



where each is a (well defined) polynomial moment of a Gaussian integral. 
It is denoted by: 



formal 



This last notation comes from tlie excliangc of tlie Taylor expansion of exp — Tr 5Vj(Mj) 

i 

and the integral. Once this commutation is performed, the integral obtained corre- 
sponds to the formal expansion of the integral J e~T^^'^(^^ dM around a saddle point 
M solution of V'{M) = 0, which we choose as: 

ni rii nj 

M = diag {'^i,...,^i,...,^i,...,^i,...,^d,---,^d\ ■ (5-6) 



The integers nj thus correspond to choosing the number of eigenvalues of the saddle 
matrix located at a particular solution of the saddle-point equation V'{^i) = 0. 

However, in general, the Taylor expansion and the integral do not commute, and 
the formal matrix integral is different from the usual convergent matrix integral: 

J formal J 

In fact, typically, convergent matrix integrals are obtained for y > 0, whereas formal 
matrix integrals have combinatorical interpretations for y < 0. 

Remctrk 5.1 The definition of the matrix integral does not depend only on the potential 
(i.e. the coefficients of this polynomial) but also on the filling fractions 

ej := , i = l,...,d , = (5-8) 

i=l 
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The formal logarithm of Zimm is also a formal power series in t of the form: 

oo 

lnZiMM = 5Zt'^ (5-9) 

fc=0 

and it can be seen from general properties of polynomial moments of Gaussian matrix 
integrals (observation first made by 't Hooft jllO| ). that each coefficient N~'^Ak is a 
polynomial in 1/iV^: 

(k) 

Ak = N^ J2 ^i'^N-^'- (5-10) 

9=0 

Therefore, we may collect together the coefficients of given powers of N, and define a 
formal power series: 

oo 

= ^2-29 ^^(9) (5.11) 
k=l 

We have: 
Theorem 5.1 

oo 

In ZiMM = ^(iV/t)2-23 Fg. (5-12) 

9=0 

This theorem is an equality between formal power series of t. This means the coefficients 
in the small t power series expansions of both sides are the same. And for a given power 
of t, the sum in the RHS is in fact a finite sum. 



5.1.1 Loop equations 

We may also define the following formal correlation functions: 

Wr^ix,, . . . , x„) = / Tr . . . Tr -^-^) (5-13) 

where the subscript c means the cumulant and the notation Tr stands for the 
formal series: 

Tr ^ V y Tr ~ (5-14) 

to be inserted in the integrand of eg. fl5-4p . Again, Wn{xi, . . . , x„) is defined as a formal 
power series in t, whose coefficients are polynomial moments of Gaussian integrals. 
Moreover, one may notice that the coefficient of is a polynomial in 1/A^, and is a 
rational fraction of with poles at the ^'^s. 

Remark 5.2 Each coefficient of Wn{xi, . . . rational fraction of the Xj 's with poles 

at the ^j's, and from eq. (j5-14p . one sees that simple poles can appear only when A; = in 
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eq. (|5-14p . i.e. terms independent of Mj. This implies that the cumulants Wn{xi, . . . ,rc„) can 
have no simple poles when n > 1, and for Wi, the only residue is: 



and when n > 1: 



Res Wi{x)dx = Ui (5-15) 

Res Wn{xi^X2, ■ ■ ■ Xn)dxi = (5-16) 
Again those two equalities are equalities between the coefficients of formal series of t. 

We may collect together coefficients with the same power of A^. That allows to 
write: 



Ty„(xi, . . . , x„) = ^ ( - ) W^^\x,, . . . , x„) (5-17) 



where each Wn'^ is a formal power series in t, whose coefficients are rational fractions 
of Xi, . . . , a;„ with poles at the ^[s. Eq. fl5-17p is an equality of formal power series of t. 
For further convenience we also define in a similar manner: 

\ Xi — M X2 — M Xn — M / ^ 

which is a polynomial in the variable Xi. Again, we may collect together coefficients 
with the same power of A^~^, and define Pn '^ such that: 

Pn{xi;x2,...,Xn) = ^i—j Pjf\xi; X2, ■ ■ ■ , Xn) ■ (5-19) 

Then, the Schwinger-Dyson equations imply: 
Theorem 5.2 We have the loop equations, ^n,g: 



h=0 ICJ 



d M-t'(a-,J\fe})-tyi'"(J) 



= rix)wi%{x,J)-Pii,ix;J) (5-20) 
where J = {xi, . . . , a;„}. 
proof: 

This theorem can be proved by integrating by parts the gaussian integrals for each 
power of t, see |17]. It is called Schwinger-Dyson equations, or loop equations, or 
sometimes Ward identities... cf |96l [80l [36l [38] . □ 

Remark also that this theorem corresponds to the global Virasoro constraints of 
theorem 14.161 

Loop equations were initially used to find the topoloical expansion of the 1-matrix 
model in the special case of 1-cut [11] and also to the first orders 2-cuts [5]. 
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5.1.2 Spectral curve 

For n = and g = the loop equation eq. fl5-20p reduces to an algebraic equation for 
wl^\x) sometimes known as the master loop equation: 

{Wi°\x)f = V'{x) Wi''\x) - (5-21) 

where pI^\x) is a polynomial of x of degree at most d — 1. The one point function is 
thus given by 

= m,yZ|)!.p<«)M. (5.22) 

We define: 



y = ^- = - F<°'(x), (5-23) 



and the master loop equation implies that the function y is solution of the algebraic 
equation (hyperelliptical) HiMM{x,y) = where 



HiMM{x,y) := y' - J + Pi'ix) (5-24) 

which is called the spectral curve associated to the one hermitian matrix model. 

From eq. (15-151) . we see that, if one chooses the ^j-cycle to be a circle, independent 
of t, around ^j, then we have (order by order in t): 

^ I W['\x)dx = ^ = , i = l,...,d. (5-25) 

This last equation gives d constraints, i.e. the same number as the coefficients of the 
polynomial Pf'' (which is of degree d — 1), therefore it determines Pf'\x). Just by 
looking at the first terms in the small t expansion, one has: 

V'ix) 



Pi'\x) = J2e.^ + 0{n (5-26) 

i=i ^* 

The data rii/N are thus equivalent to the data of P^^\ Notice that = trii/N is of 
order 0{t) in the small t expansion. 

In case some e^'s are vanishing, we define g + 1 = the number of non- vanishing e^'s, 
and we assume that ei, . . . , e^+i are non- vanishing and 6^+2 , ■ ■ ■ are zero. 
Order by order in t we have: 



y= H {x-^^-A{t)) 

i=g+2 



9+1 



A lliix - ^^ - B,{t)y - AQ{t)) (5-27) 

\ i=l 



where Ai{t), Bi{t), Ci{t) are formal power series of t. To the first orders: 



4 ^ £, 
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^-^Ei:^ -^^-Oi,, (5-29) 

C. = :^ + 0(t^). (5-30) 

Let us study the specificities of tlie 1-matrix model spectral curve {Cimm, x,y). 
Genus of £imm 

The Riemann surface Cimm has genus g lower than d — 1: 

g<d-l. (5-31) 



Sheeted structure 

The polynomial HiMM{x,y) has degree 2 in y. This means that the embedding of 
'^iMM is composed by 2 copies of the Riemann sphere, called sheets, glued hj g + 1 
cuts so that the resulting Riemann surface Cimm has genus g. Each copy of the 
Riemann sphere corresponds to one particular branch of the solutions of the equation 
HiMM{x,y) = 0. Since there are only two sheets in involution, this spectral curve is 
said to be hyperelliptic. It also means that the application z — > 2 is globally defined 
since it is the map which exchanges both sheets: 

y{z) = -y{z). (5-32) 



Pole structure 

The function x{z) on the Riemann surface Cimm has two simple poles (call them a+ 
and a_), one in each sheet. Near a±, y{z) behaves like: 

y{z) ~ ±lv'{x{z))T^^+0{l/x{zn (5-33) 



5.1.3 The 2-point function 

For n = 1 and (7 = 0, the loop equation eq. (15-201) reads: 



d Wf\x)-Wf\xx) 



dx\ X — X\ 

i.e.: 



{y\x) - 2Wi^\x))wi''\x,x^) - Pi°^(x;xi) (5-34) 



9 , p(0). 



1 _d_ V'{x)-V'{x,)+2y{x^) p(0)/ . > 
-'- -|- ^ '^^^ x—Xi 2 V 1 ly 



2{x — Xi)^ 2y{x 
(5 - 35) 
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This equation shows that W2^\x,Xi) is a meromorphic function of z and Zi on the 
spectral curve Cimm- It is a multivalued function of x,Xi, but it is a monovalued 
function in the variables z, z\. Therefore let us write: 

Wf\x{z),x{zx))dx{z)dx[zx) = Gjf[z,z-v) (5-36) 

where ujf^i^z^Zx) is a meromorphic form of z and Z\. It is clear from the first line of 
eq. fl5-35p . that [0^2'^ {^z, z\) has no pole at z = Z\, but it has a pole oX. z = z\. Moreover, 
since dx{z) / y{z) has no pole at branchpoints, we see that uj^\z^ zi) has no pole when 
z approaches a branchpoint. By looking at the behavior at large x, we see, from 
degP^^ < d — 2, that uj^\z, zi) has no pole at the two infinities a±. 

From eq. (15-351) . it may seem that u)^2\^j^i) could have simple poles at the zeroes 
of y{z), but the residues are computed by eq. (15-161) and they vanish. Thus, we see that 
the only possible pole of 0)2°'' (-^5 ^1) can be at z = Zi. 

Then, notice that the second line of eq. fl5-35p is the sum of a term which is even 
under z —>■ z, and a term which is odd under z z. Since the sum of those two terms 
must have no pole at 2; = 21, we see that the pole at 2; = ^1 must be twice the pole of 
the even part. Therefore we find that i^^\z, zi) has a double pole at z = zi, with no 
residue, and no other pole. Moreover, eq. (15-161) implies that on every ^-cycle we have: 

u;^^\z,zi) = 0. (5-37) 



The only meromorphic differential having all those properties is the Bergmann kernel: 

(z,zO = -B{z,z,) = B{z,z,) - ^^(^)^^(^i) (5.38) 

{x{z) — X{Zi)Y 

(we choose k = 0). 

5.1.4 Higher correlators 

Similarly to what we just did with W2 \ we are going to compute every Wn^ and relate 
it to the symplectic invariants of the curve y{x). 

First, notice that the loop equations eq. (l5-20p imply recursively, that each Wn'' is 
in fact a meromorphic function on the spectral curve, and thus we prefer to rewrite: 

Definition 5.2 

L0\^\zi, . . . , Zn) = Wjf\x{Zi),...,x{Zn))dx{Zi) ...dx{Zn) 

dx{zi)dx{z2) . . 

{x{zr) - x{z2)r 

Indeed, correlation functions wlf'^ are multivalued functions of the complex variable 
X whereas the Un^ are monovalued n-forms on the spectral curve C Somehow, these 
forms are built to choose one particular branch of solution of the master loop equation 
eq.dEZH). 
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Structure of the form Un^: 

One clearly sees from loop equation eg. (15-201) . that ujI^I-^{z, zi, . . . , Zn) can possibly 
have poles only at branchpoints, or at coinciding points x{z) = x{zj), i.e. at z = Zj or 
z = Zj, or also at the zeroes of y{z). 

From the degrees one can see that there is no poles at the infinities 

a±. Also, there is manifestly no pole at 2; = Zj, and one may notice that cj^+i is an 
odd function when z ^ z, thus it can also have no pole when z = Zj. 

The zeroes of y{z) are either the branchpoints, or double points, which are, order 
by order in t, of the form (see eg. (15-2 7p ): 

^, + Ai{t)+Oit^) , Ai{t)=Oit). (5-40) 

Let us assume by recursion on 2g + n, that each Un^ has no poles at those double 
points. This is true for ujf'^ and uj^\ Assume that it is true for every tu^^ ■* with 
2g' + n' < 2g + n, let us prove it for 2g + n + 1. From the loop eguation eg. (l5-20p . one 
sees that ljI^I-^^{z, zi, . . . ,Zn) could have at most a simple pole at such double points. 
But because of eg. (15-161) and eg. (15-151) . the residue must vanish, and thus there is no 
pole. 

Therefore we obtain the following structure for the Un^^s: 

Lemma 5.1 Un'^ can have poles only at branchpoints when 2g + n > 3. 

Moreover, we see from eq. Ii5-16\) and eq. Ii5-15\) . that if 2g + n > 3 we have: 

(f ui'\xi,...,Xn) = 0. (5-41) 

5.1.5 Symplectic invariants 

Let: 

dS,,,,,{z) = r B{z,z') (5-42) 

J Z2 

be the 3rd kind differential in z, having a simple pole at z = zi with residue -|-1, and a 
simple pole of residue —1 at z = Z2, and no other poles, and normalized on ^-cycles: 



dS,,,,, = 0. (5-43) 



The fact that it has only simple poles with residue -|-1 at 2:1 allows to write the Cauchy 
formula on the spectral curve (o being an arbitrary base point on Cimm)'- 

ujnliiz, zi,...,Zn) = Res dSz',o{z) uj'fliiz' , Zi, . . . , z„). (5-44) 

On the other hand, the differential form ujI^I-^^{z, Zi, . . . , Zn) has poles only at the branch 
points z ^ ai, and thus the Riemann bilinear identity tells us that: 

Res dS,'^o{z) ^^nliiz', zi,. . . ,Zn) + y^^ Res dS^'^o{z) ^1, • • • , Zn) 
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V (b B{Z, Z) i Uj'fliiz', Zi,..., Zn) 



(5 - 45) 

Due to eq. fl5-43p and eq. fl5-4ip . the right hand side vanishes and thus: 

Res dSz',o{z) tol^liiz', zi, . . . , z„) = - Res dSz',o{z) ^^nli{z , zi, . . . , Zn). (5-46) 

■> ' — * y ' * y ' — * n - 



z'^ai 



Finally, one can plug in the loop equation eq. fl5-20p and remind that the polynomial 
has no pole at the branch points, and thus we find: 

^1, • • • 5 2;„) 
V Res K{z, z') z', zl 

' z'^ai L 



, . . . , A^nJ 



h=0 lc{zi,...,z„} 

(5 - 47) 



where -^'(-2, -z') is the kernel: 



Kiz z') = ^^-"'-'(^^ (5-4^ 



In other words: 



Theorem 5.3 ujn\zi, . . . , z„) are t/ie correlators and symplectic invariants 

of the spectral curve Simm of equation HiMM{x,y) = 0. 



To get this theorem, we also recover the F^'s from theorem I4.3[ or just by homo- 
geneity, see [32] for details. For the 1-matrix model, Fq has been known from the origin 
of random matrices, Fi was first found in [11] for the 1-cut case, and in [31] for the 
multicut case. The other F^'s were first found in |32j . 



5.2 2-matrix model 

The method of loop equations can also be used to solve the formal 2-matrix model. 
One of the main applications and reasons for introducing the 2-matrix model, was the 
problem of counting Ising model configurations on random discrete surfaces, or in other 
words bi-colored maps (see section [7]), it was first introduced and solved by V.Kazakov 
|81] . It corresponds to a formal 2-matrix integral. It can be rephrased in terms of 
symplectic invariants too. 

For this purpose, one generalizes the notion of formal matrix integral to integrals 
over two normal matrices. 
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Definition 5.3 Let N be an integer and V\ and V2 two polynomial potentials: 

di + l d2 + l 7 

^i(^) = -E|^' ' V,{y) = -J2jy'- (5-49) 

k=2 k=2 

Let d = did2, and let n he a d-partition of N: 

d 

:= {ni,n2, . . . ,nd} such that ''^^Ui = N. (5-50) 



n 



1=1 



Let {(^i, ''7i)}^li^ be the d = did2 solutions of the system of equations 

r K(e.) = V^ 
1 v^ivi) = 6 ■ 



(5-51) 



One defines the non-quadratic part of the Taylor expansions of the potentials around 
these saddle points 

6V,,{x) = V^{x) - V,{^,) - ^^{x - (5-52) 

and 

SV2M = V2iy) - V2iv^) - ^^{y - V^Y- (5-53) 
For all I, one defines the polynomial in t: 
Id 

k=l/2 

= JdMi... dMddMi . . . dMd{E^ Tr SVi,,{M,) + SV2,^{m)y 



n- 

i=l 

Yl det{Mi ® 1„^, - ® Mj-) n ^^*(^* ® " ® 



i>j i>j 



as 



a gaussian integral over hermitian matrices Mi and Mi of size ni x rii. 
The formal 2-matrix model partition function is then defined as a formal power 

t (cfmw- 



series m ^ ^ ^ , ^ 

00 2k 



Z2MM ■■= E^'(E^'^'^)- (^-^^) 

fc=o i=o 

As in the 1-matrix model, one uses the notation 

Z2MM = I ^-fTrV,(M,)+V,iM,yM,M, ^5_5q) 
J formal 
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One is also interested in the formal logarithm of the partition function: the free 
energy 

-^2MM := In 2'2MM (5-57) 
which has a topological expansion (due again to 't Hooft's observation [IIU] ) 

/'N \ ^"^^ 

where each Fg is a formal power series of t. 
5.2.1 Loop equations and spectral curve 

As in the 1-matrix model (cf eq. fl5-14p ). one also defines correlation functions by 

W,,i{x,, . . . , Xfc, 1/1, . . . , yi) := ( n Tr J] Tr -^-jr-) (5-59) 



denoted as the non-mixed correlation functiona^°l. These correlation functions also 
admit a topological expansion 

Wk^i{xi,...,Xk,yi,...,yi) = ( y ) wlfjixi, ■ ■ ■ ,Xk,yi, ■ ■ ■ ,yi). (5-60) 

One also needs the polynomials in xi and y 



\ V xi- Ml y -M2 J fj^ Xi - Ml 

(5-61) 



as well as 



^/t^' (^MM^)nTi. 1 \ (5.62) 

\ \xi- Ml y- M2 ) ^J-^ Xi- Ml I ^ 



which are polynomials in y only. 



5.2.2 Loop equations 



Loop equations proceed from integration by parts in the formal matrix integral (i.e. 
integration by parts in each gaussian integral for each power of t), or by writing invari- 
ance under changes of variables, i.e. they are Schwinger-Dyson equations. 



-'^'^There exists more general correlation functions mixing the two types of matrices Mi and M2 
inside the same trace, for example < Tr {AI1M2) >, but their study is too far from the main topic of 
this review to be treated here. It is studied in l58l [59l. 
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The loop equations for the 2-matrix model were first studied by M. Staudacher 



jl08j . and then written in a more concise form in [IHl [51]. The loop equations for the 
2-matrix model are (where J = {xi, . . . , 

y (y - Vl{x)) Un+i{x, y; J) + Un+2{X, v, X, J) 
+ Wi+\iio{x, I) Ui+n~\i\{x, y, J 1 1) 

ICJ 

^ d Un{x,y;J/{xj})-Unixj,y;J/{xj}) 

j = l J J 

N N'^ 
= --P„+i(x,|/; J) + — . (5-63) 

And then, identifying the coefficients of polynomials of for each power of t, we 

have: 

_ Vl{x) + W['i (x)) Uf^ (x, y) = (V^iy) - x) <J (x) - (x, y) + 1 (5-64) 
and for {g,n) ^ (0,0): 

{y - Vi{x) + W^i{x)) uifUx,y; J) + wii,^,{x, J) ui'\x,y) 
+Ui%~,'\x, y;x,J) + j2i2 <i/|,o(^' ^) utn\{x, y; J/ 1) 

h=0 ICJ 



= -~Pii,ix,y;J). (5-65) 
where J2h Yl'i means that we exclude the terms {h, I) = (0, 0) and {g, J). 

5.2.3 Spectral curve 

Consider the first loop equation: 

_ v;{x) + W^S(x)) Ui'\x,y) = (V^iy) - x) wi'i{x) - Pf\x,y) + 1. (5-66) 
It is valid for any x and y, and in particular we may choose: 

y = y{x) = Vi{x)-W[%). (5-67) 

Since uf\x,y) is a polynomial in y, it cannot have a pole at this value oi y = y{x), 
and thus, for this value of y = V({x) — Wi'^q{x), we have the algebraic equation: 

H2MM{x,y{x)) = (VM^)) - x) {V;{x) - y{x)) - Pi'\x,y{x)) + 1 = (5-68) 
known as the 2-matrix model spectral curve. 
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Let us study the specificities of tlie 2-matrix model spectral curve {C2MM, x, y) (see 
[82]). 

Genus of C2MM 

Tlie Riemann surface C2MM has genus g lower than did2 — 1: 



g < did2 - 1. (5-69) 



Sheeted structure 



The polynomial H2MM{x,y) has degree ^2 + 1 (resp. di + 1) in y (resp. x). This 
means that the embedding of £.2mm as a branch- covering of the x-plane (resp. y- 
plane) is composed by (i2 + 1 (resp. di + 1) copies of the Riemann sphere P^, called 
x-sheets (resp. ^/-sheets), glued by cuts, so that the resulting Riemann surface C2MM 
has genus g. Each copy of the Riemann sphere corresponds to one particular branch 
of the solutions of the equation H2MM{x,y) in y (resp. x). 

Since there are d2 + 1 x-sheets, this means that there are d2 + 1 points in C2MM 
corresponding to the same value of x. We write: 

x{z') =x{z) , i = 0,...,ci2 (5-70) 

We will take the convention that z*^ = z. 

Branchpoints are zeroes of dx, and they are also places where 2 sheets merge z ^ z'^ 
for some i. By convention, we call this other point z*, and thus z is one of the z*'s. In 
general, z is not globally defined but only defined locally around the branch points. 

Pole structure of the functions x{z) and y{z) 

The function x{z) (resp. y{z)) on the Riemann surface C2MM has two poles: one of 
degree 1 (resp. degree di) at a pole called oo^ and one of degree ^2 (resp. degree 1) 
at a pole called oOj,. It means that d2 x-sheets merge at oOy, and only one x-sheet 
contains oOx alone. 

Near 00^, a local parameter is 1/x, and we have: 

y{z) ~ V;{x{z))-^ + 0{l/x{zy). (5-71) 

And near ooj^, a local parameter is 1/y, and we have: 

x{z) ~ VMz))--^ + 0{l/y{zr). (5-72) 



According to section 14. 7[ the fact that we have two poles, means that the tau- 
function built from the symplectic invariants of this curve is the tau-function of the 
1 + 1-KP hierarch\^ 



Actually, one should take di and ^2 arbitrary large to obtain the 1+1-KP tau-function. The times 
of the hierarchy being given by the coefficients of the Laurent expansions of ydx and xdy around oo^, 
and OOy respectively. 
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5.2.4 Preliminaries to the solution of loop equations 

As in the preceding section, we promote the correlation functions to differential forms 
on the spectral curve to make them monovalued: 

Definition 5.4 

UjI^\zi, . . . ,Zn) = W^'^Q{x{zi),...,x{Zn))dx{zi) ...dx{Zn) 

dx{zi)dx{z2) , . 

{x{zi) - x{z,)r 

Exactly like in the 1-matrix formal model, the very definition of the model as a 
formal series in t, implies that each with 2(7 + n > 3 is a meromorphic n-form 
with poles only at the branchpoints (i.e. the zeroes of dx), and with vanishing ^-cycle 
integrals. The only exception is u!^\zi, Z2) which can have a pole only at zi = Z2, and 
which is found to be the Bergmann kernel (see [3^): 

J^\zuZ2) = B{zuZ2). (5-74) 

Before proving that the solution of loop equations for the ujn^^s are the symplectic 
invariant's correlators, we need a small lemma. Consider the "full spectral curve": 

Eix^y) = (Vlix) - y)iV^{y) - x) - ^Pi{x,y) + 1 (5-75) 

where Pi{x,y) = ^^{N/ty^'^^ Pl^\x,y). We also consider its descendants: 

En+iix,y; zi, . . . , z„) = 5„,o(^i'(a;) - y)iV^iy) - x) - ^P„+i(x, zi, . . . , Zn) + 5„,o- 

(5-76) 

We have: 
Lemma 5.2 



En+i{x{z),y;zi,. . . ,Zn} 

■•(n(y-vrw.')) + >'^_ 

j=0 ^ V y i 

n ^ 



(5 - 77) 



and 



. x(Zj) - Ml/ C,{xi,...,Xn} 



^Un+l{x{z), y) - Snfl{V2{y) - x{z)) 

-Wi"(n(l/-V;'(x(^^)) + ^Tr 
i=i ^ 



x{z^) - Ml 
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n ^ 

(5 - 78) 



. xiZj) — Mil c,{xx,...,x„} 

J=l J' 



wheretd2+i is the leading coefficient o/^'lz/); ^* ^'^s the preimages ofx{z) (see eg. ^5-7(j^ ), 
the subscript nieans that we take the connected part with respect to the Tr^^^j-y^^g^ 

terms, but not the Tr terms, and the inverted comas " (.)" mean that, every 

time one encounters a two-point function in the cumulant expansion, one replaces it 

^2,o(x, x') := / Tr Tr + , ^ (5-79) 

\ X — Ml X — Ml I [x — x Y 

For example formula eq. (15-77p . for n = 1, reads to the first subleading order in tjN: 
Pi'\x,y;zi) = i,,^iJ2w^iiz^,Zi) J] iv - vi^')) 

<J(z^zO<J(^') n iy-y{zn) 

d2 d2 

£ <J(^^z^Zl) n (y-yiz^)) 

jj^k=0 ij^j,k,i=0 

d2 d2 

+u,^i J2 w^^{z^,zi)w^i{z\z') n iy-yi^'))- 

j^k^l=0 ij^j,k,l,i=0 

(5 - 80) 

Formula eg. (15-781) would be almost the same, but with the indices k,l > 1 instead 
of > 0. 

This lemma was proved in [23] and relies on the fact that the loop equation eq. (15-651) 
has a unique solution admitting a topological expansion. The way to prove this lemma 
mostly follows from Lagrange interpolation formula for polynomials, as well as Cauchy 
residue formula on C2MM- 

Since En+i{x, y; zi, . . . , Zn) is a polynomial of y of degree > d2 + I, hj expanding 
En+i{x, y; Zi, . . . , Zn) in powers of y, this lemma gives d2 + I equations. In particular, 
the term in y'^^ gives (if 2g + n > 2): 

d2 

Y,u^^^li,oi^\zi,...,Zn) = 0. (5-81) 

1=0 

And the term in y'^^-i gjyes a bilinear equation in the correlation functions: 

i^j h=0 ICJ 



^Remark that this notation reminds the notation "det" in th l4.1l 
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= J2 2/(^>?ii,o(^^ J) + y{z')uji'l,o{z\ J) - f!f\x{z), J) dx{zf (5-82) 

where fn\x{z), J) is a rational function of x{z) with no pole when z approaches a 
branchpoint. 

5.2.5 Solution of loop equations and symplectic invariants 

Let us write Cauchy formula, exactly like for the 1-matrix model eg. fl5-44p : 

t^nliizo, J) = - Res dS^fl{zQ) uj^J!li{z, J) 
= Res dS^fl{zo) Jnli{z, J) 

(5 - 83) 

where we have moved the integration contours using Riemann bilinear identity like for 
the 1-matrix model eq. fl5-45p . 

Now, notice that near a branchpoint Oj, uJ^+i o(-2> ^) ^i+i o(^' J) have a pole at 
z = ai and all the uj^+i o(^*' J) such that z^ ^ z,z have no pole at 2; — > aj. According 
to eg. (15-811) . we have: 

4'+i,o(^, J) + ^nlxA^. J) = regular (5-84) 
and from eq. (l5-82p . we have near af. 



9 ' 



h=0 ICJ 

+regular. (5-85) 
Inserting this last equation into the Cauchy formula eq. (l5-83p . we find: 



u;i5n(^o, J) = y~] Res K{zo, . 

< * 7. — >r?.,- 



9 I 



(9-1) 

l^? ^, ^; -r 

/i=0 /CJ 



(5-86) 

where 

2(y(2;) - 2/(2;)) rfx(z) 

This gives the theorem: 

Theorem 5.4 co'i^''(2;i, . . . , z„) and Fg are the correlators and symplectic invariants of 
the spectral curve {C2MM,x,y). 

Here we have only briefly sketched the proof of [33], and we refer the reader to 
details there, in particular for finding the Fg^s. 

Fq was found for example in [881 120] , Fi was found in [HI HOj [50] , and the other 
Fg's in [33]. 
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Remark 5.3 The correlation functions ojn which we consider here, are expectation values 
of traces of only the matrix Mi (in terms of combinatorics of maps of section [71 they are 
generating functions for bicolored maps, whose boundaries are of color 1 only). There exists 
a generalization of symplectic invariants for all other possible expectation values, i.e. all 
possible boundary conditions, but this is largely outside of the scope of this review. We refer 
the reader to [59j for further details. 

An obvious remark, is that color 1 and 2, i.e. functions x and y play similar roles. We 
can obtain generating functions for bicolored maps, whose boundaries are of color 2 only, by 
just exchanging the roles of x and y, i.e. by computing residues at the zeroes of dy. 

In particular we may compute generating functions for bicolored maps with no boundaries 
(i.e. the Fg^s), with either x or y. In other words the -F^'s are unchanged if we exchange x 
and y. This is a special case of the symplectic invariance property Fg{£,,x,y) = Fg(jC,y,x). 

In fact, the general proof of symplectic invariance consists in defining some mixed gener- 
ating functions for bicolored maps, whose boundaries are bicolored, it was done in [61j . 



5.3 Chain of matrices in an external field 

Another matrix model which can be solved with the same technics is the chain of 
matrices matrix model. 

Consider the model of an arbitrary long open chain of matrices in an external field, 

which includes the one and two matrix models as particular cases. 

dk+i 

Consider m potentials \4(a;) = — ^^^x-' , k = 1, . . . ,m. The formal chain of 

i=2 

matrices matrix integral is: 

Cm m \ 

.enam , ^ ^==1 ''=^ ' dMi . . . dM^ = e ^ 

J formal 

(5-88) 

where the integral is a formal integral in the sense of the preceding sections and M^+i 
is a constant given diagonal matrix Mm+i = A with s distinct eigenvalues Aj with 
multiplicities If. 



(5-89) 



with = N. 

i 

Note also that we may choose Cm,m+i = 1 since it can be reabsorbed as a rescaling 
of A. 

Once again, in the definition of the formal integral, one has to choose around which 
saddle point one expands. Saddle points are solutions of: 

\fk = 1, . . . ,m, Vl{ik) = Ck-i,kik-i + Ck^k+iik+i , 3j, C.m+1 = Aj. (5-90) 
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This system is an algebraic equation with D = sdid2 ■ ■ - dm solutions. 

Therefore the choice of a saddle point is encoded in the choice of a set of filling 
fractions 



e, = T- 



N 



(5-91) 



for i = 1, . . . , D with D = did2 ■ ■ ■ dmS and rii arbitrary integers satisfying 

^ni = N. (5-92) 



5.3.1 Definition of the correlation functions 

The loop equations of the chain of matrices were derived in jSTl [52] . and they require 
the definition of several quantities, as follows: 

For convenience, we introduce in the sense of eq. fl5-14p : 



1 



We also consider the minimal polynomial of A, such that 5'(A) = 0, i.e.: 



(5-93) 



(5-94) 



i=l 



and we introduce the following polynomial in z 

e„).^^(£)^. (5.95) 

Cn,n+1 Z — 1\ 

We also define the polynomials fij{xi, . . . ,Xj) by fij = if j < i — 1, = 1, and 



fi,ji-^i: ■ ■ ■ ,Xj^ det 

\ 

if j > i. They satisfy the recursion 



Cj-ijXj 



(5-96) 



(^i—l,ifi,ji_-^i: • • • ; -^j) ^ (-^j ) j (^^j+l , • • • ; Xj) Cj Xj Xi^i fi+2(^'^i+2y • • • ; -^j)' 

(5-97) 



Let us then define the correlation functions and auxiliary functions: 

Wo{x) = (TrG'i(x)). 

For z = 2, . . . , m, we define: 



(5-98) 



Wi{Xi,Xi, ...,Xm,z)= Po\x„...,Xmfi,m{Xi, Xm) { Tr {Gi{Xi)Gi{Xi) . . . Gm{Xm)Q{z))) , 

(5-99) 
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which is a polynomial in variables Xj, . . . , Xm, z, but not in xi. And for z = 1, we define: 

Wi{xi,X2, ...,Xm,z) = Pol^,,,„,^„/i,m(a;i, . . .,Xm) (Tr {Gi{Xi)G2{x2) . . .Gm{xm)Q{z))) ■ 

(5-100) 

which is a polynomial in all variables. 
We also define: 

(2^1 ; Xii ■ ■ ■ ^ Z] Xi) 

)(Tr (GiK)) Tr {G^{x,)G,{x,) . . .GUxm)Q{z))),. 

(5 - 101) 

All these functions admit a topological expansion, for example: 

VTo = ^(Ar/t)i-29H^^(5) ^ Wi = J2iN/ty~^'Wi'^ ' VT, = ^(Ar/t)i-2^'PV.(^) 

(5-102) 

and 

Wi.^ = ^(iV/t)-2wg\ (5-103) 
9 

5.3.2 Loop equations and spectral curve 

In this model, the master loop equation reads [511 E2] : 

W2-l{Xi, . . .,Xm+i;xi) 

+ {Ci^2X2 - Vl{Xi) + —Wo{Xi))(^ — W2{Xi, . . .,Xm+l) " S{Xm+l)^ 

= -^Wi{Xi, . . . , Xm+l) + (^(3^1) - (^1,2X2) S{Xm+l) + 
t 

+ AT (^/(^i) ~ Ci-l,ia;i-l - Q,i+iXi+i) Wi+i{xi, Xi, . . . , Xm+l). 



(5 - 104) 



N 

i=2 



This equation is valid for any set of variables Xi, X2, . . . , x^+i, however, it can be 
simplified by choosing special values for those variables, in particular values for which 
the last terms in the RHS vanishes. For this purpose, one defines some Xj(xi,X2) as 
functions of the two first variables Xi and X2, as follows: 

Xi(xi,X2)=Xi , £2(3^1,3^2) = a;2, (5-105) 

and for z = 2, . . . , m: 

Ci,i+iXi+i(xi,X2) = V;'(Xj(Xi,X2)) - Ci_i,iXi_i(Xi,X2). (5-106) 

Choosing Xj = Xj(xi,X2), reduces the master loop equation to an equation in xi and 
X2: ^ 

W2;l(a:i,X2;Xi) + — (Ci,2X2 - F(xi)) f/(Xi,X2) = £'(Xi,X2) (5-107) 
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where 

I t N 

Y {x) = Vl{x) - —Wo{x) , U{xi,X2) = W2{xi,X2,Xs, . . . ,Xm+l) - Y^i^rn+l), 

(5-108) 

W2-l{xi,X2;Xi) = W2-l{xi,X2,Xs, . . .,Xm+i;xi) (5-109) 

and 

Eix,,X2) = -^W,{Xl,X2) + {V({X,)-CI,2X2)S{X,,X2) (5-110) 

with 

S{xi,X2) ^ S{Xm+l) , Wi{xi,X2) ^Wi{xi,X2,X3, . . . ,Xm+l)- (5-111) 

Notice that Wi{xi,X2), and thus E{xi,X2) is a polynomial in both Xi and X2- 
Finally, the leading order in the topological expansion gives 

,X2) = {C1,2X2 -y^°^(Xi)) U^°\xi,X2). (5-112) 

We may notice that this equation is more or less the same as in the 2-matrix model, 
and it is solved in the same way. 

Again, this equation is valid for any xi and X2, and if we choose X2 such that 
Ci,2X2 = Y^^\xi), we get: 

Hchain{xi,X2) E^°\xi, X2) = 0. (5-113) 
This algebraic equation is the spectral curve of our model. 
Study of the spectral curve 

The algebraic plane curve Hchain{xi,X2) = 0, can be parameterized by a variable 
z living on a compact Riemann surface C chain of some genus and two meromorphic 
functions x\{z) and X2{z) on it. Let us study it in greater details. 

Genus of Cchavn 

The Riemann surface C chain has genus g lower than D — s: 

g<D-s, (5-114) 

where D — sdi. . . dm- 
Sheeted structure 

The polynomial H2mm{xi, X2) has degree 1 + ^ (resp. di + ^^) in X2 (resp. Xi). This 
means that the embedding of Cchain is composed by H- ^ (resp. + ^^) copies of the 
Riemann sphere, called Xi-sheets (resp. X2-sheets), glued by cuts so that the resulting 
Riemann surface Cchain has genus g. Each copy of the Riemann sphere corresponds to 
one particular branch of the solutions of the equation Hchain{xi, X2) = in 2:2 (resp. 

Xi). 

Pole structure 
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In the preceding cases (1 and 2-matrix models), one was interested in the pole 
structure of only two functions x and y on the spectral curve. In the case of the chain of 
matrices, the problem is slightly richer since one can consider, not only the meromorphic 
functions Xi and X2, but also all the Xj(p) := Xi{xi{p), X2{p)) as meromorphic functions 
on Cchain- Their negative divisors are given by 

s 

[x k{p)] _ = -r koo - Sk'^Xi (5-115) 

1=1 

where 00 is the only point of Cchain where Xi has a simple pole, the Aj are the preimages 
of Aj under the map Xm+i{p)- 

Xm+l{k) = K (5-116) 
and the degrees and Sk are integers given by 

ri := 1 , Tfc := did2 • • • 4-i , Sm+i := , := 1 and Sk := 4+i4+2 ■■■dmS. 

(5-117) 

Note that the presence of an external matrix creates as many poles as the number 
of distinct eigenvalues of this external matrix M^+i = 

Remark 5.4 This matrix model has also a combinatorics interpretation in terms of counting 
colored surfaces. This interpretation is discussed in chapter [71 



5.3.3 Solution of the loop equations 



The loop equations have been solved in [52] by the same method as the 2-matrix model. 
It proceeds in three steps. One first shows that the loop equations eq. (15-1041) have a 
unique solution admitting a topological expansion. One then propose an Ansatz of 
solution and prove that it is indeed right. This gives 



Theorem 5.5 

E{x{z),y) 



(5-118) 



One finally develops this expression as a polynomial in y to get the bilinear relation 



h=0 



Z)ijJ 



(9-h)(- 



(5-119) 



where, as in the preceding section. 



Uj':^\zi, . . . , Zn) = Wlf\x{zi), . . . ,x{Zn))dx{zi) . . .dx{Zn) 

dx{zi)dx{z2) (Kion\ 
+^n,2(>gfl -T-J~-^ TTTTT?- (5-120) 



{x{Zx) - x{Z'2)f 



This allows to obtain the theorem: 



^■^The cases of matrix models without external field correspond to a totally degenerate external 
matrix A = cid with only 1-eigenvalue. There are thus two poles as in the 1 or 2 matrix models 
studied earlier. 
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Theorem 5.6 ujn yZi, . . . , z„) and Fg are the correlators and symplectic invariants of 
the spectral curve {Cchain,Ci,2Xi,X2) ■ 

Since €1^2X1 + 02,3X3 = V2{x2), we may use the symplectic invariance theorem 14.11 
and homogeneity theorem eg. (14-41) . which allows also to write: 

01,2X1, X2j 
~ Fgi^chain, —C2,3X3,X2) 

02,32:3, 2:2) 
aini 3^2,02,3X3) 
~ Fgi^chain, 02,32^2, X3). 

(5 - 121) 

And by an easy recursion, for any k = 1, . . . ,m: 

Fg — Fg(Cchain,Ck,k+lXk,Xk+l)= Fg(Cchain, Xk, Ck,k+1 Xk+l) ■ (5-122) 

In other words, the F^'s can be computed by choosing the spectral curve of any two 
consecutive x^'s. It means it doesnot depends on k, it doesnot depend on where we 
are in the chain. 



5.3.4 Matrix quantum mechanics 

Matrix quantum mechanics is the limit of an infinitely long chain of matrices m — » oo. 
This model is very useful in string theory [TO] . 

In this limit, the index k of the matrix Mk, becomes a continuous time variable 
t. The coefficients are scaled in a way such that the coupling term Tr {Mk — Mk+iY 
becomes a kinetic energy Tr {dM/dty. 

More explicitly, consider the chain of m matrices 



-AfTr 

dMi . . . dMm e 



m—l 



k=l k=l 



(5-123) 



and take the rj —>■ limit, and T = mr] of order 1. The index k becomes a time t = ki], 
and in the r] —>■ limit we have: 

Z = [ D[M{t)] e-^^ Ioiy(.M(t),t)+!^{dM/dt)^]dt (5-124) 



The spectral curve is characterized as before: 

find a compact Riemann surface C, and some time-dependent function x{z, t) ana- 
lytical in the variable z on some domain of C, which satisfy eg. (15-1061) which become 
Newton's equations of motion: 

fix{z,t) = -V'{x{z,t),t) (5-125) 
and such that the initial and final impulsions 

p(x(2;, 0), 0) = /ix(z, 0) , p{x{z,T),T) = ^x{z,T) (5-126) 
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are analytical outside some cuts. 

Therefore, theorem 15.61 gives Vt G [0,T]: 

\nZ = ^A^2-29 Fg{C,x{z,t),iJx{z,t)). (5-127) 
a 

In other words, the spectral curve is here the realtionship between impulsion and 
position for the classical equation of motion. Although the spectral curve depends on 
time t, the F^'s are independent of t. 

In the case where the potential V{x,t) = V{x) does not depend on time t, the 
equations of motion eq. (15-1251) can be integrated and give the energy conservation: 

E{z) = - x\z, t) + V{x{z, t)) (5-128) 

i.e. 

fxx{z,t) = ^2^{E{z) - V{x{z,t)) (5-129) 

and we have: 

In Z = J] A^2-23 ^^^^^ ^^^^ v^2/i(E(z) - V{x{z,t))). (5-130) 

9 

5.4 1-Matrix model in an external field 

As a special example of the chain of matrices above, let us consider the special case 
m = 1, i.e. 1-matrix model with an external field. 

The formal 1-matrix model in an external field A is defined as |118j : 



nil 



^M.ext 



(A)= f e-T^(nM)-^M)^M ^ A = diag A„...,A„...,A„...,A. 

J formal \ 

(5-131) 

where formal, as usual means that we Taylor expand near a critical value and then 
exchange the order of Taylor expansion and gaussian integral. A critical point is a 
matrix Mq solution of V'{Mq) = A. Let us assume that A has s distinct eigenvalues Aj 
of multiplicities rrij. Its minimal polynomial is: 

s 

S{y) = l[{y-Ai). (5-132) 



For each Aj, let = 1, . . . ,degV^' be the d = degl^' solutions of V'{$,ij) = Aj. A 

critical point Mq is characterized by s partitions of the mj's into at most d = deg V 
parts 



iTij = (5-133) 
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It is of the form: 

Mo = diag . (5-134) 

The parameters: 

= ^ (5-135) 

are called the filling fractions, and they parametrize which formal integral we are 
considering. 

The 1-matrix model in an external field is of course a special case of the chain of 
matrices described in the previous section 15.31 and therefore one finds that it has a 
topological expansion given by the symplectic invariants of a spectral curve: 

oo 

ln(ZM.ext(A)) = Y,{N/Tf-^^Fg{£M.e.t) (5-136) 

g=Q 



• • ■ ) - Ewtr--'-" ■ ■ ■ . (5-137) 



and 



x{zi)-M"' x{zn) - M / ^ 
However, let us make the results of section 15.31 a little bit more explicit in that case. 

5.4.1 Spectral curve 

The spectral curve SM.ext, obeys the equation: 

= £MU^,y) = (V'ix) - y) - 1 ^ ^M. (5.138) 

where Pj{x) is a polynomial of degree at most d — 1. Such a spectral curve is typically 
of genus g < sd — s. 

All the coefficients of all the P/s are fixed by the filling fractions requirement that 
(order by order in t): 

/ ydx = ei (5-139) 

where / = («, j) runs through the values 2 = 1, . . . , s, j = 1, . . . , c?, and Ai is a small 
circle around ^jj. We remind that the are not independent, only sd — s oi them are 
independent: 

E'« = ^- (5-140) 

The two functions x{z) and y{z) are characterized by the fact that x{z) has a simple 
pole at cxD, simple poles at some Aj such that y{\i) = Aj, and y{z) has a pole of degree 
d at oo. And we have: 

V'{x{z))-y{z) ~ -l^ + 0{l/x{zr) (5-141) 

z^oo X[Z) 
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and near Aj, the residues of x are such that: 



and the cycle integrals: 



5.4.2 Rational case 



Res xdy^ (5-142) 



1 f tn- ■ 
(h ydx = ei = — ^. (5-143) 



It is interesting to study the case of a rational spectral curve. The two rational functions 
x{z) and y{z) are of the form: 

J x{z) = 2; - Tr Q/(AV2_A) ('c, 1/1/1^ 

- 1 y(.) = g(.) 

where A is a diagonal matrix determined by: 

Q(A) = A (5-145) 
and (5 is a polynomial of degree d — deg V , determined by: 

V\x{z)) = Q{z) + ^ + 0(z-2). (5-146) 
5.5 Convergent matrix integrals 

So far, we have been discussing formal matrix integrals, which consist in exchanging 
integration and the small t Taylor series of e~T'^^(^). Formal matrix integrals always 
have a "topological expansion" of the type: 

InZ = p^^ty (5.147) 
9 

Now, let us consider a "convergent" matrix integral: 

Z^f dMe-T-^ym (5-148) 



where the integration domain Hj^{'y) is the set of normal matrices with eigenvalues on 
a path 7: 

Hj^i^y) = {M / M = UKU\ U G U{N) , A = diag(Ai, . . . , A^v) , A^ e 7} (5-149) 
equipped with the complex U (N) invariant measure 

dM = Y[{Ai - Ajf dU W dki (5-150) 

i>j i 
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where dU is the Haar measure on U{N), and dAi is the curvihne measure along 7. 

• For example Hn(R) = is the set of hermitian matrices. 

• For example H]\f{Si) = U{N) is the set of unitary matrices (5*1 is the unit circle). 
In eigenvalues, the convergent matrix integral is: 



1 f ^ 
— / dxi...dxN JJ(xi-Xj)^ JJe 

i>j i=l 



'-Vixi) 



(5-151) 



Imagine that \^ is a polynomial of degree d + 1 (the present discussion can be 
extended easily to a 2-matrix model, or a chain of matrices, and to all cases where V 
is a rational fraction [H]). There are d homologically independent paths on which the 
integral / e"^'-^-' dx is convergent, let us call 71, . . . , 7^, a basis of such paths (a choice 
of basis is not unique). See the figure for the example of a quartic potential {d = 3): 







The path 7 in the matrix integral eq. fl5-148p is a linear combination of such paths: 



7 



E 

i=l 



(5-152) 



and the convergent matrix integral eq. (l5-15ip can be written: 



^(7) = E 



'^1 ■ ■ - ^d 

ni\...nd\ J^^i 



N 



dxi...dxN ]^(2;j — a;j)^ Y\.^ 



'-V(x,) 



(5-153) 

This leads us to define the convergent matrix integral with fixed filling fractions nj, 
as: 



def 



^ni,...,na 



nil 



r / dXi . . . dXN TT( 



N 

>2 TTg-fvCx,) 



1=1 



(5-154) 
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(5-155) 



and thus we have: 

n-i+...+nci,=N 

This holds for any choice of basis 71, . . . , 7^. 

There is a conjectur^ that there exists a "good" basis of paths 71, . . . , 7^ such 
that ^rai,...,ra^ is a formal matrix integral (and thus it has a topological expansion in 
powers oit/N)\ 

The "good" paths 71, ... ,7^ can be seen as "steepest descent" paths, they should 
be such that the effective potential: 



Kfr(x) = V{x) - ^ (In (det(x - M))) 



(5-156) 



is such that along each 7j, the real part of the large leading order of Ves is decreasing, 
then constant and finally increasing, and at the same time, the imaginary part of V^g 
is constant, then increasing, then constant: 




If such paths exist, then we can write: 

Zni,...,na 

where e, are the filling fractions: 



t rij 
N 



(5-157) 



(5-158) 



and the coefficients Fg in the expansion, are the symplectic invariants Fg of the corre- 
sponding formal matrix integral. Since the F^'s are analytical functions of the filling 
fractions, we have, for any "background" filling fraction 77 = {rji, . . . ,?7d-i) 



k=0 



(5-159) 



"'^'^This conjecture is proved in some cases, and in particular proved for the 1-matrix model with 
arbitrary 7 and arbitrary polynomial V (the proof follows from M. Bertola's work [12]), but, at the 
time this article is being written, it is not proved in more general cases, for instance not proved for 
the general 2-matrix model. 
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where we assume tensorial notations and sums over indices. For simplicity, the unfa- 
mihar reader may assume d = 2, i.e. rj and e are scalar. 
Eg. (15-1 5 7p thus becomes: 

Z, = e 9 {n-Nri/t)NF^/t^l(n-Nrj/tf F^' 

g g k>2-2g 

^ e 9 ^{n-N,j/t)NF^/t^l{n-Nri/tf F^' 

I 9l,---,9l ki,...,ki I 

(5-160) 

where we have separated the terms with a positive power of from those with a 
negative power of N (^ means that we consider only terms with ki > and 2gi + ki — 
2 > 0). Then, writing: 

Ci = e2*"^% (5-161) 
we perform the sum over filling fractions in eq. (I5-155P , and we get [211 ES] : 



^(n-Nv/t)NF^/t^l{n-Nr,/t)'' F^ 



Z{l) ~ e s 

I gi,--;gi ki,...,ki I 

(5-162) 

That is, we find the non-perturbative partition function of section H75l (see [551 EH): 
Z{'~f) ~ e ^ 

Eh E E ""1;..., n^-if-'w eif^'cvfoA-f^) 

/ gi,--;gi ki,...,ki i 



(5 - 163) 
where 

0(M,.)(«,i"o) = 5^ e2-"-e2-("-^''/*+'^)"e^("-^''/*)'^o'. (5.164) 

This formula is expected to give the large N expansion of convergent matrix models. 
It was proved in several cases, but a general proof is still missing. 
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To the first orders, eq. fl5-163p reads: 

Z{^) ^ e^^° e^^ [e + ^ (e' F[ + ^^-^) +•••]. (5-165) 

We have to make several remarks: 
• background independence: 

Formula eq.f l5-163"ll is independent of the background rj. This was discussed in 
section 14. 5[ and it is related to the fact that the non-perturbative partition function is 
modular. 

This also implies that eq. fl5-163p cannot be a good large N asymptotic expansion 
for all values of r]. The conjecture is that one should choose 77 as a real minimum of 
ReFg. Notice that if rj is real, then 

ReF^' = 27r Imr > (5-166) 

where r is the Riemann matrix of periods of the spectral curve, and thus ReFo is a con- 
vex function of rj, within each cell of the moduli space of spectral curves corresponding 
to the potential V . Unfortunately, this moduli space is not very well known, and it is 
not known how to find such 77 in general. 

If such ?7 can be found we have §^ydx = 77 G M, and Fg = ydx, and thus, for 
any contour C we have: 

■ 4> ; 



Re ydx = 0. (5-167) 
A spectral curve with that property is called a "Boutroux" curve (See [19j). 



• Characteristics (/i, z^): 

We see here, that the characteristics (//, z/) of the 9-function, is related to the path 
7 chosen at the beginning to define the convergent integral. It is not fully understood 
how to associate a path to a characteristic and vice-versa. 



6 Non-intersecting Brownian motions 

6.1 Dyson motions and integrability: introduction 

Let us consider Brownian motions on the real line whose positions at time t are 
denoted by Xj(t) for i = 1, . . . , A^. Let us constrain them not to intersect and fix their 
starting and ending points: the particle i goes from at t = to 6j at t = 1: 

Xi{0) = at , Xi{l) = hi (6-1) 

with 

ai < a2 < ■ ■ ■ < ajq and hi <h2 < ■ ■ ■ < h^. (6-2) 

Once these parameters are fixed, one is interested in the statistic of these Brownian 
movers at a given time t G (0, 1). They are given by the correlation functions: 

RkiXu A2, . . . , Xk\t) = ^(jlj^ <5(A. - M)^ (6-3) 



76 



where we denote M = diag(a;i, . . . , xn)- 

These correlation functions can be written under a determinantal form HS 



Rkixi, ...,Xk\t)= det [HN,t{xu a^i)]fj=i (6-4) 

for some kernel H]^^t{xi, Xj) depending on time t. In particular, the density of Brownian 
movers at time t is given by 

Ri{x) = H^t{x,x). (6-5) 

Let us now consider a particular case where some of the starting and ending points 
merge in groups: 

ni 712 7^2 

(oi, . . . , otv) = (ai, .r. ^2, .r. ,a2\ . . . , h2,.^.,a2, • • • , hp, . T. , a^) (6-6) 



and 



with 



ni n2 n2 

. . . ,6^) = . . . ,/3i,/32, . . . ,/32, . . . ,/32, • • • ,/32, • • • • • • ,/5g) (6-7) 



J]n, = J]n, = iV. (6-8) 

i=l j=l 

It means that one considers p groups of rii particles starting from distinct points at 
at t = 0, merging for intermediate times and splitting into q groups of rii movers 
reaching /5j at t = 1. Remark [T] that the kernel reduces to the kernel of the p + q 
mutli-component KP integrable hierarchy: 

HN,t{x,x') = Hi^f{x,x') (6-9) 

since it satisfies the corresponding Hirota equation. 

In the following one studies the behavior of this phenomenon as the number of 
particles goes to infinity while the ratios 

= — and Cj = — (6-10) 

are kept fixed and finite. In this case, the Brownian movers form clouds which fill a 
connected region of the complex plane describing the space time. For example, for 
1 starting point and two ending points one typically gets a configuration of the type 
depicted in figure [D all movers leave the origin and begin to flee from one another. 
It creates a larger and larger segment of the space flUed by the Brownian movers. As 
the time grows, because the movers want to reach different points, they split into two 
groups heading towards these two end-points. 

As often in the study of such integrable system, the kernel exhibits universal be- 
haviors : in any point of the space time, one can rescale the kernel so that one obtains 
a universal kernel independent of the position of the considered point; typically, one 
recovers the Sine, Airy and Pearcy kernels. The purpose of this part is to emphasize 
the role played by the spectral curve and algebraic geometry in the study of these 
universality properties. 
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Figure 1: Example of one starting point at and two ending points at +1 and -1 half 
of the movers going to +1 and the other half to -1( the particles go from the left to 
the right). In the large N limit, the Brownian movers fill the sector of the space time 
delimited by the figure. One sees that before t — ^, they are distributed along a unique 
segment which splits into two disjoint segments for i > |. 
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6.2 Particles starting from one point: a matrix model repre- 
sentation 



Let us consider the particular case of one starting point p = 1. This means that all 
the Brownian movers start from the same point which we may assume to be located at 
the origin of the real axis. It was proved [46] that this case is equivalent to a matrix 
model. More precisely, the statistic of the Brownian movers at a given time t is the 
same as the statistic of the rescaled eigenvalues of an hermitian random matrix M of 
size N X N submitted to an external field A{t). This matrix model is given by the 
partition function 

Z (Ait)) = J dMe-"""' {MHA(t)M) (6.11) 

and 



=diag ^A,{t), . . . ,A,{t),A2{t), . . . ,A2{t), . . . ,Ag{t), . . . ,A^{t) | (6-12) 
with the time dependent elements 



2t 

A{t) = ^lj^^^f3,. (6-13) 

The limit of a large number of particles corresponds to the large matrix limit. 

This model can also be expressed in terms of the eigenvalues {xi, . . . ,xiy) of the 
random matrix M using the sell known HCIZ integral formual [TOl [71|. They are 
submitted to a probability measure 

dfiixi, . . . ,xn) = YldxiA{xye~''^'y^~^'''') . (6-14) 



1=1 



After the rescaling 



Xi^ Xi^/t{l-t), (6-15) 

they have the same statistic as Brownian movers starting from the origin. 
The correlation functions 

Rk{xi,X2,...,Xk) = ^JjTr(5(xi -M)^ , (6-16) 

have also a determinantal expression 

Rk{xi, ...,Xk) = det [HN,t{^i, , (6-17) 

in terms of a kernel HN^f As A^ goes to infinity, these eigenvalues merge into dense 
intervals of the real axis and the density Ri{x) is supported by a finite number of 
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segments [z2i~i,Z2i]- A classical result of the study of random matrices states that 
these segments correspond to the discontinuity of the resolvent W{x) = ^^^^ 

W+{x) - W-{x) = Ri{x) (6-18) 

for X e [Z2i-1, Z2i]. 

6.2.1 Gaussian matrix model in an external field 

Let us now quickly remind this matrix model's loop equations and large solution. 
This model is a special case of the matrix model in an external field studied in section 

m 

The resolvent W{x) is thus the solution of the algebraic equation eq.f l5-138"|) with 
V\x) = X, i.e.: 

q 

Eix, Y{x)) = Y{x) - ^ + E Y{x)-a.{t) = ° ^^'^^^ 

where 

Y{x) = W{x) - X. (6-20) 

This equation can be seen as the embedding of a Riemann surface C into CP^ x CP^. 

First of all, one can see that this spectral curve has always genus 0: it admits a 
simple rational parametrization 

<? 



y{z) = z 



^ N{z-a,(t)) _ (6-21) 



It is composed of g + 1 sheets, one of which contains the pole z = oo. This sheet is 
called the physical sheet. 

For a fixed number q of distinct eigenvalues of the external matrix, it may have 
up to q cuts linking the physical sheet to the q others. These cuts correspond to the 
discontinuity of the resolvent giving rise to the density of eigenvalues eq. (16-181) . The 
cuts [z2i-i, Z2i] are thus the support of the eigenvalues. Each Zi is solution of x'{zi) = 0. 
To be precise, using the notations of section [21 the density of eigenvalues on the cut i 
is given by 

p(x(p))=i/(p(°))-i/(p«) , for xip)e[z2^.i,Z2i] (6-22) 

where one labels the physical sheet by and the label i corresponds to the sheet linked 
to the physical one by the cut [z2i-i, Z2i\. 

Let us now study the evolution of the structure of this spectral curve as the time 
evolves from to 1 : one is particularly interested in the time evolution of the position 
of the branch points Zi[t). 
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Figure 2: For large times, the spectral curve is composed of g + 1 sheets linked by 
q real cuts [z2i-i,Z2i]. These cuts can be seen as the section of the support of the 
Brownian movers at fixed time. In the example depicted in figure [1], it corresponds to 
times greater than i. 



For a given time t, these branch points are the simple real roots of the equation 



Let us first consider large times close to 1. In this case, the eigenvalues ai{t) become 
large and are far apart from one another. Thus, the equation has 2k distinct real roots 
in z: the spectral curve has k distinct cuts [z2i-i,Z2i] (see figure [2]). 

Then, as the time decreases, the branch points come closer from one another and 
merge for some critical time before becoming complex conjugated with an increasing 
imaginary part. It means that two real cuts merge into one and an imaginary cut 
linking two non physical sheets appear (see figure [3]). 

Finally, as the time becomes small and approaches zero, all the branch point are 
coupled complex conjugated numbers except two of them: the ones with the smallest 
and the largest real parts: there is only one real cut left. 

This time evolution of the spectral curve has a simple interpretation in terms of 
statistic of the eigenvalues. Indeed, thanks to eq. (16-181) . the real cuts are the support of 
the random matrix eigenvalues (whereas imaginary cuts just follow from the interaction 
between the different groups of eigenvalues). In terms of non intersecting Brownian 
motions, these real cuts are the segments filled by the Brownian movers at a given time 
t: it is the constant time section of the region of space-time filled by the Brownian 
movers. 

The time evolution of the spectral curve can thus be interpreted as follows. For 



dyE{x,y)\ 



y=y{x) 



0: 



k 




(6-23) 
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Figure 3: For intermediate times some of the real cuts have merged and created imag- 
inary cuts hnking two non-physical sheets. In this example, Z2 and have collapsed 
and given rise to the imaginary cut [z, l] . 



times close to one, the movers form q groups lying on segments centered around the q 
end points and whose extremities are the branch points of the spectral curve. As the 
time decreases, the branch points come closer to each other and finally some of them 
merge, i.e. two of the disjoints segments supporting the Brownian movers merge into 
one and the sector filled by the Brownian movers exhibits a cusp. Then the different 
segments keep on merging as time decreases until they give a simply connected support 
for t ^ 0. This follows the intuition that all the particles leave in one group which 
step by step splits into smaller groups to end up with q groups reaching the end points 
at t = 1 (see figure [1] for q = 2). 

Remark 6.1 For the following, it is interesting to note that the critical times when two 
disjoint segments merge correspond to singular spectral curves in the sense of definition 12. 3i 
Indeed, at this time, the spectral curve has a double branch point at the location where the 
two simple branch points merge. 

6.2.2 Replica formula and spectral curve 

Let us now follow another approach, exact for finite N, exhibiting the role played by 
the spectral curve directly in the formulation of the kernel Hn. For this purpose, we 
sketch the derivation of a double integral representation of the kernel using the replica 
method developed in this context by Brezin and Hikami pll [251 [26l [27] . 
Let us first consider the "Fourier" transforms of the correlation functions 




(6-24) 
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and, in particular, one gets the Fourier transform of the two points correlation function: 



Z{A)N 



1 ^ r f ^ 



A{x) 



N 



-rf-Xj[^aj+itlSj^ai+it2Sj.a2 ) 



01,02=1 \i=i 
One can now integrate the variables xj by noting that 



(6-25) 



N 



2 



A{b)e ^=1 



(6-26) 



and using the expansion A(x) = Ili^jl^* ~ ^i)• 



^/2(^l,^2 



AT 

E ^ 

01,02=1 



-tlt2<5a]^,Q2^ 



X 



JJ^ (ai - am + itl{5l^ai - 5m,ai) + it2{5l,a2 " (^m.oa)) 



l</<m<Af 



n 



l<Z<m<Af 



(6 - 27) 

One can see that this can be written as a double contour integral 



f/2(tl,t2 



6 2 



tit 



1^2 



dudv ^m(tiu+t2v) {u - V + iti - it2){u - v) 



X 



(2i7r)2 



(« — !' + iti)(u — V — it2) 



X 



1 + 



^to 



(6 - 28) 



or 



e 2 



tit2 



dudv 
(2i^ 



tlt2 



X 



n 1+ 



zti 



u- at 



{u — V + iti){u — V — it2) 
it2 



V - ak 



(6 - 29) 



where the integration contours encircle all the eigenvalues ak and the pole v = u — it '^^' 
We can now go back to the correlation function 



/?2(A,/i) 



00 pCO 



dtidt 



00 J —00 



47r^ 



l^-iN(tiX+t2fJ.) 



U{h,t2 



(6-30) 



15 



see for instance [53] for more details around eq.(2-20) and eq.(4-40). 
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By first integrating on ti and t2 with the shifts ti ^ ti — iu and t2 ^ t2 — iu, we can 
show that: 

i?2(A, /i) = K^iX, A)ir^(/i, /x) - i^Ar(/i, A)irjv(A, /i) (6-31) 
where the kernel is defined by 

J 2lT J 2tTC ^^^\V - ttk J V - it 

where the integration contour for v goes around all the points and the and the 
integration for t is parallel to the real axis and avoids the v contour. Moreover, it 
is straightforwardly proven that any /c-point function can be written as the Fredholm 
determinant: 

i?fc(xi, . . . , Xfc) = det [KN{xi, Xj)]\ .^-^ . (6-33) 
By Wick rotating the integration variable t it, one gets 

-«(^-) - /l^/l^n (^) ^e-e^---) (e-34) 

where the integration contour for t is now parallel to the imaginary axis. One can then 
rewrite it under a more factorized form: 

i^iv(A,/i) = l^i ^,-NiSM-six,t))J_ (g.35) 
J Im J liTX V — t 

where 

,2 



S{x,y) = — - xy + ^e, ln(y - a^) (6-36) 

i=l 

with the "filling fractions" given by 

:= ^. (6-37) 

How to compute such an integral ? Let us use the saddle point method for both 
integrals. The saddle points of the first exponential are given by the y solutions of 

k 

dyS{x, y) = y-x + Y, = (6-38) 

i=i y 

which is nothing but the equation of the spectral curve eg. (16-191) ! 

In this setup, the spectral curve can thus be seen as the location of the saddle points 
of the action S{x, y) in this formulation of the kernel. Remark also that this formulation 
of the kernel is very similar to the formulation of th J4.4l in terms of symplectic invariants. 
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6.2.3 From the replica formula to the symplectic invariants formalism 



Let us now start from the expression of th J4.4l for the kernel. For this purpose, one has 
to compute the one form ydx. Using the global parameterization 



one gets 



x{z) = z + y. — ^ 

^ ' z—ai{t) 

y{z) = z 



7 / \ \ — ^ €i(iz 

ydx[zj = zdz + 



z - aiif)' 



Integrating by parts, one can see that 



ydx = — x{z)dy{z) + [x{z)y{z)] 



(6-39) 



(6-40) 



(6-41) 



Z2 



Moreover, the function y{z) can also be considered as a global coordinate on the spectral 
curve since it coincides with the z coordinate. The previous equation can hence be 
written 



yi 



ydx{y) 



y2 



x{y)dy + [x{y)y]l\ . 



(6-42) 



S/2 



On the other hand, the spectral curve has a particular form: it has only one x-sheet 
and takes the form 

H{x,y) = x-x{y) = Q (6-43) 

with the function 

Thus the action, i.e. the integral of the spectral curve wrt reads 



S{x,y) =xy- x{y)dy 



(6-45) 



and further 

S{,xi,yi) - S{x2,y2) 



xiyi - X2y2 

"Vl 



y2 
yi 



x{y)dy 

ydx{y) + yi [xi - x{yi)] + y2 [x2 - x{y2)] . 



(6 - 46) 



y2 



It is now possible to compare the kernel built from the symplectic invariants and the 
spectral curve 



1 N r,^''iKdx{y) 

-g ■'y(^2) 



Hn{xi,X2) = 

Xi - X2 

and the kernel following the replica formula 



1 + 



Hn{Xi,X2) 



dyidy2 



1 



MS{xi,yi)-s(x2,y2)] 



Vl - 2/2 
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(6-47) 



(6-48) 



Indeed, the later reads 

Hn{xi,X2) = j y"rfyirfi/2^^^(x(yi),x(i/2))e^^^^["^-"'(^^)l+^^t"^-"(^^^]^ (6-49) 

The saddle point equation directly then states the equality of both kernels in the large 
N limit if the integration contours are good steepest descent contours. 

Remark 6.2 Note that this particular representation of the kernel under the form of a 
double integral on the complex plane (or more precisely the Riemann sphere) is a direct 
consequence of the specific parameterization of the spectral curve x = a;(y). 



6.2.4 Critical behaviors and singular spectral curves 

Since both kernels coincide, one can use the singular limits derived in the symplectic 
invariant setup to obtain some critical universal behaviors of the exclusion process. 

• Universality on the Edge: the Airy kernel 

Let us study the statistic of Brownian movers around the edge of their limiting 
support. It means that one is interested in a point (x, t) of the space time such that 
(x, y{x)) is close to a branch point (xc, ?/c) = (x(zc), l/(^c)) of the spectral curve at time 
tc- Following the study of section H7g| let us rescale the global coordinate z = Zc + -^Z 

and one gets the blown up spectral curve 

x{z) = x{z,) + Ax'{z^) + -^x"{z,) + O (i) 

V '7 1 \ ■ (6-50) 

,(.)=,(..) + ^,'(..) + 0(^) 

Since the point Zc is a branch point, one has x'{zc) = 0, and the blown up curve reduces 
to 



(6-51) 



to leading order as ^ oo. Remark that the scaling is such that ydx{Z) = O (-^). 
This curve is the Airy curve described in the second example of section I2.1.1[ The- 
orem U]22] states that, in terms of the rescaled variable Z, i.e. the distance from the 
considered critical point, the kernel and the correlation functions reduce to the one of 
the Airy curve, independently of the position of the critical point. 

Note that this Airy curve has also the form x = x{y) where the function x{y) = 
n y^li v"^- The kernel can thus also be written under a double integral form 

HAiry{xuX2) = I i ^^^^^2 ^SMry(x,,y,)-SMry{xuy.) (6-52) 
/ / 1/1-1/2 

with SAiry{x,y) = xy - ^r^Y' ^-^^ ^^^y kernel. 
• Universality at the cusp: the Pearcy kernel 
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Let us finally consider a point where two groups of Brownian movers merge. This is 
obtained when two edges merge or, in the spectral curve formalism, when two branch 
points merge as time decreases. At this critical time the corresponding spectral 
curve is singular since the merging of two simple branch points gives rise to a double 
branch point. Let us be more specific: one considers a cusp at the position (xc,tc) in 
the space time. Let us blow up the space time around this point using the rescaling 



(6-53) 



where z is the global parameter of the spectral curve. The rational parameterization 
of the spectral curve thus reads 

x{z, t) = x{zc, tc) + ^d^x{zc, tc) + (d^)^ x{zc, tc) + (d^)^ x{zc, tc) 

+ ^^M2l^d,d,xiZc,tc)+0{j^) 



y{z) = y{zc) + ^y'izc) + O 



(6-54) 

Since the critical point is a double branch point, the blown up curve reduces to 



x{Z, T) = ^ (dz)^ x{zc, tc) + ayatZTdzdtx{zc, tc 



y{Z) = ^dMzc 



(6-55) 



which could be called the Pearcy curve. Indeed, since this curve also has the form 
X = x{y), the associated kernel has the double contour integral representation 

HpearcviXl,X2) = I / ^^^^^^ ^Sp,^^^,{x,,y,)-Sp^^^^y{x,,y,) .g_5gN 
/ / Z/l - 1/2 

With Spcarcy{x. y) = Xy - |1 - (^-^^) y^T . This is UOtiug but the 

Pearcy kernel once the right integration contour is found and the rescaling coefficients 
are fixed by: 



al {dzf x{zc,tc 
(dzVizc))^ 

and 

atdzdtx{zc,tc) 



6 (6-57) 



2d,y{zc) 

7 Enumeration of discrete surfaces or maps 



(6-58) 



The symplectic invariants provide a solution to Tutte's equations for counting dis- 
crete surfaces (also called maps), for arbitrary topologies. Indeed, it was found by 
Brezin-Itzykson-Parisi-Zuber |28], and further developed by [12], [351 EH], that gener- 
ating functions for discrete surfaces can be written as formal matrix models. Let us 
review how to enumerate various ensembles of discrete surfaces. 
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7.1 Introduction 

Definition 7.1 Let Mn^ be the set of connected orientable discrete surfaces of genus 

g obtained by gluing together polygonal faces, namely triangles, quadrangles, ... 
Uk k-angles, as well as n marked polygonal faces of perimeters /i, . . . , /„, each of the 
marked faces having one marked edge on its boundary. Let us call v the number of 
vertices of a discrete surface, and let M.n\v) be the set of discrete surfaces in M.n\ 
with V vertices. 

We require that unmarked faces have perimeter > 3, whereas marked faces are only 
required to have perimeter > 1. 

Marked faces are also called "boundaries". 

Notice that nothing in our definition prevents from gluing a side of a polygon, to 
another side of the same polygon. 



Theorem 7.1 M)f'{v) is a finite set. 
proof: 

Let e be the number of edges of a discrete surface in M^^^(v). The total number of 
half edges is: 

n 

'2e = J2jn, + J2h. (7-1) 

j>3 i=l 

The Euler characteristics is: 



1 1 " 

^ = 2-2^ = ^;-e + n + + J^(j-2)n,--^Z,. (7-2) 

j>3 j>3 i=l 

This implies: 



1 1 

^J2^j-'^)nj + -J2h^2g-2 + n + v (7-3) 

j>3 i=l 

and therefore the n^-'s are bounded, and the /j's are bounded. There is then a finite 
number of possible discrete surfaces having a finite number of faces, edges and vertices. 
□ 

In order to enumerate discrete surfaces, we define the generating functions: 
Definition 7.2 The generating function is the formal power series in t: 

Wjf\xi,...,Xn;t3,...,td;t) = —Sn,iSg,o 

oo 1 +n3{S) .na{S) n 

^Z^' #Aut(5) x^^iS) IniS) ii^ 

(7-4) 

Most often, we will write only the dependance in the Xi 's explicitly, and write: 

Wjf\x,,...,x^;ts,..., t,- 1) = Wjf\x,, ...,x^). (7-5) 
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The generating functions counting surfaces with marked faces of given perimeters 
li, . . . Jn are by definition: 

Tj\^^ ^ = (-1)" Res ... Res x[' ...x^;: Wlf\xi, . . . , x„) c/xi . . . dx^. (7-6) 

Notice that rooted discrete surfaces, i.e. surfaces with only 1 marked edge, have no 
non-trivial automorphisms, and thus #Aut = 1 when n = 1. 



7.2 Tutte's recursion equations 

Tutte's equations are recursions on the number of edges |112[ 1113] . If one erases the 
marked edge on the 1st marked face whose perimeter is /i-|-l, several mutually exclusive 
possibilities may occur: 

• the marked edge separates the marked face with some unmarked face (let us say 
a j-gon with j > 3), and removing that edge is equivalent to removing a j-gon 
(with weight tj). We thus get a discrete surface of genus g with the same number 
of boundaries, and the length of the first boundary is now /i + j — 1 . 

• the marked edge separates two distinct marked faces (face 1 and face m with 
2 < m < n, ), thus the marked edge of the first boundary is one of the Im 
edges of the m}^ boundary. We thus get a discrete surface of genus g with n — 1 
boundaries. The other n — 2 boundaries remain unchanged, and there is now one 
boundary of length li + 1^ — ^■ 

• the same marked face lies on both sides of the marked edge, therefore by removing 
it, we disconnect the boundary. Two cases can occur: either the discrete surface 
itself gets disconnected into two discrete surfaces of genus h and g — h, one having 
I J| -|- 1 boundaries of lengths j, J, where J is a subset oi K = {I2, . . . , /„}, and 
the other discrete surface having k — \ J\ boundaries of lengths Zi — 1 — j,K/ J, 
or the discrete surface remains connected because there was a handle connecting 
the two sides, and thus by removing the marked edge, we get a discrete surface 
of genus g — 1, with n + 1 boundaries of lengths j, h — j — 1, K. 

This procedure is (up to the symmetry factors) bijective, and all those possibilities 
correspond to the following recursive equation: 



d-i g n 

j=0 h=0 JCK m=2 
d 

i=3 
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This equation is illustrated as follows (where the 1st marked face is the "exterior face"): 




7.3 Loop equations 

Rewritten in terms of the wlf^^s, Tutte's equations eq.( ]7-7p read: 

Theorem 7.2 Loop equations. For any n and g, and L = {x2, ■ ■ ■ we have: 



h=0 JCL 

^yA_ L/{j}) - i^S(L) 



i=2 



30 ^ OC rj 



where 



V'{x) = X — tiX- 
i>3 



(7-8) 
(7-9) 



where Pn\xi, L) is a polynomial in Xi, of degree d — 3 (except P^^^ which is of degree 
d-2): 

d-l j-1 oo ,-j-{g) 

Pi^\xi,X2, ...,Xn) = -Y^ tj+i ^ X\ 



i=2 



i=0 l2,-,ln=l ^2 



^i^l^ (7-10) 
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proof: 

Indeed, if we expand both sides of eg. (17-81) in powers of Xi — > oo, and identify the 
coefficients on both side, we find that the negative powers of the Xj's give precisely the 
loop equations eg. (17-71) . whereas the coefficients of positive powers of Xi cancel due to 
the definition of Pn \ which is exactely the positive part of V {xijWn'' '■ 

Pi^)(Xi,X2,...,X„) = Pol {V'{Xi)Wi^\xi,X2,...,Xn)) (7-11) 

where Pol means that we keep only the polynomial part, i.e. the positive part of the 
Laurent series at xi — > oo. □ 

We see that Tutte's eguations eg. (17-80 are identical to the matrix model loop egua- 
tions eg. (15-201) . The solution is thus the same, and it is expressed in terms of symplectic 
invariants. One only has to find the corresponding spectral curve. 

7.3.1 Spectral curve and disc amplitude 

The spectral curve is given by the function W^\x). 

With n = 1 and (7 = 0, the loop eguation eg. (l7-8p reads: 

{W[''\x)f = V'{x) Wf'\x) - P[''\x) (7-12) 

which implies: 

Wf'\x) = i (v\x) - \Jv'{xy -AP[^\x) ^ (7-13) 
where Pf\x) is a polynomial of degree deg\^' — 1 in x, namely: 

j=2 i=Q 

Notice that eg. (17-31) implies that discrete surfaces with g = Q and n = 1 must have 
f > 2, and thus: 

rP =t5,,, + 0{e). (7-15) 
Therefore P/°^ {x) is a formal series in t such that: 

pf)^:,)=tYl^ + 0{e). (7-16) 

X 

A general spectral curve with filling fractions e, = ^ §^ W^\x)dx, would corre- 
spond to (see eg. (l5-26p ): 

Pr\^) = Y.e, ^ + 0{e) (7-17) 

i=l * 

where Xi, i = 1, . . . ,d — 1 are the zeroes of V'{x). 
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In other words, the spectral curve counting discrete surfaces, has only one non- 
vanishing filling fraction, it is a 1-cut spectral curve, or equivalently, it is a genus g = 
spectral curve. 

More precisely, eq. fl7-16p implies that the zeroes of V'{xf - APi^\x) have the 
following small t behaviour: 

• two zeroes a, b are of the form: 

ar^2Vi+0{t) , b 2Vi + 0{t); (7-18) 

• there are d — 2 double zeroes of the form Xi ±0(t), where V'{Xi) = 0, and Xi ^ 0. 

Therefore there exists a, b and a polynomial M{x) such that: 

V'{xf - 4P/°^ (x) = {x-a){x- b) M{xf (7-19) 

and thus: 

W}°\x) = ^ (v'{x) - M{x)>y{x -a){x -b)^ (7-20) 
where a = 2Vi + 0{t), b = -2^1 + 0{t), M{x) = ^ + 0{t). 

7.3.2 Rational parametrization 

Since the spectral curve has only one cut [a,b], it has genus g = 0, and thus it is a 
rational spectral curve, and it can be parameterized by rational functions of a complex 
variable. Here, this can be done very explicitly. 

We parameterize X as: 

a + b a — b f 1\ ( 1\ a^b a — b 

x{z) = — \ — [z + -\=a + 'y[z + -] , a = , 7 



2 4 \ zj ' \ zj ' 2 ' ' 4 

(7-21) 

This parametrization is convenient because we have: 



^{x-a){x-b) = 7 ( z - i ) (7-22) 



and therefore, from eq. fl7-20p . we see that Wi:^\x{z)) is a rational fraction of z. 

Since x{z) = x{l/z), we can find some complex numbers Uq, Ui, . . . , Ud-i such that: 

d-i 

V'{x{z)) = Y,nk{z'' + z-^) (7-23) 

k=0 

and similarly: 

d-l 

M{x{z)) ^{x{z)-a){x{z)-b) = Y.^kiz'' - z'^). (7-24) 

k=l 

Thus we have: 

W['\x{z)) =uo + \ 5^(n, - n,) ^'^ + ^ ^(m^ + u^) z'K (7-25) 

k>l k>l 
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Since, by definition, Wi'^\x{z)) contains only negative powers of X, Wi^\x) ~ i + 
0(l/a;^), it must contain only negative powers of z, and therefore we must have Uq = 
and Uk = i.e. W^\x{z)) is a polynomial in 1/ z: 



w['\x{z)) = Y,UkZ-\ (7-26) 

k=\ 

Since W^\x) ~ | at large a;, the coefficient oi 1/z must he ui = t/'-). Therefore, a 
and h are determined by the two equations: 

iio = , = (7-27) 

7 

All this can be summarized by the theorem: 

d-l 

Theorem 7.3 Let V'{x) — x — ^^ij^+ix'^ where is the Boltzmann weight for k-gons. 
For an arbitrary a and 7, we write: 



d-l 



V'{a + 7(z + 1/z)) = J2 M^'' + z'^)- (7-28) 

fc=0 

The coefficients are thus polynomials of a and 7. We determine a and 7 hy: 

uq = Q , ui = - (7-29) 

7 

and hy the conditions that a — 0{t) and — t + 0{t^) at small t. 
Then the spectral curve £ — (P^, x, y) is: 

x{z) =a + ^{z + l/z) , y{z) = -Wi°\x{z)) = -J^^kZ^''. (7-30) 



7.3.3 Generating function of the cyhnder, annulus 

For n — 2 and g — 0, Tutte's equations give: 

2Wr\x^)Wr\x„X,) + A Wi'\x^)-Wr\x,) ^ y:^^^^^iO)^^^ ^^ _ p(0)(^^ 

0x2 Xi — X2 

(7-31) 

and one can prove that: 

W^^\x{z^,x{z2)) = 7 , , , - , , , , (7-32) 

' ^ ^ " {Z^-Z2)^X'{ZI)X'{Z2) {X{zi) - X{Z2))^ ^ ' 

i.e. 

Wf + - Bi.,...) ^ ^ (7-33) 

which is the Bergmann kernel on P^, i.e. the Bergmann kernel of the spectral curve 
£ = {¥\x,y). 



93 



7.3.4 Generating functions of discrete surfaces of higher topologies 

For arbitrary n and the generating functions Wn"* counting discrete surfaces of genus 
g with n boundaries, are obtained from the cj^f^'s of the spectral curve £ = (P\x, y), 
by: 



-- Wlf\x{zi), . . .,x{Zn)) dx{zi) . . .dx{Zn) 

dx{zi)dx{z2) , . 

{x{z,) - x{z,)r ^^-^^^ 

In particular, with n = 0, the generating function for counting surfaces of genus g 
and no boundary, is given by the symplectic invariants of S: 

7.4 Example quadrangulations 

If we count only quadrangulations, we choose ^4 7^ 0, and all other tk = 0, i.e.: 

V'{x) = x-Ux^. (7-36) 
The spectral curve is found from theorem eq. fl7.3p . We write: 

x{z) = a + -f(z + -) (7-37) 



z 



and 

V'{x{z)) 

= x{z) — t4^X^{z) 

= a + -f{z + z'^) -ti(^a^ + 3a^-f{z + z"'^) 

+3a-f\z^ + 2 + z^^) + 7=^(^3 + 3z + 3z-' + z^^) ) . 

(7-38) 
In other words: 

2uo = a - U{a^ + Qa-f^) , = 7 - 3t4(a^7 + 7^) (7-39) 

U2 = -3t4«7^ , U3 = -ti-f^. (7-40) 
The condition Uq = implies: 

= a(l -t4(a^ + 67^)). (7-41) 

Since we must choose a solution where a = 0(t) and 7^ = 0(t) at small t, we must 
choose a = 0. Then the condition ui = t/7 gives: 
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I.e. 



2 1 - VI - 12^4 



7^ = ^ (7-43) 



The spectral curve is then: 



x{z)=^(z+-] , y{z) = -± + ^-^. (7-44) 

\ Z J '-)Z Z'^ 

7.4.1 Rooted planar quadrangulations 

The number of planar {g = 0) quadrangulations with only 1 boundary of length 21 is 

T^i^^ = Resx(2;)^' y{z)dx{z) and thus: 

oo 

rr=7^'y^^ (2(/ + l)t-/7^). (7-45) 

In particular, if we require all faces, including the marked face of the boundary, to be 
quadrangles, we choose 21 = 4, and we find the generating function counting planar 
quadrangulations with one marked edge (rooted quadrangulations)!^- 

(0) _ \^-i> \^ +"4(5) 



4 



SeM[°\v),l{S)=4: 



(3t - y ) 

2 3"( 
n! (n + 2)! 

Thus we recover the famous result of Tutte |113j that the number of rooted planar 
quadrangulations with n = + 1 faces is: 

2 3"(2n)! 



n\ in + 2) 



7.4.2 Quadrangulations of the annulus 



(7-47) 



The generating function counting quadrangulations of the annulus = 0, n = 2 is given 
by the Bergmann kernel: 



^h{^)^'2{.) ^Aut(S) 



(7-4^ 



{Zi - Z2y X'{zi)x'{z2) {Xi - X2)'^ 

and thus, if we fix the perimeter lengths of the 2 boundaries as 2li and 2/2, we have: 

,n4{S) 

%f\i = Tf y ' 



S&d-^'^ {v),h {S)=2h ,l2 (5)=2/2 



#Aut(5) 



^^Notice that discrete surfaces with n = 1 marked face and marked edge can have no non-trivial 
symmetry conserving the marked edge, and thus #Aut = 1. 
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= Res Res — -. ^ — dzidz2 

Zl->-O0Z2^0O l^Zi — Z2) 

= 72(...) 2, (2^0K2/.)! 

(7 - 49) 

If we require the marked boundary faces to be quadrangles, we choose 2/i = 2/2 
and thus: 

^n4(5) 

^i? = E*" E 

" 5eM^°\i;),h(5)=4,Z2(5)=4 
= 367' 

= 36 y ^.^r^^l', (3^4)"-' (7-50) 
4^n!(n + 2)!^ ^ ^ ^ 

i.e. the number of annulus quadrangulations with n = + 2 faces, where all faces 
including the 2 marked faces are quadrangles, is 4 x 3" n?(^"+2^)! • 

7.4.3 Quadrangulations on a pair of pants 



'(0) _ Y^+r; ^4 

#Aut(5) 



The generating function for quadrangulations on the pair of pants o = 0, n = 3, is 
ven by : 

(0), . 1 / dzi dz2 dz^ dzi dz2 dz^ 



u (0) 
given by : 

(7-51) 

and for instance we find: 

T^, = (12)3 2!! (^1 + ^ ^ =t''y2' r , ^^"".t^^' ^. {tur-' (7-52) 
^'^'^ ^ ^ 2i V VI - 12^4/ ^ (n + 2)!(n- 1)! ^ ^ ^ ^ 

i.e. the number of quadrangulations on a pair of pants, where all n faces, including the 
3 marked faces, are quadrangles, is 2^ 3" (^n+2)\^{n-iy. ■ 

7.4.4 Quadrangulations on a genus 1 disc 

The generating function for quadrangulations on the genus 1 disc = 1, n = 1, is given 

(1) -z + 8z^-z^ + ^^iz-5z^ + z^) 

= 3(72 - 2)2(.2 _ 1)4 (7-53) 

r« = = ^ ( ^ - ^ ] (7-54) 

^ (7^-2)2 6t4Vl-12t4 VI - 12^4/ ^ ^ 

ri^) = 2 (1 + ( - V2 ) (-1)"-) (12^4)"-^ = 2 (^1 - ^^^^) (12^4)"-^ (7-55) 

i.e. the number of rooted quandrangulations of genus 1 with n faces is: 

n\2-^ n ' 6 V n\n\ ^ ' 



by J?: 



and: 
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7.5 Colored surfaces 



Exactly like the 1-matrix integral is related to the enumeration of discrete surfaces, the 
2-matrix integral and the chain of matrices are also related to enumeration of discrete 
surfaces carrying colors (1 color per matrix). 

7.5.1 The Ising model on a discrete surface 

The Ising model is a problem of enumeration of bicolored discrete surfaces, and it is 
related to the 2-matrix model. It was introduced by Kazakov [8T] . 

Consider M^f^ = the set of connected orientable discrete surfaces of genus g obtained 
by gluing together polygonal faces of two possible colors (or spin) ±, namely 77,3+ 
triangles of color +, 71,4+ quadrangles of color +, ... nk+ fc-angles of color +, and also 
n3_ triangles of color — , n^^ quadrangles of color — , ... Uk- fc- angles of color — , as 
well as n marked polygonal faces of color + of perimeters Zi, . . . , /„, each of the marked 
faces having one marked edge on its boundary. Let us call v the number of vertices of a 
discrete surface, and let M.n\v) be the set of discrete surfaces in Mli^'', with v vertices. 

We call ?T,++ the number of edges with + on both sides, n the number of edges with 

— on both sides, and the number of edges separating faces of different colors. 

We require that unmarked faces have perimeter > 3, whereas marked faces are only 
required to have perimeter k > 1. Moreover, we recall that marked faces are required 
to have color +. 

We define their generating functions as follows: 

the generating function is the formal power series in t: 

Wjf\xi, . . .,Xn;t2,t3, . . . ,td;i2,i3, ■ ■ • Jrf-; t) 

Xi ^ 



^ ^ #Aut(5) 



.n3+{S) ,nd+{S) rn3-{S) Y^d-(S) .n—{S)rn++{S) 

tg • • • f-rf • • • ''d (-2 ^2 



(7-57) 



^^l+'iCS") ^ ^ ^ (^2^2 - l)"++{'S')+n++{5')+n+_{S) ' 



For instance if we want to consider the Ising model on a random triangulation, we 
choose tk = for > 4, and tk = for k > 4. 

One may write Tutte-like recursion equations on the number of edges, which are 
identical to the loop equations of the 2-matrix model {108\ |48] , and for which we may 
use the results of section 15.21 Like in the 1-matrix model, the spectral curve is found 
by the requirement that Wi^{x) be a 1-cut solution of some algebraic equation, and 
have a good small t expansion. The recipe for finding the correct spectral curve is 
summarized in the following theorem: 
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Theorem 7.4 We define the potentials: 

d-l d-1 

v;{x) = t2x-Y, tj+ix^ , v^{y) = hy - Yl ij+^y'- (7-58) 

j=2 j=2 

For arbitrary coefficients ■j,ai,Pi, we define the two rational functions: 

d-l d-1 

x(z) ^-fz + J^O'i^'' ' y(z) ^-fz'^ + (7-59) 

i=0 i=0 

The coefficients 7,Q;j,/9j are uniquely determined by the conditions: 

V;{x{z)) - y{z) ~ - + 0{z-') , Vi{y{z)) - x{z) ~ - + 0{z') (7-60) 

7^ 7 

and such that 7^ as well as ai,Pi are power series in t which behave like 0{t) for t — >• 0. 
Then, the spectral curve is: 

Sluing = {CP\x,y), (7-61) 

the function W^^^ {x) is given by: 

wi'\xiz))^Vlixiz))-yiz), (7-62) 
the function W2^\xi,X2) is given by: 

and all the stable Wn^ 's are given by the ojn^ 's of the spectral curve falsing by: 

Wlf\x{zi), . . . , x{Zn)) dx{zi) . . . dx{Zn) = OO^f^Zu ...,Zn). (7-64) 

In particular, the generating function for counting bicolored maps with no boundaries 
are the symplectic invariants Fg(^^ising)- 

Example: Ising model on quadrangulations. 

We have 

V;{x)^t2X-t4X^ , V^{y)=i2y-Uy\ (7-65) 

We find: 

x{z) = 72; + aiz^^ + a-iZ^^ , y{z) = ^z~^ ii^z + i^zz^ (7-66) 
with the equation: 

/?3 = -^47^ , A = ^27 - 3t47^"i , "3 = -^47^ (7-67) 

ai = ^27 - 3t47^/3i , aiA + Sag/^s = 7^ + ^- (7-68) 
That gives an algebraic equation for 7^: 



r ,4 I (^2-3t2t47^)(t2-3t2^V) t 



and we choose the unique solution such that: 



l'--j^ + 0{t'). (7-70) 
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7.5.2 The chain of matrices discrete surfaces 



A chain of matrices (see section 15.31) . with m matrices Mi, . . . , Mm can also be in- 
terpreted as a generating function for enumerating discrete surfaces with m possible 



colors. The "colors" are labeled 1, 



m. 



We are going to consider discrete surfaces, whose unmarked faces can have any 
color, and are at least triangles, and marked faces have color 1. 

Consider M^f^ be the set of connected orientable discrete surfaces of genus g ob- 
tained by gluing together polygonal faces of m possible colors k = 1, . . . , m (let nj^k be 
the number of faces of size j and color k, we assume j > 3), and n marked faces of 
color 1, with a marked edge, and of size /j > 1, i = 1, . . . , n. 

Let n^i^jy be the number of edges such that the two sides are faces of color i and j. 

Let C be the following Toeplitz matrix of color couplings: 



/ t2,l -10 ... \ 

-1 t2,2 -1 

■•. ■•. ■•. 

\ 0-1 t2,m ) 



(7-71) 



We define their generating functions as follows. The generating function is the 
formal power series in t: 

— ^n,\^Qf) 

oo TT rrm ."j,j(S') -r-rm ^"<i,j>(5') 

h h,,. xl+'^^'^-.x^'"^^^ #Aut(5) ■ 



(7 - 72) 



Again, one may write Tutte-like recursion equations on the number of edges, which 
are identical to the loop equations of the chain of matrices model |5T], and we may 
apply the results of section 15. 3[ Like in the 1-matrix model, the spectral curve is found 
by the requirement that wf{x) be a 1-cut solution of some algebraic equation, and 
have a good small t expansion. The recipe for finding the correct spectral curve is 
summarized in the following theorem: 

Theorem 7.5 We define the 'potentials Vi, . . . ,Vm by: 

Vkix) = t2,kX - ^ tj+i^kx\ k = l,...,m. (7-73) 

i=2 
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For arbitrary coefficients ai^k, we define the rational functions: 

rk 

Xk{z) = ^ aj^kZ^ (7-74) 



with: 

ri = 1, rfc+i = rfc.(4 - 1), (7-75) 

Sm = 1, = Sfc.(4 - 1). (7-76) 
The coefficients 7, ai^k o'^e uniquely determined by the conditions: 

= = 7, (7-77) 

\f k = 2, . . . ,m - 1 , Vl.{xk{z)) = Xk^i{z) + Xk+i{z), (7-78) 

V/c = l,...,m-l, y^J aj^k+ia-j,k = t (7-79) 

i 

and such that 7^ as well as ai^k are power series in t which behave like 0{t) for small 
t. 

Then, the spectral curve is: 

^ch.mat = (CP\X1,X2), (7-80) 

the function W^^'' (x) is given by: 

W^'\x^{z)) = V;ix,{z)) - X2{z), (7-81) 
the function W2^\xi{zi),xi{z2)) is given by: 

Wf{x^{z^),x^{z,))x[{z^)x',{z,) = , ^ - f^^^^^^i^^^^^^ (7-82) 



(2:1-^2)^ ixi{zi) - Xi{z2)) 



and all the stable Wn^ 's are given by the ujn^ 's of the spectral curve i^ch.mat by: 

Wjf\x^{zi), xi(2;„)) dxi{z,) . . . dxi{zn) = u;i^\z^, z^), (7-83) 
i.e. Fg is independent of k. 

Remark 7.1 Because of eq. (j7-78p . we see that the spectral curves {CP^,Xk,Xk+i) are all 
symplectically equivalent for any A; = 1, . . . , m — 1, and thus, we have: 

Vfc = l,...,m-1, Fg = Fg{iCP\xk,Xk+i)) = Fg{{CP\xk+i,Xk)). (7-84) 
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8 Double scaling limits and large maps 



8.1 Minimal models and continuous surfaces 

In the preceding section, we explained how the symplectic invariants can be used to 
count discrete surfaces. The theorem 14.231 allows to find various limits of symplectic 
invariants, and here, it can be used to find the asymptotics of generating functions of 
large discrete surfaces [78] . 

The conjecture |116j was that large discrete surfaces tend towards continuous sur- 
faces weighted with the Liouville theory action, possibly coupled to some conformal 
matter fields. 

The idea is to count discrete surfaces made of large numbers of polygons, and send 
the size of polygons to 0, so that the total area remains finite [291 1^ 169] . 

In section [7] we found the generating function for counting maps of genus g as the 
symplectic invariants Fg{S) of some rational spectral curve S. In this enumeration of 
maps, the expectation value of the number of fc-gons in the considered maps is given 
by 

K)=4^^ (8-1) 
atk 

and the expectation value of the number of vertices is: 

<.)^*^. (8-2) 

Large discrete surfaces are obtained when those numbers diverge, i.e. when the pa- 
rameters tk (or t) approach a singularity, for which the spectral curve is singular (in 
the sense of eq. fl2.3l) ). 

Let us now study the blow up of the matrix models' spectral curves around these 
singularities. 

8.2 Minimal model (p, q) and KP hierarchy 

For discrete surfaces (formal 1-matrix model), or bicolored discrete surfaces with an 
Ising model (2- matrix model), or for the formal chain of matrices, the spectral curve 
depends on the parameter t in an algebraic way. We choose potentials Vk as well as a 
critical t = tc, such that at t = tc the spectral curve has a cusp singularity of the type 
y{z) ~ {x{z) — x{a)Y/'^ (see [31] for a list of critical potentials of minimal degrees, for 
the 2- matrix model). 

We expand 8 in the vicinity of the critical branchpoint z ^ a and t — > tc- 

z = a + {t-tcyC (8-3) 

r x{z) = x(a) + (t - t^r Qic) + o((t - t^r) , . 

I viz) = via) + it- t^r pio + oiit - ur) ^ ' 

where Q and P are polynomials of the complex variable (, of respective degrees q and 
p, and where the exponent u is a scaling exponent such that the blown up spectral 
curve is regular. 
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We define the double scaling limit [78] spectral curve as the blow-up of the singu- 
larity: 

_ r xio = QiO 

In order to find the exponent u, one may notice that for the chain of matrices (and 
thus also 1-matrix and 2-matrix model), the derivative with respect to t of the form 
ydx is: 

±(ydx{z)) = -r^{t-tX-dC. (8-6) 
at z a 

Comparing this with eg. (18-41) . it is easy to see that this implies: 

(8-7) 



p + q - I 
and: 

pPiOQ'iC) - qQiOP'iC) = -■ (8-8) 

ct 

Theorem 14. 231 of section HSl imply that the symplectic invariants defined in eq J2-46l 
give the double scaling limit: 

F,{S) r^it- t^f~'^^^ F,{Sm) (1 + o{l)). (8-9) 

The spectral curve eg. (18-51) . is the spectral curve of the {p, q) minimal model in con- 
formal field theory [2S], coupled to gravity [S3]. It corresponds to a finite dimensional 
irreducible representation of the group of conformal transformations, it has a central 
charge: 

c=l-6^^. (8-10) 
pq 

The exponent z/ = p^g_i or more precisely the exponent "7— string": 

7 = 2-2(p + g)z/ (8-11) 

is given by the famous KPZ formula [39l [83] ■ Notice that the symplectic invariance of 
the Fg's under x y, is related to the (p, q) {q,p) duality [83] . 

The corresponding tau function is: 

r{N) = exp (^(p,.))) • (8-12) 

It is a tau function of the {p, q) reduction of the integrable hierarchy of Kadamtsev- 
Petviashvih (KP) [H]. The function 



00 



m = Y.t^'-''^^^F, = ln(r(t^)) (8-13) 

9=0 
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and its second derivative 

u{t) = F"{t) (8-14) 

can be found as follows: find two differential operators P = dP + pu(P'~'^ + . . . and 
Q = dfl + ud'^~'^ + . . . (where d = d/dt) satisfying the string equation: 

[p,g] = i. (8-15) 

This equation implies the differential equation of KP for the function u{t). 

For example, for pure gravity (p, q) = (3, 2) (with central charge c = 0), we have: 

Q = d'^-2u , P = d^-3ud--u' (8-16) 

and [P, Q] = 1 imply the Painleve equation for u = F": 

3u^ - ^u" = t. (8-17) 

8.3 Minimal model {p, 2) and KdV 

As a special case, we consider the 1-matrix model. In that case, the spectral curve is 
always hyperelliptical y'^ = Pol(x), and thus the only cusp singularities must be half 
integers y ~ x^/^, i.e. q = 2 and p = 2k + 1. 

The operators Q and P are such that: 

Q = rf^-2u(t), (8-18) 

2fc-2 

P = _ (2A; + l)u{t)d^''-^ + J2 Vj{t)d^. (8-19) 

j=0 

The string equation [P,Q] = 1 implies a differential equation for u(t): 

Rk+iiu)=t (8-20) 

where Rk is the k^^ Gelfand-Dikii differential polynomial (see [391 [381 [H] ) . They obey 
the recursion: 

Ro = 2 , = -2u R!j - u' Rj + - R'^'. (8-21) 



The first few of them are: 



Ro ^ 2 
Ri = -2u 

R2 = 3u^ - -u" 

2 

Rs = -5^=^ + ^uu" - ]u'^ - lu" 

2 4 8 

(8-22) 
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We have seen in the previous section, that the spectral curve eg. (18-51) of the {2k + 
1,2) scahng hmit of the 1-matrix model, is: 



x{z) = z — 2uo 

k 



(2fc+l,2) 



j=0 



(8-23) 



It is the classical limit of (P, Q), where d = d/dt becomes a complex number z: 



Q = d^ -2u 



(8 - 24) 



X = z"^ — 2uq 

y = ^2fc+i _ (2A; + 1)moz''=-i + . 



The symplectic invariants defined in eq J2-46l give the double scaling limit: 

(2-2g)(2fc+3) (2-23)(2fc+3) 

^^(^imm) ~ (t - Q (2.+2) Fg{S(^2k+i,2)) + o((t - tc) (^'=+2) ) (8-25) 



and, ii zi, . . . , Zn all lie in the vicinity of 1 + 0((t — tc) ^'=+2 ), we have: 



(2-2g-Ti)(2fc+3) 



(2-2g-n)(2fc+3) 



(8-26) 



where = 1 + (t - tc) 2'=+^ 
The formal function: 



.^-^ (2-2g)(2fc+3) 



(8-27) 



is such that its second derivative 



satisfies the k + 1*^* Gelfand-Dikii equation: 



Rk+i{u) = C 



(8-28) 
(8-29) 



9 Partitions and Plancherel measure 

In many different problems of mathematics or physics, one needs to count partitions 
with the Plancherel Weight. 




h h h 

3 2 1 



104 



Given a partition Ai > A2 . . . > Aat > 0, we define: 

• its weight |A| = Ai + A2 + . . . + Atv- 

• its length n(A) = Aj 7^ 0}, we have n(A) < N. 

• its Plancherel weight: 

where the hiS correspond to the down-right edges of the partition rotated by 7r/4: 

X^-i + N , hi> h2> ...> hN >0. (9-2) 

• The Casimirs: 

^ / TKT ^ /-, \ k 



(9-3) 



Ck{X) = J2{h^-^] -E o-M +(i-2-')C(-fc)- 

j=i ^ / i=i V / 
For example 

1 ^ 

Ci(A) = |A|-- , C2(A) = J]A.(A,-2i + l). (9-4) 



i=l 



9.1 The partition function 

The partition function one would like to compute is: 



n(A)<Af 

and particularly its large Q expansion: 

\nZr,{Q) = J2Q^'^'Fg. (9-6) 
9 

Indeed, since Q is coupled to the weight |A| of partitions, the large Q limit corresponds 
to the limit of large partitions. 

More generally, one is also interested in expectation values of moments of Casimirs, 
and one whishes to compute the partition function: 

d 

Z.{Q,t,)= Y: (=) Q^'^' (9-7) 

where we have taken into account the scaling of the Casimirs in the large Q limit. 
Again, we are interested in the large size expansion: 

lnZN{Q,tk) ~ J2Q'~''^9itk,N). (9-8) 
9 
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In the case where tk = for all k, the answer has been known for a long time, and 
we have the following identity: 



f dim(A)V ^2|A| 



Q^l^l = e'"^ (9-9) 



A 

and therefore: 

F, = 5,,o. (9-10) 

When some of the t^'s are turned on, the answer can be written in terms of a matrix 
model and symplectic invariants of an appropriate spectral curve. 

Some applications of this model include the statistics of the longest increasing 
subsequence of a random sequence, which is equivalent to the statistical physics of 
growing 2D crystals, indeed the Plancherel measure V{X) is precisely the probability 
to obtain shape A by dropping square boxes from the sky at random times and random 
positions, see |103[ 1105] . 

Another application concerns algebraic geometry and topological string theory. 
The partition function Zjv(Q,tfc) is also the generating function for counting ramified 
branched coverings of P^, i.e. the Hurwitz numbers, see [92l [99l llOOl 11021 1103] . 

It has been observed in many works fiSi [771 1105j . that locally, in the large Q limit, 
the Plancherel statistics of partitions, shows many similarities with universal statistical 
laws observed in various limits of matrix models. In fact, it was found in [56], that the 
similarity is much stronger, and in fact Z]^{Q,tk) is a matrix model for all Q (not 
only large Q). In particular this shows that the F^'s are again symplectic invariants. 

9.2 Matrix model for counting partitions 

Consider the contour C which surrounds all positive integer points. 



+-H \ \ \ \ h 

n 1 2 3 4 5 



The following matrix Integra 0: 



Z= / rfMe-'3T'-^(^) (9-11) 

JHn{C) 



can be written in eigenvalues: 



/ ^^i---^^^ \{{^^-^if n^^'^"'^"'^ (9-12) 



where we choose: 



sin (ttQx) r{Qx + 1)2 ■ ^ ' 



Hn{C) stands for the set of normal matrices whose eigenvalues lie on a contour C. 
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The integration over C picks up residues at the poles of e 



-QV(x) 



I.e. 



and thus forces Qxi G N, and the factor Y[ 



.X, 



— Xj)"^ forces x. 



at the poles of 

sinTrQi' """"^ ^""""^ --^-^^^ ^-^i ^ ^n, ^..^ v..^ .^^v^. LLi>j\•^^ -^JJ ^wx^v.^ ^ Xj, thuS We 

have Qxi = /ij G N, and since they are dum variables, we can always reorder them so 
that hi > h2 > ■ ■ ■ > > 0- Therefore we have: 

Z = [ dMe-«^^^(^) 

Jhm{C) 



hi>...>hM>0 
OC ZN{Q,tk), 



Q 



k i 



(9-14) 



i.e. we recover the Plancherel generating function Zi^{Q^tk) for partitions eg. (19-71) . 
Therefore, Z]^{Q,tk) can be written as a normal matrix integral with eigenvalues sup- 
ported on C. Since the loop equations are independent of the integration contour 
(provided there is no boundary term when integrating by parts), we have the same 
loop equations as for any other matrix model, and thus the same solution in terms of 
symplectic invariants. 

The spectral curve, i.e. the equilibrium density of eigenvalues of the matrix, is 
computed like for any other matrix model, like in section O Here, the potential V{M) 
may look quite complicated, and the corresponding spectral curve also looks at first 
sight quite complicated. However, due to the special properties of the T functions, it 
simplifies considerably, and reduces to a rather simple expression, which coincides with 
a "naive" large Q limit. The naive large Q limit has been computed in many works 
[931 11021 11031 1114] . and has been known for some time. It has the property that it is 
nearly independent of N. Let us just state the final result, and refer the reader to [56] 
for details. 

The spectral curve is obtained from the following recipe: 



Theorem 9.1 Define: 



d-l 



U{x)=J2tk, 

k=l 

Then define the coefficients Uq,Ui, . . . ,Ud-i by the equation: 

d~l 

[/(e-"°(z + 1/z - a)) = J^M^'^ + 
Then the spectral curve is: 



(9-15) 



(9-16) 



k=0 



N- 



Plancherel 



XiZ] 



Q 



-he-"o(2 + 1/2 - Ml) 



d-l 



\n.z + \^^Uk{z^ - z ^] 



k=l 



(9-17) 
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and the large Q expansion of Zi\f{Q,tk) is given by: 



In Zm{Q, h) ~ ^^(^Plancherel) (9-18) 

3=0 

where Fg are the symplectic invariants defined in section\^ 

N-- 

Remark 9.1 Since symplectic transformation, Fg is independent of A^, 

and therefore eq. (|9-18|) seems to be independent of N. In fact, the N dependence is smaher 
than any power of Q, it is exponentially small. This phenomenon is known as the arctic 
circle property. All partitions with size n(A) > Q seem to be "frozen" and contribute only 
exponentially to the partition function. 

9.2.1 Example no Casimir 

When all = we have U = 0, i.e. Uq = and Ui = 0. The spectral curve is thus: 



N-- 

x{z) = + z + \/z 

N— - 

y{z)=\nz = Arcosh((x ^)/2) 



x{z) = + e-"o(z + l/z- Ml) 
y{z) = ln2;+ ^Ui{z - i) 

For that spectral curve one finds: 



g-2M0 



2 ' 
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(9-19) 



For that spectral curve one finds: 

Fo = 1, Fi = F2 = F3 = . . . = (9-20) 

which is in agreement with 

ZiQ)=e'^\ (9-21) 

9.2.2 Example: Plancherel measure with the 2nd Casimir 

If we choose ^2 7^ and all other = 0, we have from eq. fl9-16p : 

^26"''° (^z + ^-ui^ = 2uo + ui(^z+ (9-22) 

and thus: 



tt, = -2uoe^"° , Ml = V^2^. (9-23) 

That gives: 

2uo = L{-tl) (9-24) 

where L{x) is the Lambert function, solution of Le^ = x. 
The spectral curve is thus: 



(9-25) 



l + Mo)(2-Mo), (9-26) 



Fi = lln(e-2"«(l + 2no)), (9-27) 

e2"« uf,(l- 12%) 

F2 = ^ 9-28 

180 (1 + 2mo)^ 

and so on... 

9.3 g-deformed partitions 

The g-deformed Plancherel weight is obtained by replacing integer numbers /ij G N by 
the g— numbers [hj\ = ^ , and thus: 

<j5 -g 5 



dim,(A)V_ (U^>^h^~h,]\ _U^>JiQ^-q 



n(A) = (^=F) =i^^F7^i = ::7:.., 0-29) 



which can also be written: 

n(A) = n(^'' - ^"'f n^^""^^'' ( ^ ) (9-30) 

where g{x) is the g-product: 

nfl-ig-y (9-31) 



00 

9[x) 



1 X ^ 
n=l X 



Again, our goal, for various applications in physics and mathematics, is to compute 
the following sum: 

ZN{q;tk) = ^ P,(A) e k (9_32) 

n{X)<N 

and one would like to compute it in the g ^ limit, i.e. 

oo 

lnZ;v~$^(lng)'-'^F,. (9-33) 

3=0 

Again, we will find that the Fg's are the symplectic invariants of some spectral curve. 

The main application concerns algebraic geometry and topological strings. ZN{q] tk) 
is the partition function for the Gromov-Witten invariants of some family of Calabi-Yau 
3-fold, see [92]. 
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9.3.1 Matrix model 

Again, tlie idea is to represent the sum eq. fl9-32p as a matrix integral. 

Consider the contour C which surrounds all points of the form 1, q, g^, g^, . . ., i.e. C 
is a circle of radius 1 < r < |g~^|: 




The following matrix integral: 



Z= / ^iMeT^^^^^^^) (9-34) 

can be written in eigenvalues: 

[ dx^...dxN llix,-x,y He^^""^^^^ (9-35) 



where we choose: 



9(1 



with 



i7ri2^ /^N9 ('"^) 

e ii"? oil) e 2intj 

/(^) = - ^ ; ^ . . ^ • (9-37) 

^ ^ v^(l-x)(7(x)(7(l/x) ^ ^ 

The integration over C picks up residues at the poles of e'^^^^\ i.e. at the poles 
of /(x), and thus forces G N, and the factor ]^(a;j — XjY forces Xi ^ Xj, thus we 

i>j 

have Xi = g'^' where hi G N, and since they are dum variables, we can always reorder 
them so that hi > h2 > ■ ■ ■ > hj^ > 0. The residues are such that we have: 

Z^(g;tfc)oc / dM e^^'^^'''\ (9-38) 

Therefore, as in the preceding case, Z]\f{q, tk) can be written as a normal matrix integral 
with eigenvalues supported on C. Once again, the small q expansion of the free energy 
is thus given by the corresponding symplectic invariants. 
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9.3.2 Spectral curve 

It remains to compute the spectral curve, i.e. the equilibrium density of eigenvalues 
for the matrix potential: 

V{x) = V^(lnx- (A^- -)lng)^' + (lnx)2 + i7rlna; 

k 

ad/x) 

+ (1 - A^)lng Inx + lng ln ^\\^ - In g In (1 - x). (9-39) 

9{x) 

Again, the spectral curve, i.e. the equilibrium density of eigenvalues of the matrix, 
is computed like for any other matrix model, like in section[5l Here, the potential V{M) 
may look quite complicated, and the corresponding spectral curve also looks at first 
sight quite complicated. However, due to the special properties of the g-product g{x), 
it simplifies considerably, and reduces to a rather simple expression, which coincides 
with a "naive" large Ing limit. The naive large Ing limit has been computed in many 
works [92], and has been known for some time. It has the property that it is nearly 
independent of N. Let us just state the final result, and refer the reader to [56] for 
details. 

In principle, the spectral curve could be found for all t^'s, but for simplicity, we 
compute it explicitly for ta = ^4 = . . . = 0, i.e. for only ti and ^2- We write: 

ti = t , t2=p-l , t3 = U = ... = 0. (9-40) 

In that case the spectral curve is: 



£={ " , e-* = - ( 1 - - ) (9-41) 

^(-ln. + |ln(jf^ 



^0 V ^0 



J2 (^^^^) e-'\'\ q"^^'^^'^ ^e^^^'^^^'-^'^^^'K (9-42) 



yy^) x{z) \ "--I 2 ""\^ 1-1/220 

and thus we have: 

The large Ing expansion of \nZ is given by the symplectic invariants of curve S, and 
this expansion is independent of N, provided that N > n, where: 

n=l + \nq (p In (1 - l/^o) + (p - 2) In (1 + l/^o)) • (9-43) 

In fact, this means that the dependence is in non-perturbative terms, smaller 
than any power of In q. This is again the arctic circle phenomenon, the partitions of 
size > n have an exponentially small probability, the size of the system seems to be 
frozen to n. 
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9.3.3 Mirror curve 



In topological strings (see [3l [711 [90l [22l 1115] and section [TTl). it is known that Gromov- 
Witten invariants of some Toric-Calabi-Yau manifolds, can be written as sums over 
partitions, this is called the topological vertex method. 

In particular, for the Calabi-Yau manifold Xp = 0{—p) © 0{p — 2) — > P^, which 
is a rank 2 bundle over P^, the Gromov Witten invariants Mg^diXp) which count the 
number of Riemann surfaces of genus g and degree d which can be embedded in Xp 
and going through given points, are given by: 

J](lng)2-2^ 5^e-*'^Ar,,,(X,) = \nZ (9-44) 

9 d 

where 

A 

Therefore, we have found in the previous section that the Gromov- Witten invariants 
of Xp are given by the symplectic invariants of £: 

Y,^-"'MUX,) = F,{S). (9-46) 

d 

This result is interesting by itself, since it already gives a practical way of computing 
the Gromov- Witten invariants of Xp. 

But we can go further. Notice that x{z) = u{z) is a rational function of z, and 
y{z) = \nv{z) where v{z) is also a rational function of z, and notice that: 

dx A dy = dlnu A dlnv (9-47) 

thus, S is symplectically equivalent to the following spectral curve: 

x(z) =ln f(l- ^)(1- ^ 



1 ( 1-z/zf, \ 2\ ^0 ^0 

1-1/zZQ 



e-* = ^ (1 - 1)^(^-2) (9-48) 



£ = 

y{z) = In 
and thus we have: 

F,{S) = F,{S). (9-49) 

Notice that u = and f = e^ are both rational fractions of z, and thus, by 
eliminating z, there exists an algebraic relationship between them, i.e. there exists a 
polynomial H such that 

H{u,v) = 0. (9-50) 
This curve is known in the context of topological strings [71] : 

H{e^,ey) =u+uj^ (9-51) 

is a Calabi-Yau 3- fold Xp, which is mirror of Xp under mirror symmetry , and if (e^, e^) = 
is the singular locus of Xp. Therefore, we have obtained that, in agreement with the 
"remodelling of the B- model proposal" (see [22] and section [TT]) . we have: 

the Gromov- Witten invariants of Xp, are the symplectic invariants of the 
singular curve of its mirror Xp. 
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10 Intersection numbers and volumes of moduli spaces 

10.1 Kontsevich integral and intersection numbers 
10.1.1 Matrix integral 

Let A be a diagonal matrix A = diag(Ai, . . . , Xjy). 

Kontsevich's integral p5] is the following formal matrix integral (defined as a formal 
power series in large A) 

^Kontsevich (A) = J ciM e"^ ^ (10-1) 

and we consider its topological expansion: 

oo 

In ZKontsevich(A) = A^^'^^ F,. (10-2) 

9=0 

Notice that, upon shifting M — > M — A we have 

^Kontsevich ( A) = e3 ^"^^ dM e L3 J (io_3) 

and thus, this integral is a special case of 1-matrix integral with an external field (see 
section [5^ . which implies that the coefficients Fg are the symplectic invariants of the 
corresponding spectral curve. 



10.1.2 Kontsevich's Spectral curve 

We have seen in section 15.41 that the topological expansion of a matrix integral with 
external field, of the type: 

Z = j rfMe-^Tf^MKA^Af (^Q_4) 

is given by the symplectic invariants of its spectral curve: 

\nZ = ^N^-^3Fg{8) (10-5) 

9 

where the spectral curve £ is characterized by the algebraic equation: 



V'{x) - y = 




(10-6) 



where Pi{x) is a polynomial of x of degree at most degV^", which behaves at large x 
like V'{x)/x. Here we have V{x) = ^, hence V'{x) = x^, and thus Pi{x) = x + Pi{0). 

The fact that Z is to be understood as a formal power series at large A, means that 
we look for a rational spectral curve, and this determines all the polynomials Pi{x). 
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Here we find that the rational spectral curve of the form eq. fll0-6p is: 

S^=S'^['l=\-^J'^) (10-7) 

I y{z) = + ti 

where t^ = -i: Tr i and 

A'^ = A^ + ti. (10-8) 
From now on, for simplicity, we shall assume that ti = 0: 

*i=0 = 1tV^ (10-9) 



and therefore we have 
and the spectral curve is: 



A = A (10-10) 



x(z) = z — ^Ti 



Ek = \ ^ . (10-11) 

I y{z) = z' 

10.1.3 Symplectic invariants 

To compute the F^'s of the spectral curve £k, we need to consider the branchpoints, 
i.e. the zeroes of x'(z), and they are quite complicated. 

However, we may use symplectic invariance, and compute the F^'s after exchanging 
the roles of x and and thus, consider the spectral curve: 

^^=1 yS = f-iTr^— • (10-12) 

t y\Z) Z ^ LL 2A{A-z) 

This spectral curve has now only one branchpoint solution of x'{z) = 0, which is located 
at 2; = 0. Since the F^'s are obtained by computing residues near z = 0,we may Taylor 
expand y{z) near z = 0, and we have: 

I x(^z) = z'^ 

£k — { / \ 1 v^oo , k (10-13) 

I y{z) = Z-lJ2k=0^k+2z' 

Now, it is rather easy to compute the first few symplectic invariants: 



,,(0). x_ 1 dzi dz2 dz3 nni^i^ 

^ ~ ^3 ^1^2 ^3 



U^^2\ZU Z2) = f(2 - h)\^zt + 54 + ^Zfzl) 

ol^ ~ ^3j Z^Z2 I 

+Qtlz^z^ + (2 - t3){6t5Z^zl + Qt^zlzl + 5^72^4) 
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(10-16) 



;io- 17) 



dz 



+/(2 - ^3)^(252 tl + 348 UJrz^ + 145 1?/ + 308 UJqz^) 
+z\2 - t3)(203 h + 145 + 105 z% + 105 /tn) 

+105 (2 - ta) 



And so on... 

For example, the first and second order free energies are: 



Kontsevitch 



24 V 2 



;io-i8) 



and 



^(2) 

Kontsevitch 



1 252 tl + 435 ^5^7(2 - is) + 175 ^9(2 - t-sf 
1920 (2 - tsY 



(10-19) 



Remark 10.1 The fact that the symplectic invariants depend only on odd t^s can be 
understood in terms of symplectic invariance: indeed, adding to y{z) any rational function of 
x{z) (i.e. any rational function of z^) is a symplectic transformation. It does not change the 
Fg^s, and therefore the F^'s depend only on the odd part of y{z), i.e. only on the odd t^s. 



10.1.4 Intersection numbers 

The Kontsevich integral is important because it computes intersection numbers of 
Chern classes of line bundles over the moduli spaces of Riemann surfaces |116l 185] . 

Let Aig^n be the moduli space of Riemann surfaces S of genus g, with n marked 
points pi, ■ ■ ■ ,Pn- This moduli space is a complex manifold of dimension: 

dim A^g,„ = dg,n = 3g-3 + n. (10-20) 

This moduli space can be compactified into a compact space A4g,n by adding all stable 
nodal surfaces (stable means that each component has a Euler characteristics < 0). 

The cotangent bundle £j is the bundle over Mg,n, whose fiber is the cotangent 
space of S at the point pi. Let 

i^i = c,{Ci) (10-21) 

be its first Chern class. Ci and ipi can be extended to Aig^n- 

Chern classes ipi provide useful information on the topology of Aig^n- Indeed if one 
computes the integral of ipi over a cycle in Aig,n, this integral tells how many times 
the cotangent space rotates. The intersection numbers are defined as: 

< r,, . . . r,„ >= / ^|Jf^^Pt • • • ^n" (10-22) 

JMg,n 
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and are non-zero only if ^^rfj = dg^n = Sg — 3 + n. 

i 

The Kontsevich integral is a generating function for those numbers: 

In ^Kontsevich (A) 

= j:(iV/2)-^^ Y: E 11(2^^ - 1)!! ^ < r,, . . . > 

g n di + ...+dn=dg,n i 

(10-23) 

where the sum is restricted to rfj > because we have assumed ti = 0. 
For example we have: 

(10-24) 

10.1.5 Correlators and unmarked faces 

The symplectic invariants of the spectral curve eg. (110-131) are the Fg^s. They generate 
intersection numbers of ■i/'-classes, i.e. Chern classes of cotangent line bundles over 
marked points, i.e. they generate the intersection numbers of the type: 

<^f\..^t">- (10-25) 

The correlators Un^s of that spectral curve, also have some interpretation in terms 
of integrals of some classes over moduli-spaces. 

Notice that Kontsevich integral eq. (ll0-23p contains a summation over n, i.e. over 
the number of marked points. One may whish to distinguish some of those marked 
points, and fix marked faces around them, and perform the sum over the other marked 
points, in some sense forget the other marked points. 

The forgetful map, is the map from Jlig^n+m to Mlg,n, which consists in forgetting 
m marked points. Under this map, ip classes project to the Mumford ^-classes: 

/_ • • • e n = /_ • • • ^n" E n '^(E... a,). (10-26) 

JMg,n + ni = l J Mg,n creS™ C=CycleS of (T 

For examples with m = 1 and m = 2: 



ri'... C" c+i = /_ • • • c" «^a, (10-27) 

^t--- ¥n rn+1 rn+2 = [ ■ ■ ^n" i^a,+a, + /^a./^aj- (10-28) 



Mg,„+2 JMg.n 



One finds [57], that the correlators cun (-^i, . . . , -z^), are the generating functions for 
tv classes, coupled to n T/^-classes: 



[^ly ■ ■ ■ y^n 



2-«^"(t3-2)-^- E E 



do 

(3)1. ,.„, _ 

k\ 

do+di+...+dri=dg^„ k=l bi + ...+bi^=do,bi>0 
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i=i 1=1 1=1 i=i 



(10-29) 

where the coefficients if, are the Schur transform of the i^'s: 

' 2a,- + 1! t2a,+3 



f' = E^^ E n 



/ ^ -"--L a,! ^3-2 

ai+...+a;=6,af>0 j 



Their generating function is obtained by: 



CXD 



2 



;io-3o) 



/W = E^^T-" ' - In (1 -/(.)) = $: 4. ^ (10-31) 

a=l 0=1 

Example: 



ii = -6 , i2 = -60 + 18 , , . . . . (10-32) 

h-2 ' ^3 - 2 (^3 - 2)2 ' ^ ^ 

a;i^^ is the Laplace transform of: 

n Y ^'^J 1 k k 

(io+di+---+tin='i9,n i=l k b\+b2-^—+bk=dQ,bi>0 i=l i=l j 

(10-33) 

Vg^n can be interpreted as the generating function for counting intersection numbers of 
K-classes, on the moduh-space of Riemann surfaces with n boundaries (discs removed 
from the surface), of perimeters Li, . . . , L„. 



10.2 Application: Weil-Petersson volumes 

Consider the stable moduli space M.g,n of Riemann surfaces of genus g with n bound- 
aries (stability means 2 — 2g — n < 0). Every surface in A4g,n has a negative Euler- 
characteristics, and thus a negative average curvature. It can be equipped with a unique 
constant negative curvature metric, called Poincare metric, such that the boundaries 
are geodesies. Let Li, . . . , L„ be the geodesic lengths of the boundaries. 

Every Riemann surface in Adg^n can be decomposed into 2g — 2 + n pairs of pants, 
whose boundaries are geodesies (such a decomposition is not unique). Conversely, 
2g — 2 + n pairs of pants with fixed given boundary perimeters can be glued together 
to form a Riemann surface of Aig,n provided that the lengths of boundaries which are 
glued together match. The gluing is not unique, because 2 circles of the same perimeter 
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can be glued in many ways, twisted by an arbitrary angle. 




The 3g — 3 + n geodesic lengths of the glued boundaries h, . . . , l^gs+n, together 
with the 3g — 3 + n twisting angles 9i, . . . , O-^g^^^n, provide a local set of coordinates 
parameterizing Aig^n- It turns out that, although we don't have uniqueness of the 
decomposition, the corresponding symplectic form, called Weil-Petersson symplectic 
metric is well defined on Air 



<-g,n- 



3(/— 3+n 
i=l 



; 10-34) 



and it can be extended to the compactified M.g^n- The Weil-Petersson volume of M.g^T. 
is then defined as: 



Vol3,„(L] 



, . . . , j^nj 



;i0-35) 



where the n external boundaries are restricted to have fixed geodesic lengths Lj's and 
fixed angles (i.e. some marked points). 

It can be proved |117j that the Weil-Petersson metrics comes from the Kahler 
metrics: 

271^^1 (10-36) 

i.e. we have: 

Vol<;,„(Li, . . . ,L„) = 



(2vrV + ^X:4^^.)^- 



i=l 



) dg^n 



E 



d^] di\ 

dQ+...+dn=Zg—Z+n i=l 

(10-37) 

It can be made to coincide with eq. fll0-33p . provided that we choose ti = 4??^ and 
^3 = 3. Doing the reverse transform using eq. fll0-3ip . it corresponds to — In (1 — f{z)) = 
Air'^z, i.e. f{z) = 1 — e""^"^^ and thus: 

(2z7r)2'^ 



t 



2d+3 



{2d+l)\ 



+ 25dfl 



(10-38) 
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i.e. it corresponds to the spectral curve 



£wp = { i ■ (0 ^ ■ (10-39) 

In other words, the Laplace transforms of the Weil-Petersson volumes Volg^„(Li, . . . , 
are the Wn^^s of the spectral curve Swp- 

Taking the inverse Laplace transform of the recursion relations eg. (12-441) which 
define the wi^^^s for the spectral curve eg. (110-391) . we recover Mirzakhani's recursion 
[971 EH] relations for the Weil-Petersson volumes (see [621 ET] ior the proof). 



10.3 Application: Witten-Kontsevich theorem 

Witten's conjecture [116J was the assertion that the limit of the generating function of 
large discrete surfaces was indeed the generating function of intersection numbers for 
continuous Riemann surfaces. In other words the double scaling limit of the F^'s of 
maps, should coincide with the F^'s of the Kontsevich integral. 

We have seen that the F^'s of maps are given by symplectic invariants, and thus 
their limits as t ^ tc, are given by the symplectic invariants of the blown up curve. 
Thus the F^'s of the double scaling limit of maps, are the F^'s of the minimal {p, 2) 
model, i.e. the F^'s of the spectral curve: 

F,(^^aps) ~ (t - Q^^-^s^^ Fg{S^p,2)) (10-40) 

where /i = and: 

{x{z) = z'^ — 2u 
k=0 

On the other hand, the F^'s of the Kontsevich integral, are also given by the sym- 
plectic invariants of a spectral curve, which is the eguilibrium density of eigenvalues in 
the Kontsevich integral's Matrix Airy function. The spectral curve (see section 110.11) 
is (again we assume for simplicity that ti = 0): 



Sk = \ )_i 2 ^""^""''^ ■ (10-42) 



y{z) = 

We have seen (theorem 14. ip . that if two spectral curve are eguivalent under sym- 
plectic transformations, then they have the same F^'s. In particular the F^'s don't 
change if we change x ^ y and y —x, thus: 

In that new formulation, there is only one branchpoint located at 2; = 0, and since all 
guantities computed are residues at this branchpoint, we may Taylor expand y{z) in 
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the vicinity of z = and thus: 




x{z) = 
y{z) = 



k=0 



oo 




(10-44) 



If we choose the diagonal matrix A such that: 




-fc-2 



(10-45) 



we have: 



(10-46) 



and therefore: 



(10-47) 



which proves Witten's conjecture. 

In fact this conjecture was first proved by M. Kontsevich [85]. Kontsevich's method 
consisted of two parts: first prove that the Airy matrix integral now known as Kontse- 
vich integral was the generating function of intersection numbers, and then prove that 
both Fg K and -Fg,{p,2) obeyed the same set of differential equations, namely KdV hierar- 
chy. Witten's conjecture has received several proof since then, in particular Loijeenga's 
[89], or Okounkov-Pandaripande |104j . 

Here we see that the spectral curve of the Kontsevich model and the spectral curve 
of the {p, 2) model are obtained from one another by the symplectic transformation 
X ^ y,y ^ —X. That symplectic transformation (the tt/2 rotation in the (x, y) plane), 
leaves Fg invariant, but it changes the cul^^'s. 

11 Application: topological strings 

The counting of surfaces with particular weights, or enumerative geometry, has been 
very important in physics since the arising of string theories. Indeed, those theories 
consist in replacing the point-like particles of usual theories such as classical mechanics 
(0 dimensional objects) by strings (1 dimensional objects): a state is now given by a 
string state instead of a point state. A string evolving in time describes a surface in 
the space-time, the world-sheet, instead of a line for a point evolving in space-time. 
Hence the usual path integral corresponding to sum over all possible stories from one 
initial state to one final state is now a "sum" over all possible surfaces in space-time 
(target space), linking the initial strings to the final ones. This could explain why 
these theories are related to the symplectic invariants since the latter already appear 
in many problems of "surface enumeration". So far, no proof is available but there 
exist many hints that the symplectic invariants of some specific spectral curve should 
be the partition functions of some particular string theories: type IIB toplogical strings 
on some special target space. Many checks have been made that this is indeed the case 
[92], so that Bouchard, Klemm, Marino and Pasquetti [22] proposed to define the 
topological string partition function and observables as the symplectic invariants and 
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correlation functions, considering the spectral curve as the target space of the string 
theory. 

Even if this conjecture is not proved yet, some new clues were recently given by 
Dijkgraaf and Vafa [H] who already conjectured that some matrix models are dual to 
some particular topological string theories [3D]. Those new hints rely on the study of 
an effective field theory for the string theory: the Kodaira-Spencer theory. We review 
those advances in the second part of the present section. 

11.1 Topological string theory 
11.1.1 Introduction 

Type IIB topological string theories are obtained by twisting N = 2 superconformal 
sigma model in 2 dimensions. The precise description would lead us too far away 
from the main topic of this review, and the interested reader is invited to consult 
[90| [9T| 11151 [71] for details as well as [23] for the particular topic of toric geometries. 
There exist two ways of twisting this theory leading to two different models referred 
to as A and B models. We will be particularly interested in this paper in the B-model 
and won't develop the A-model too far. Nevertheless, since the special geometries of 
the target space of the B-model we are interested in, are inherited from the A-model, 
we will say a few words about the latter and especially the possible geometries of its 
target space and of the objects one wants to compute in the next section. 

First of all, let us mention that the topological string theories can be thought of, 
as theories of the maps from a Riemann surface S (the world-sheet) to a Calabi-Yau 
manifold M (the target space) of arbitrary dimension. More precisely, in the A model 
side, the amplitudes to be computed are related to the Gromov-Witten invariants as 
follows. Let us consider a worldsheet S^^^ of genus g with k holes (or boundaries). One 
wants to count maps which map the boundaries to a Lagrangian submanifold of M 
denoted as the branj^ L. Such maps are characterized by two additional parameters 
prescribing how the boundaries are mapped to the brane L: a bulk class /? G H2{M, L)0 
and winding numbers Ui E 1^, i = 1, . . . , k, telling how many times the boundaries wrap 
around the brane L. One sums up those "number of maps" called Gromov-Witten 
invariants Njf^ in generating functions: 

/3eH2(A/,L) 

and one also considers the open string amplitudes 

k 

Af\z^,...,z,):=Y,F^^^l[zr (11-2) 

where the open string parameters Zi are parameters of the moduli space of the brane 
L as well as the closed string amplitudes: 

T^s^ := f!^'\ (11-3) 

^^The geometry of such submanifolds is studied in the next section. 
^^One assumes for simphcity that 6i(L) = 1. 
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Let us also precise the moduli spaces of both models, since they play an important 
role in the link between topological string theories and the symplectic invariants. In 
the A model side, the moduli are the Kalher parameters of the target space M whereas 
the B model depends on the complex structure of M. 

11.1.2 Mirror symmetry, branes and toric geometries 

One of the most fascinating features of topological string theories is the existence of a 
duality linking the A and B models: the mirror symmetry. This symmetry states the 
equivalence of the A model on a target space M and the B model on a mirror target 
space M obtained from M by a mirror map which exchanges the Kalher structure of 
M with the complex structure of 

In the following, we will only be concerned with a special but interesting class of 
target spaces M: Toric Calabi-Yau threefolds for the A model and their image under 
mirror symmetry. We now precise the structure of those geometries as well as the 
geometry of their Lagrangian submanifolds. 

Let us start with a Toric Calabi-Yau on the A model side. A Toric Calabi-Yau 
threefold M can be built as a submanifold of C'^"'"^ as follows. Consider + 3 complex 
scalars Xi = e*^% i = 1, . . . , k + 3, transforming under the action of U{1)^ as 

Xi e'^'^'-Xi (11-4) 

for some integers Q", a = 1, . . . , fc. One then considers the 3-dimensional submanifold 
of M^"^^ obtained by constraining the |Xjps to satisfy 

Y,Qt\X^?=ra , « = l,...,fc. (11-5) 

i=l 

The CY 3- fold M is the bundle of Tori generated by the ^j's modulo the action of 
U{lY, over this real submanifold. 

The parameters are the Kahler moduli of the Toric threefold M. The Calabi-Yau 
condition is then: 

fc+3 

Va = l,...,fc , 5^Q? = 0. (11-6) 

i=l 

Remark that one can see the coordinates X^ as a Si-fibration (coordinates e*^') over 
(coordinates \Xi\) giving to M a structure of fibration over the subspace of 
defined by the constraints eg. (11 1-51) . It is then interesting to note that this fibration 
has singular loci when one or several \Xi\ vanish. Indeed, the fiber defined by the 
corresponding 9i shrinks. These loci will be important in the folowing study of the 
type A branes and they are encoded in the so-called toric graph of the threefold M 
i21E2l[23]. 

^''The works on the subject of mirror symmetry are numerous in physics and mathematics. Entering 
this subject would quickly lead us too far away from our main topic. For a nice review of the subject, 
one can refer to the excellent book [7Tj . 
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The boundaries of the worldsheet must be mapped to special Lagrangian subman- 
ifolds of X called branes. For this purpose, let us remind the notation 

X, = \X,\e''\ (11-7) 

The canonical symplectic form on M is then given by 



3 

i , 



1 ^ 

-Y,d\X,? ^de,. (11-8) 



2 



One can cancel this form and obtain a special Lagrangian submanifold L by fixing the 
^j's with the equation 



as well as constraining the moduli of Xj by 

fc+3 



|^^l' = c" (11-10) 



4 = 1 

for a = 1, . . . , r, where the satisfy 

A:+3 



E^" = o- (11-11) 



A special submanifold in M is then given by the set of complex numbers qf and c". 

Moreover, one can consider such manifolds L passing through some singular locus 
of the manifold M. In this case, it splits into two submanifolds and . This is 
one of the latter submanifolds that we consider as A brane, e.g. L^^\ 

Let us now describe the mirror geometry of this target space and of the branes in 
the B model. ^ ^ 

The mirror map transforming M into M can be built as follows. M has homogenous 
coordinates Xi := e^' G C* with i = 1, . . . , A; + 3 whose moduli are constrained by 



Xi 



The mirror geometry M of M is then given by 

fc+3 

Y^X, (11-13) 



i=l 



for two complex scalars {uj^^uj ) G and non vanishing complex homogenous coor- 
dinates Xi G C* satisfying 



fc+3 



J]xf'" = e-*" = g„ (11-14) 



i=l 



One can chose or L as this brane without changing anything in the following. 
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for any a = 1, . . . , /c, where 

ta:=ra + iOa (11-15) 

are complexified Kahler parameters of the threefold M. The equation of the mirror 
geometry M reduces to 

H{X, Y\ta) = uj+uj- = i/(e^, e^|t«) (11-16) 
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where X = and Y = are two non- vanishing coordinates chosen among the X, 
The holomorphic volume form on M is then given by 

Q = ^r^c— = — dxdy. (11-17) 

uXY ^ 

Under this mirror map, one can easily characterize the B model branes, i.e. the 
image of the special Lagrangian submanifolds of M under the mirror map [1]. The 
constraints eq. (lll-lUp are translated into constraints on the Xf. 

fc+3 



l[xf=e-'^- (11-18) 



i=l 

for a = 1, . . . ,r. Moreover, for r = 2, if one considers the singular brane L^, its image 
is a one dimensional complex submanifold described by the algebraic equation 

H{e^, e^) = = H{X, Y) (11-19) 

on C* and can thus be obtained by fixing = and considering uj~^ as a parameter 
of this brane. In the following, we will consider this equation as the spectral curve. 



11.1.3 Embedding of the spectral curve and open string parameters 

In the preceding section, we showed that the spectral curve corresponds to the moduli 
space of the open string boundaries, i.e. of the B-branes. In particular, it was shown 
that the description of this moduli space depends on a choice of local coordinates. Let 
us make this statement more precise by studying the whole target space of the B model 
and its projection to local patches. 

Remember that the moduli space of the B-branes is a Riemann surface £ given by 
a set of coordinates Xj and uj± constrained by 

fc+3 

= cj_ = (11-20) 

1=1 

as well as the constraints 

fe+3 

Y[xf=e-'" , a = l,...,k (11-21) 

i=l 

^^The choice of such coordinates Xi and Xj depends on the sector of the moduh space that we are 
studying. We exphcit this choice in the next section. 
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and 

k+3 

Y[xf=e-''- , a = l,...,r. (11-22) 

i=l 

In order to describe this moduli space, one has to choose as set of coordinates char- 
acterizing it. Indeed, one can describe it as an embedding of the Riemann surface C 
into C* X C* (resp. to C x C) through the coordinates Xi (resp. the coordinates x. 

e.g. the choice of two coordinates Xi and Xj (resp. Xi and Xj) among the set 

(resp. {xa}aZi) allowing to describe C by an equation 



k+3 



a=l 



H,,j{X„X,) = i7,,,(e"%e"0 = (11-23) 

obtained by elimination of all the other coordinates Xa in eg. (111-201) . eq. (lll-2ip and 
eg. (11 1-221) . In the following, one generically denotes this embedding of the spectral 
curve by the equation 

i7(e^,e^) =0 = i/(X,F). (11-24) 

If all these equations represent the embedding of the same surface, they correspond 
to different description of the branes, i.e. different types of boundary conditions for 
the worldsheet. The choice of such an embedding is not random: depending on the 
regime we are considering, i.e. the sector of the moduli space we are studying, some 
embeddings are more appropriate (see for example the discussion in section 2.2 of [22j). 
Remark that the reparameterization group of the spectral curve C is 

Gc = SL{2,Z) X ( J J ) (11-25) 

in terms of the variables {X,Y) of H{X,Y) = 0. 



Remark 11.1 Note that these transformations preserve the symplectic form \dxAdy\ = 
. This reminds of the symplectic transformations section 14.21 which do not change 
ectic invariants. In this topological string setup, these transformations acting on 



dX A dY 
X Y 

the symp^ 



the open string parameters should preserve the closed string amplitudes. 

These transformations, i.e. changing the open string parameters, are important 
in the study of the open string amplitudes. Indeed, the whole open string moduli 
space, or moduli space of branes, exhibits different phases. In each of these particular 
regimes, one can use a specific embedding to describe the brane moduli space in ap- 
propriate coordinates (see [TTl |22] for a review on the subject). The usual methods of 
computation allow to know the open string amplitudes in some very particujar regime. 
This means that one can compute these amplitudes in a specific patch of M and not 
on the others. It is thus interesting to be able to go from one patch_to the others. 
These "phase transitions", corresponding for example to blow ups of M, are elements 
of Gc encoding the transition from the open string parameters of one embedding to 
the others. 
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The choice of an embedding, i.e. the choice of a coordinate x, does not only fix 
the location of a brane but also the last remaining ambiguity known as the framing. 
Roughly speaking, the framing consists in a discrete ambiguity and can be fixed by 
choosing an integer /. This ambiguity correspond to the elements of Gc'- 

(11-26) 

or, in X and y coordinates, 

{x,y)^{x + fy,y) (11-27) 
for integer /. Note that this is also a symplectic transformation considered in section 

El 

We have thus shown that fixing an embedding of the spectral curve £ in C x C (or 
C* X C*), one fixes the open string parameter space which can be seen as the coordinate 
X. Going from one patch in the parameter space to another is obtained by changing 
embedding thanks to an element of Gc- 



11.1.4 Open/closed fiat coordinates 

As it was already mentioned in the preceding sections, the moduli space of the B model 
(resp. A model) is given by the complex (resp. Kahler) parameters of the target space 
M (resp. M). Let us precise the fiat coordinates describing these complex and Kahler 
structures which are mapped to each other by the closed string mirror map. 

These flat coordinates T", a = 1, . . . ,g are given by the periods of the meromorphic 
one form ydx = InY^ on the spectral curve C: 



T"" = ^ (p ydx (11-28) 

'Ac. 



2m 



where {Aa, Ba) is a canonical basis of one cycles on C. This also ensures the existence 
of a holomorphic function F{Ta) such that the dual periods can be expressed as 

What about the open flat coordinates? If the closed coordinates are given by closed 
integrals of A over the cycles, the open flat coordinate is expected to be given by chain 
integrals 

U = — ydx (11-30) 



2m _ _ 

where is an open path over which y jumps by 2iTT. 

Moreover, it is interesting to note that the open string disc amplitude Af'\u) can 
be computed explicitly: it is also a chain integral of the one form 

aP{x) = [ ydx. (11-31) 



23 



See [22] for explanations on this phenomenon. 
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Remark 11.2 Once again, it is interesting to note the similarities between the theory of 
symplectic invariants and B model. Let us summarize this correspondence in a short array: 



Symplectic invariants 


B model 


Spectral curve 


Brane moduli space 


Symplectic transformations 


phase transition 


filling fraction ej 


closed flat coordinates Tq, 


genus zero free energy F^^^ 


prepotential F 


variations - f^. ydx 


variations §^ = ^ 


Genus one point function ydx 


Disc amplitude / ydx 



11.1.5 Symplectic invariants formalism: a conjecture 



In [22j, following some checks of |92] and the seminal paper [40j, Bouchard, Klemm, 
Marifio and Pasquetti proposed to define the open and closed string amplitudes of 
the A model as the symplectic invariants and correlation functions computed on the 
spectral curve of the mirror B model branes. 

The conjecture of [22j simply states that the open string amplitudes Af'\zi, . . . , Zk) 
and closed string amplitudes J-'^^^ of the A model whose mirror background gives rise 
to the spectral curve H{e^,e^) = are given by correlation functions and symplectic 
invariants built from the equation if(e^,e^) = 0: 




;il-32) 



and 



Fh{tc 



(11-33) 



Remark 11.3 In [22], the authors seem to slightly change the recursive rules defining the 
correlation functions and symplectic invariants. However, this apparent transformation re- 
sults from their choice to work with the coordinates X = and Y = instead of x and y 
in order to start from an equation 

H{X,Y) = Q (11-34) 

which is algebraic. It appears that the algebraicity of the spectral curve is not essential and 
that one can directly work with the coordinates x and y, avoiding this change of coordinates. 



Let us emphasize a few important points. In order to get the A model amplitudes, 
one first has to compute the correlation functions and symplectic invariants from the 
B model spectral curve and then plug in the mirror map to obtain the result in terms 
of the A model parameters. It should be underlined that the choice of coordinates 
X and y out of the choice of embedding of £ in C x C, corresponds to 

a choice of brane in the A modeQ Thus, it is interesting to study the behavior of 

^''That is to say, the choice of the location of the Brane as well as a choice of framing. This topic 
is well developed in [22j . 
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the amplitudes when one moves in the brane moduh space, which corresponds to the 
space of parameterizations of the spectral curve. In particular, the closed amplitudes 
should not depend on this choice of parametrization since only the boundaries of the 
worldsheet are sensitive to the definition of the branes. This independence of the closed 
amplitudes on the embedding of the spectral curve is indeed true and follows directly 
from the symplectic invariance of the F^s)^s. 

11.1.6 Checks of the conjecture 

There have been many checks of this conjecture before the definition of the symplectic 
invariants and after. 

First of all, most of this conjecture was inspired by the idea of Dijkgraaf and Vafa 
who conjectured that the partition function of the type B topological string on some 
special backgrounds is given by a random matrix integral [IQ]. Now that the symplectic 
invariants extend the notion of random matrix partition function to any spectral curve, 
it seemed natural to conjecture that these symplectic invariants do coincide with the 
partition function of B model topological string with more general backgrounds. 

Another further property of the symplectic invariants points in the same direction. 
As it is reminded in section 14.4.21 using the modular properties of the symplectic in- 
variants, it was proved that one can promote the latter to modular invariants whose 
non-holomorphic part is fixed by the same set of equations as the B model topologi- 
cal string partition function [63]: the holomorphic anomaly equations of BCOV [17] . 
This means that the non-holomorphic part of these functions coincide. To prove the 
conjecture, one thus have to prove it only for the holomorphic part. 

Further studies were led by Marino and collaborators [921 122]: they checked for many 
explicit examples of possible backgrounds for the B model topological strings that the 
partition function and open string amplitudes indeed coincide with the correlation 
functions and symplectic invariants computed on the associated spectral curve. Every 
single check indeed works, giving more weight to this conjecture! 

Another general check can be made by looking at the short summary made in 
the array in section 111.1.41 It is also interesting to note that the disc and cylinder 
amphtudes can be computed independently from this conjecture for any background 
in the B model: they satisfy the relation conjectured by |22j . 

Moreover, the computation of the sum over large partition with respect to the q- 
deformed Plancherel measure makes the link with the topological vertex approach and 
proves the conjecture in a particular family of target spaces. The extension of this 
method could lead to a direct proof of this conjecture. 

Finally, a last clue has been added recently by Dijkgraaf and Vafa [H], using an 
effective field theory conjectured to be equivalent to B model: the Kodaira-Spencer 
theory. This check is the subject of the next section. 
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11.2 Kodaira-Spencer theory 

11.2.1 Introduction: an effective field theory for the B-model 

The six dimensional Kodaira-Spence theory is the string field theory for the B model 
on Calabi-Yau threefold. Consider the case of non-compact Calabi-Yau threefold M 
defined by 

H{x,y) = uj^uj- (11-35) 

whose holomorphic volume form is 

n = ^dxdy. (11-36) 

The Kodaira-Spencer theory^ is the quantization of the cohomologically trivial varia- 
tions of the operator 5 on M with fixed complex structure. 

The setup dual to the preceding section corresponds to the local surface 

H{x,y) = 0. (11-37) 

Since one can see that the periods of Q can be reduced on this local surface to the 
integrals of the one form ydx, the two dimensional reduction of the Kodaira-Spencer 
theory is defined by the pair [d, ydx) on the spectral curve: it means that it is the 
study of the deformations of d keeping the cohomology class of ydx fixed. For this 
purpose, one looks at the variations of the this operator under the form 

d^d--^d (11-38) 
ydx 

where is a scalar field satisfying 

9(90 = 0. (11-39) 

Indeed, under this variation, the cohomology class of ydx is not changed since it trans- 
forms as 

ydx — > ydx + dcj). (11-40) 
Finally, we are left with a field theory on the spectral curve C given by the actioij^: 

S= f dcfdcj) + + ^90 mf (11-41) 

Jc ^ ydx 

where we rescale the differential ydx by the string coupling constant A = ydx — 
Let us just explain the three different terms of this action. The first term is a simple 
kinetic term whereas the second term corresponds to the coupling to a holomorphic 
background gauge field The most interesting term is the third one: this cubic 
interaction encodes the perturbative corrections and is the fundamental ingredient of 
this action. 

Let us now move to the observables of this theory. First of all, the partition function 
can be written as 

Z = e-^ (11-42) 
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For a very pedagogic construction of this action, see [iT] . 
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where the free energy T has a topological expansion in terms of the string coupling 
constant 

^= ^A^f-^F^^). (11-43) 

One also defines the correlation functions 

W^{z^, . . . , ^fc|A) := (90(zi) . . . 90(zfe))^ (11-44) 

where the subscript c denotes the connected part. These correlation function also have 
a topological expansion 

l^fc(^i, . . . , z^\\) = J2 X''^''-'wt\z^, ...,z,) (11-45) 

coming from the mteraction term e •'^ j"*^ . 

11.2.2 Recursive relations as Schwinger-Dyson equations 

One can remark that the integrant in this interaction can be written as a total derivative 
and does not give any contribution except at the zeroes at the denominator, i.e. the 
zeroes of ydx. These zeroes are the zeroes of y{z) and dx{z). However, one can show 
that only the zeroes of dx do give no n- vanishing contributions (see [H]). Thus, the 
interaction term can be written as follows: 

^'^(^'^)' =yl (11.46) 

ydx ^ J^^ ydx 

This means that the interaction vertex is localized at the branch points. 

In order to compute the correlation functions, one thus has to compute terms of the 

form (^d(f){zi) '^^^]j^z)dx(z)^^'' . . A first step consists in the computation of the 
two point chiral operator {d(f){z)d(f){zi)) which is known to be the Bergmann kernel. 
From this point, one can easily compute the contraction of d(j){zi) with 0(z) 



i{z')=-(i{z) 



{<P{z)d<P{z,)),^^^, = - I Biz',z,) (11-47) 



where the subscript twist refers to the fact that one constrains the scalar field to be 
an odd function of a local variable ^(z) as z approaches a branch point: 

H-a^)) = -m^))- (11-48) 

In other terms, using the notations of section [2], it implies, thanks to De rHopital's 
rule: 

hm = lini I ^ /^f ;) ^ ^ = - hm K{z,, z). (11-49) 

z^ai ydx[z) z^ai 2 [y[z) — y[z))dx[z) z^a^ 
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Finally, taking into account the normal ordering of the cubic interaction term, the 
Schwinger-Dyson equations of this theory give the recursion relation for the correlation 
functions 

z L 
i h=0 IgJ 

(11-50) 

The correlators of the Kodaira-Spencer theory on C are thus the correlation functions 
computed from the latter curve. 

What about the partition function? On the one hand, from the topological expan- 
sion, one easily finds that 

>^^ = J2(^9-2)X''~'F('^- (11-51) 

On the other hand, let us reexpress the LHS directly in terms of the correlators of the 
theory. Indeed, thanks to the expression of the action eg. (11 1-411) . one gets 

In order to compute these terms, one proceed as in the case of the correlation functions 
by localizing these expressions around the branch points. One can first remark that 
the second term vanishes since it corresponds to the interaction operator with no field 
inserted. Let us thus compute the first term. Since the integrant can be written as a 
total derivative ydx dcp = d {ydx 0), this term can be localized around the poles of the 
integrant which are nothing but the branch points, i.e. 

^i = lE(i/''-*)- (11-53) 

Consider now a primitive $ of ydx 



d^ = ydx. (11-54) 



Integrating by parts, it implies 

A 



In terms of the topological expansion, this equation coincides with the definition of the 
symplectic invariants 

F^<^^ = J2 / ^iz)wf'\z). (11-56) 

2 — 2(7 ■ Jai 

This means that the partition function of the Kodaira-Spencer theory is the tau func- 
tion Tjv defined from the symplectic invariants in section [2J 
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12 Conclusion 



In this review article, we have presented overview of the recent method introduced 
for solving matrix models loop equations, and its further extension to a more general 
context. We have defined the notion of symplectic invariants of a spectral curve, and 
we have studied its main applications, as in the present state of the art. In some sense, 
starting from the spectral curve of a classical integrable system, we have proposed a 
way to reconstruct the full quantum integrable system. 

The study of applications to enumerative geometry and integrable systems of those 
notions is probably only at its beginning, and in particular the consequences for topo- 
logical string theory are still mostly to be understood... 
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